<p>This thesis concerns the approximation of optimally controlled partial differential equations for applications in optimal design and reconstruction. Such optimal control problems are often ill-posed and need to be regularized to obtain good approximations. We here use the theory of the corresponding Hamilton-Jacobi-Bellman equations to construct regularizations and derive error estimates for optimal design problems. The constructed Pontryagin method is a simple and general method where the first, analytical, step is to regularize the Hamiltonian. Next its stationary Hamiltonian system, a nonlinear partial differential equation, is computed efficiently with the Newton method using a sparse Jacobian. An error estimate for the difference between exact and approximate objective functions is derived, depending only on the difference of the Hamiltonian and its finite dimensional regularization along the solution path and its<em> L</em>² projection, i.e. not on the difference of the exact and approximate solutions to the Hamiltonian systems. In the thesis we present solutions to applications such as optimal design and reconstruction of conducting materials and elastic structures.</p>

<p>This thesis concerns the numerical simulation of creeping rays and their contribution to high frequency scattering problems.</p><p>Creeping rays are a type of diffracted rays which are generated at the shadow line of the scatterer and propagate along geodesic paths on the scatterer surface. On a perfectly conducting convex body, they attenuate along their propagation path by tangentially shedding diffracted rays and losing energy. On a concave scatterer, they propagate on the surface and importantly, in the absence of dissipation, experience no attenuation. The study of creeping rays is important in many high frequency problems, such as design of sophisticated and conformal antennas, antenna coupling problems, radar cross section (RCS) computations and control of scattering properties of metallic structures coated with dielectric materials.</p><p>First, assuming the scatterer surface can be represented by a single parameterization, we propose a new Eulerian formulation for the ray propagation problem by deriving a set of <i>escape </i>partial differential equations in a three-dimensional phase space. The equations are solved on a fixed computational grid using a version of fast marching algorithm. The solution to the equations contain information about all possible creeping rays. This information includes the phase and amplitude of the ray field, which are extracted by a fast post-processing. The advantage of this formulation over the standard Eulerian formulation is that we can compute multivalued solutions corresponding to crossing rays. Moreover, we are able to control the accuracy everywhere on the scatterer surface and suppress the problems with the traditional Lagrangian formulation. To compute all possible creeping rays corresponding to all shadow lines, the algorithm is of computational order O(<i>N</i>³ log <i>N</i>), with<i> N</i>³ being the total number of grid points in the computational phase space domain. This is expensive for computing the wave field for only one shadow line, but if the solutions are sought for many shadow lines (for many illumination angles), the phase space method is more efficient than the standard methods such as ray tracing and methods based on the eikonal equation.</p><p>Next, we present a modification of the single-patch phase space method to a multiple-patch scheme in order to handle realistic problems containing scatterers with complicated geometries. In such problems, the surface is split into multiple patches where each patch has a well-defined parameterization. The escape equations are solved in each patch, individually. The creeping rays on the scatterer are then computed by connecting all individual solutions through a fast post-processing.</p><p>We consider an application to mono-static radar cross section problems where creeping rays from all illumination angles must be computed. The numerical results of the fast phase space method are presented.</p>

<p>Ill-posed sets of linear equations typically arise when discretizing certain types of integral transforms. A well known example is image reconstruction, which can be modelled using the Radon transform. After expanding the solution into a finite series of basis functions a large, sparse and ill-conditioned linear system arises. We consider the solution of such systems. In particular we study a new class of iteration methods named DROP (for Diagonal Relaxed Orthogonal Projections) constructed for solving both linear equations and linear inequalities. This class can also be viewed, when applied to linear equations, as a generalized Landweber iteration. The method is compared with other iteration methods using test data from a medical application and from electron microscopy. Our theoretical analysis include convergence proofs of the fully-simultaneous DROP algorithm for linear equations without consistency assumptions, and of block-iterative algorithms both for linear equations and linear inequalities, for the consistent case.</p><p>When applying an iterative solver to an ill-posed set of linear equations the error typically initially decreases but after some iterations (depending on the amount of noise in the data, and the degree of ill-posedness) it starts to increase. This phenomena is called semi-convergence. It is therefore vital to find good stopping rules for the iteration.</p><p>We describe a class of stopping rules for Landweber type iterations for solving linear inverse problems. The class includes, e.g., the well known discrepancy principle, and also the monotone error rule. We also unify the error analysis of these two methods. The stopping rules depend critically on a certain parameter whose value needs to be specified. A training procedure is therefore introduced for securing robustness. The advantages of using trained rules are demonstrated on examples taken from image reconstruction from projections.</p> / <p>Vi betraktar lösning av sådana linjära ekvationssystem som uppkommer vid diskretisering av inversa problem. Dessa problem karakteriseras av att den sökta informationen inte direkt kan mätas. Ett välkänt exempel utgör datortomografi. Där mäts hur mycket strålning som passerar genom ett föremål som belyses av en strålningskälla vilken intar olika vinklar i förhållande till objektet. Syftet är förstås att generera bilder av föremålets inre (i medicinska tillämpngar av det inre av kroppen). Vi studerar en klass av iterativa lösningsmetoder för lösning av ekvationssystemen. Metoderna tillämpas på testdata från bildrekonstruktion och jämförs med andra föreslagna iterationsmetoder. Vi gör även en konvergensanalys för olika val av metod-parametrar.</p><p>När man använder en iterativ metod startar man med en begynnelse approximation som sedan gradvis förbättras. Emellertid är inversa problem känsliga även för relativt små fel i uppmätta data. Detta visar sig i att iterationerna först förbättras för att senare försämras. Detta fenomen, s.k. ’semi-convergence’ är väl känt och förklarat. Emellertid innebär detta att det är viktigt att konstruera goda stoppregler. Om man avbryter iterationen för tidigt fås dålig upplösning och om den avbryts för sent fås en oskarp och brusig bild.</p><p>I avhandligen studeras en klass av stoppregler. Dessa analyseras teoretiskt och testas på mätdata. Speciellt föreslås en inlärningsförfarande där stoppregeln presenteras med data där det korrekra värdet på stopp-indexet är känt. Dessa data används för att bestämma en viktig parameter i regeln. Sedan används regeln för nya okända data. En sådan tränad stoppregel visar sig fungera väl på testdata från bildrekonstruktionsområdet.</p>

<p>In this masters’ thesis report we present a graph theoretical method used for automatic bilingual lexicon acquisition with parallel texts. We analyze the concept of graph similarity and give an interpretation, of the parallel texts, connected to the vector space model. We represent the parallel texts by a directed, tripartite graph and from here use the corresponding adjacency matrix, A, to compute the similarity of the graph. By solving the eigenvalue problem ρS = ASAT + ATSA we obtain the self-similarity matrix S and the Perron root ρ. A rank k approximation of the self-similarity matrix is computed by implementations of the singular value decomposition and the non-negative matrix factorization algorithm GD-CLS. We construct an algorithm in order to extract the bilingual lexicon from the self-similarity matrix and apply a statistical model to estimate the precision, the correctness, of the translations in the bilingual lexicon. The best result is achieved with an application of the vector space model with a precision of about 80 %. This is a good result and can be compared with the precision of about 60 % found in the literature.</p>

In underground construction and tunnelling, the strive for a more time-efficient construction process naturally focuses on the possibilities of reducing the times of waiting between stages of construction. The ability to project shotcrete (sprayed concrete) on a rock surface at an early stage after blasting is vital to the safety during construction and function of e.g. a tunnel. A complication arises when the need for further blasting affects the hardening of newly applied shotcrete. If concrete, cast or sprayed, is exposed to vibrations at an early age while still in the process of hardening, damage that threatens the function of the hard concrete may occur. There is little, or no, established knowledge on the subject and there are no guidelines for practical use. It is concluded from previous investigations that shotcrete can withstand high particle velocity vibrations without being seriously damaged. Shotcrete without reinforcement can survive vibration levels as high as 0.5−1 m/s while sections with loss of bond and ejected rock will occur for vibration velocities higher than 1 m/s. The performance of young and hardened shotcrete exposed to high magnitudes of vibration is here investigated to identify safe distances and shotcrete ages for underground and tunnelling construction, using numerical analyses and comparison with measurements and observations. The work focuses on finding correlations between numerical results, measurement results and observations obtained during tunnelling. The outcome will be guidelines for practical use. The project involves development of sophisticated dynamic finite element models for which the collected information and data will be used as input, accomplished by using the finite ele­ment program Abaqus. The models were evaluated and refined through comparisons between calculated and measured data. First, existing simple engineering models were compared and evaluated through calculations and comparisons with existing data. The first model tested is a structural dynamic model that consists of masses and spring elements. The second is a model built up with finite beam elements interconnected with springs. The third is a one-dimensional elastic stress wave model. The stress response in the shotcrete closest to the rock when exposed to P-waves striking perpendicularly to the shotcrete-rock interface was simulated. Results from a non-destructive laboratory experiment were also used to provide test data for the models. The experiment studied P-wave propagation along a concrete bar, with proper­ties similar to rock. Cement based mortar with properties that resembles shotcrete was applied on one end of the bar with a hammer impacting the other. The shape of the stress waves travelling towards the shotcrete was registered using accelerometers positioned along the bar. Due to the inhomogeneous nature of the rock, the stress waves from the blasting attenuate on the way from the point of explosion towards the shotcrete on the rock surface. Material damping for the rock mass is therefore accounted for, estimated from previous in-situ measurements. The vibration resistance of the shotcrete-rock support system depends on the material properties of the shotcrete and here were age-dependent properties varied to investigate the behaviour of young shotcrete subjected to blast loading. The numerical simulations require insertion of realistic material data for shotcrete and rock, such as density and modulus of elasticity. The calculated results were in good correspondence with observations and measurements in-situ, and with the previous numerical modelling results. Compared to the engineering models, using a sophisticated finite element program facilitate modelling of more complex geometries and also provide more detailed results. It was demonstrated that wave propagation through rock towards shotcrete can be modelled using two dimensional elastic finite elements in a dynamic analysis. The models must include the properties of the rock and the accuracy of the material parameters used will greatly affect the results. It will be possible to describe the propagation of the waves through the rock mass, from the centre of the explosion to the reflection at the shotcrete-rock interface. It is acceptable to use elastic material formulations until the material strengths are exceeded, i.e. until the strains are outside the elastic range, which thus indicates material failure. The higher complexity of this type of model, compared to the engineering models, will make it possible to model more sophisticated geometries. Examples of preliminary recommendations for practical use are given and it is demonstrated how the developed models and suggested analytical technique can be used to obtain further detailed limit values. / Inom undermarks- och tunnelbyggande leder strävan efter en mer tidseffektiv byggprocess till fokus på möjligheten att reducera väntetiderna mellan byggetapper. Möjligheten att projicera sprutbetong på bergytor i ett tidigt skede efter sprängning är avgörande för säkerheten under konstruktionen av t.ex. en tunnel. En komplikation uppstår när behovet av ytterligare sprängning kan påverka härdningen av nysprutad betong. Om betong, gjuten eller sprutad, utsätts för vibrationer i ett tidigt skede under härdningsprocessen kan skador som hotar funktionen hos den hårdnade betongen uppstå. Kunskapen i ämnet är knapphändig, eller obefintlig, och det finns inga etablerade riktlinjer för praktisk användning. Slutsatsen från tidigare undersökningar visar att sprutbetong kan tåla höga vibrationer (partikelhastigheter) utan att allvarliga skador uppstår. Oarmerad sprutbetong kan vara oskadd efter att ha utsatts för så höga vibrationsnivåer som 0,5–1 m/s medan partier med förlorad vidhäftning till berget kan uppträda vid vibrationshastigheter högre än 1 m/s. Funktionen hos ung och hårdnande sprutbetong som utsätts för höga vibrationsnivåer undersöks här för att identifiera säkra avstånd och sprutbetongåldrar för undermarks- och tunnelbyggande, med hjälp av numeriska analyser och jämförelser med mätningar och observationer. Arbetet fokuserar på att finna samband mellan numeriska resultat, mätresultat och observationer från tunnelbyggande. Det slutliga resultatet kommer att vara riktlinjer för praktisk användning. Projektet omfattar utveckling av sofistikerade dynamiska finita elementmodeller för vilka insamlad information och data kommer att användas som indata för det finita elementprogrammet Abaqus. Modellerna utvärderades och förfinats genom jämförelser mellan beräknade och uppmätta resultat. Först jämfördes befintliga enkla mekaniska, ingenjörsmässiga modeller vilka utvärderades genom beräkningar och jämförelser med befintliga data. Den första modellen är en strukturdynamisk modell bestående av massor och fjäderelement. Den andra är en modell uppbyggd av finita balkelementet sammankopplade med fjädrar. Den tredje är en endimensionell elastisk spänningsvågsmodell. Spänningstillståndet i sprutbetongen närmast berget, utsatt för vinkelrätt inkommande P-vågor simulerades. Resultat från icke-förstörande laborationsprovningar användes också som testdata för modellerna. Experimentellt studerades P-vågsutbredning i en betongbalk med egenskaper likvärdiga med berg. Cementbruk med egenskaper liknande sprutbetong applicerades på balkens ena ände medan en hammare användes i den andra. Formen hos den genererade spänningsvågen som propagerade mot sprutbetongänden registrerades med accelerometrar utplacerade längs balken. På grund av bergets inhomogena karaktär kommer spänningsvågorna från sprängningen att dämpas ut på vägen från detonationspunkten till sprutbetongen på bergytan. Materialdämpningen hos bergmassan måste därför beaktas och kan uppskattas utifrån resultat från fältmätningar. Vibrationståligheten hos förstärkningssystem av sprutbetong och berg beror av materialegenskaperna hos sprutbetongen och i den här studien varierades de åldersberoende egenskaperna för att undersöka beteendet hos nysprutad och hårdnande sprutbetong utsatt för sprängbelastning. De numeriska simuleringarna kräver realistiska materialdata för sprutbetong och berg, som t.ex. densitet och elasticitetsmodul. De beräknade resultaten var i god överensstämmelse med observationer och mätningar i fält plats, och med de tidigare numeriska resultaten. Jämfört med de mekaniska modellerna kan ett sofistikerat finit elementprogram underlätta modellering av mer komplexa geometrier och också ge mer detaljerade resultat. Det visades att vågutbredning genom berg och fram emot sprutbetong kan beskrivas med tvådimensionella elastiska finita element i en dynamisk analys. Modellerna måste beskriva bergets egenskaper och noggrannheten i dessa parametrar kommer att ha stor inverkan på resultaten. Det kommer att vara möjligt att beskriva vågutbredningen genom bergmassan, från detonationspunkten till reflektionen vid bergytan, det vill säga gränssnittet mellan sprutbetong och berg. Det är acceptabelt att använda elastiska materialformuleringar tills materialet elasticitetsgräns överskrids tills töjningar utanför det elastiska området nås, vilket därmed indikerar materialbrott. Den högre komplexiteten hos den här typen av modell, jämfört med de mekaniska modellerna, kommer att göra det möjligt att analysera mer komplexa tunnelgeometrier. Exempel på preliminära rekommendationer för praktiskt bruk ges och det visas hur de utvecklade modellerna och den föreslagna analysmetodiken kan användas för att fastställa ytterligare detaljerade gränsvärden. / QC 20120427

This dissertation is focused on theoretical and experimental studies of optical properties of materials and multilayer structures composing liquid crystal displays (LCDs) and electrochromic (EC) devices. By applying spectroscopic ellipsometry, we have determined the optical constants of thin films of electrochromic tungsten oxide (WOx) and nickel oxide (NiOy), the films’ thickness and roughness. These films, which were obtained at spattering conditions possess high transmittance that is important for achieving good visibility and high contrast in an EC device. Another application of the general spectroscopic ellipsometry relates to the study of a photo-alignment layer of a mixture of azo-dyes SD-1 and SDA-2. We have found the optical constants of this mixture before and after illuminating it by polarized UV light. The results obtained confirm the diffusion model to explain the formation of the photo-induced order in azo-dye films. We have developed new techniques for fast characterization of twisted nematic LC cells in transmissive and reflective modes. Our techniques are based on the characteristics functions that we have introduced for determination of parameters of non-uniform birefringent media. These characteristic functions are found by simple procedures and can be utilised for simultaneous determination of retardation, its wavelength dispersion, and twist angle, as well as for solving associated optimization problems. Cholesteric LCD that possesses some unique properties, such as bistability and good selective scattering, however, has a disadvantage – relatively high driving voltage (tens of volts). The way we propose to reduce the driving voltage consists of applying a stack of thin (~1µm) LC layers. We have studied the ability of a layer of a surface stabilized ferroelectric liquid crystal coupled with several retardation plates for birefringent color generation. We have demonstrated that in order to accomplish good color characteristics and high brightness of the display, one or two retardation plates are sufficient.

<p>Detta examensarbete presenterar en utvärdering av konsekvenser som kan uppstå vid Införandet av Multiplex Power-mättmetoden för att kontrollera FM-deviation vid radiosändningar i Sverige.</p><p>Metoden är beskriven i en publicerad standard från ITU: ITU-R BS.412-9. Metoden, även kallad Multiplex Power, är en begränsningsstandard för analog (terrestrial) FM sändning på VHF band. Skälet att införa Multiplex Power mätning är att man genom metoden både kan minska grannkanalinterferenserna och samtidigt jämna ut de stora variationerna i de upplevda ljudnivåerna som finns mellan de olika FM-stationerna. Denna mätmetod går ut på att man samplar värden på deviation, lagrar samplen och räknar ut ett medelvärde av energin under ett intervall av 60 sekunder. Detta energivärde jämförs med en referensnivå, som i ITU-normen definieras som +/-19kHz toppdeviation vid 400 Hz modulationssignal. Skillnaden detta referensvärde kan anges i dB av instrumentet, som en enskild siffra vilkets värde skrivs dynamiskt över tid.</p><p>För att kunna utvärdera metoden genomfördes ett lyssningsexperiment. I experimentet användes ljudfiler som är processade enligt hur dagens normer tillämpas i Sverige och ljudfiler som är processade för att följa BS.412-9-normen. Som tillvägagångssätt användes ett webbformulär och en databas. Testpersoner kunde utföra lyssningen och notera resultatet som sedan lagrades för vidare beräkningar. Ljudfilerna och tillhörande mätningsvärden från lyssningsexperimentet användes sedan för att utvärdera resultatet.</p> / QC 20100707

<p>In this thesis we study hybrid numerical methods to be used in computational electromagnetics. The purpose is to address a wide frequency range relative to a given geometry. We also focus on efficient and robust numerical algorithms for computing the so called Smooth Surface Diffraction predicted by Geometrical Theory of Diffraction (GTD). We restrict the presentation to frequency domain scattering problems.</p><p>The hybrid methods consist in combinations of Boundary Element Methods and asymptotic methods. Three hybrids will be presented. One of them has been developed from a theoretical idea to an industrial code. The two other hybrids will be presented mainly from a theoretical perspective.</p><p>To be able to compute the Smooth Surface Diffracted field we introduce a numerical method that is to be used with surface curvature sensitive meshing, complemented with auxiliary data taken from a geometry database. By using two geometry representations we can show first order convergence and we then achieve an efficient and robust numerical algorithm. This numerical algorithm may be an essential part of an GTD implementation which in its turn is a component in the hybrid methods.</p><p>As a background to our new techiniques we will also give short introductions to the Boundary Element Method and the Geometrical Theory of Diffraction from a theoretical and implementational point of view.</p>

<p>In this report, we show how to implement direct ray tracing of B´ezier surfaces on graphics processing units (GPUs), in particular bicubic rectangular Bézier surfaces and nonparametric cubic Bézier triangles. We use Newton’s method for the rectangular case and show how to use this method to find the ray-surface intersection. For Newton’s method to work we must build a spatial partitioning hierarchy around each surface patch, and in general, hierarchies are essential to speed up the process of ray tracing. We have chosen to use bounding box hierarchies and show how to implement stackless traversal of such a structure on a GPU. For the nonparametric triangular case, we show how to find the wanted intersection by simply solving a cubic polynomial. Because of the limited precision of current GPUs, we also propose a numerical approach to solve the problem, using a one-dimensional Newton search.</p>

<p> The inverse heat conduction problem also frequently referred as the sideways heat equation, in short SHE, is considered as a mathematical model for a real application, where it is desirable for someone to determine the temperature on the surface of a body. Since the surface itself is inaccessible for measurements, one is restricted to use temperature data from the interior measurements. From a  mathematical point of view, the entire situation leads to a non-characteristic Cauchy problem, where by using recorded temperature one can solve a well-posed nonlinear problem in the finite region for computing heat flux, and consequently obtain the Cauchy data [u, u_x]. Further by using these data and by performing an appropriate method, e.g. a space marching method, one can eventually achieve the desired temperature at x = 0.</p><p>The problem is severely ill-posed in the sense that the solution does not depend continuously on the data. The problem solved by two different methods, and for both cases we stabilize the computations by replacing the time derivative in the heat equation by a bounded operator. The first one, a spectral method based on finite Fourier space is illustrated to supply an analytical approach for approximating the time derivative. In order to get a better accuracy in the numerical computation, we use cubic spline function for approximating the time derivative in the least squares sense.</p><p>The inverse problem we want to solve, by using Cauchy data, is a nonlinear heat conduction problem in one space dimension. Since the temperature data u = g(t) is recorded, e.g. by a thermocouple, it usually contains some perturbation in the data. Thus the solution can be severely ill-posed if the Cauchy data become very noisy. Two experiments are presented to test the proposed approach.</p>