The derivation of the standard model from a higher-dimensional action suggests a further study of the fibre bundle formulation of gauge theories to determine the variations in the choice of structure group that are allowed in this geometrical setting. The action of transformations on the projection of fibres to their submanifolds are characteristic of theories with fewer gauge vector bosons, and specific examples are given, which may have phenomenological relevance. The spinor space for the three generations of fermions in the standard model is described algebraically.

In this thesis, we study the expected mean of the torsion angle of an n-stepequilateral random walk in 3D. We consider the random walk is generated within a confining sphere or without a confining sphere: given three consecutive vectors →e1 , →e2 , and →e3 of the random walk then the vectors →e1 and →e2 define a plane and the vectors →e2 and →e3 define a second plane. The angle between the two planes is called the torsion angle of the three vectors. Algorithms are described to generate random walks which are used in a particular space (both without and with confinement). The torsion angle is expressed as a function of six variables for a random walk in both cases: without confinement and with confinement, respectively. Then we find the probability density functions of these six variables of a random walk and demonstrate an explicit integral expression for the expected mean torsion value. Finally, we conclude that the expected torsion angle obtained by the integral agrees with the numerical average torsion obtained by a simulation of random walks with confinement.

In this thesis, we study the expected mean of the torsion angle of an n-stepequilateral random walk in 3D. We consider the random walk is generated within a confining sphere or without a confining sphere: given three consecutive vectors →e1 , →e2 , and →e3 of the random walk then the vectors →e1 and →e2 define a plane and the vectors →e2 and →e3 define a second plane. The angle between the two planes is called the torsion angle of the three vectors. Algorithms are described to generate random walks which are used in a particular space (both without and with confinement). The torsion angle is expressed as a function of six variables for a random walk in both cases: without confinement and with confinement, respectively. Then we find the probability density functions of these six variables of a random walk and demonstrate an explicit integral expression for the expected mean torsion value. Finally, we conclude that the expected torsion angle obtained by the integral agrees with the numerical average torsion obtained by a simulation of random walks with confinement.

I propose two Galerkin methods based on the transfer-of-approximation property for static and dynamic acoustic boundary value problems in seismic applications. For problems with heterogeneous coefficients, the polynomial finite element spaces are no longer optimal unless special meshing techniques are employed. The transfer-of-approximation property provides a general framework to construct the optimal approximation subspace on regular grids. The transfer-of-approximation finite element method is theoretically attractive for that it works for both scalar and vectorial elliptic problems. However the numerical cost is prohibitive. To compute each transfer-of-approximation finite element basis, a problem as hard as the original one has to be solved. Furthermore due to the difficulty of basis localization, the resulting stiffness and mass matrices are dense. The 2D harmonic coordinate finite element method (HCFEM) achieves optimal second-order convergence for static and dynamic acoustic boundary value problems with variable coefficients at the cost of solving two auxiliary elliptic boundary value problems. Unlike the conventional FEM, no special domain partitions, adapted to discontinuity surfaces in coe cients, are required in HCFEM to obtain the optimal convergence rate. The resulting sti ness and mass matrices are constructed in a systematic procedure, and have the same sparsity pattern as those in the standard finite element method. Mass-lumping in HCFEM maintains the optimal order of convergence, due to the smoothness property of acoustic solutions in harmonic coordinates, and overcomes the numerical obstacle of inverting the mass matrix every time update, results in an efficient, explicit time step.

Finite difference and finite element methods exhibit first order convergence when applied to static interface problems where the grid and interface are not aligned. Although modified and unstructured grid methods would address the issue of misalignment for finite elements, application to large models of stratified media, such as those encountered in exploration geophysics, may require not only manual mesh manipulation but also more degrees of freedom than are ultimately necessary to resolve the solution. Instead using fitted or otherwise modified grids, this thesis details an improvement to an existing upscaling method that incorporates fine-scale variations of material properties by composing standard piecewise linear basis functions with a specific type of harmonic map. This technique requires that the problem domain be discretized using two meshes: one fine mesh where the harmonic map is computed to resolve fine-scale structures, and a coarse mesh where the solution to the problem is approximated. The implementation of this method in the literature restricts these composite basis functions to triangular elements in 2D leading to a non-conforming finite element method and suboptimal convergence. However, the support of these basis functions in harmonic coordinates is triangular. I present a mesh-mesh intersection algorithm that exploits this alternative representation to determine the true support of the composite basis functions in terms of the fine mesh. The result is a conforming, high-resolution finite element basis that is associated with the original coarse mesh nodes. Leveraging this fine scale information, I develop a new finite element matrix assembly algorithm. Knowing the shape of the basis support leads naturally to an integration method for computing the finite element matrix entries that is exact up to the accuracy of the harmonic map approximation. This new conforming method is shown to improve the accuracy of solutions to elliptic PDE with discontinuous coefficients on coarse, regular grids.

Complex number is a major mathematical discovery. It can be used in many scientific fields, including engineering, electromagnetism, quantum physics, applied mathematics, and chaos theory. The aim of this paper investigates the problems of algebra, trigonometry and geometry, which can be solved cleverly by the properties of complex numbers.

Dieses Manuskript entstand aus Vorlesungen über Geodätische Berechnungen an der Hochschule für Technik und Wirtschaft Dresden. Da diese Lehrveranstaltung im ersten oder zweiten Semester stattfindet, werden noch keine Methoden der höheren Mathematik benutzt. Das Themenspektrum beschränkt sich deshalb weitgehend auf elementare Berechnungen in der Ebene. Nur im Kapitel 7 kommen einige Methoden der Vektorrechnung zum Einsatz.

The Laplace transform analytic element method (LT-AEM), applies the traditionally steady-state analytic element method (AEM) to the Laplace-transformed diffusion equation (Furman and Neuman, 2003). This strategy preserves the accuracy and elegance of the AEM while extending the method to transient phenomena. The approach taken here utilizes eigenfunction expansion to derive analytic solutions to the modified Helmholtz equation, then back-transforms the LT-AEM results with a numerical inverse Laplace transform algorithm. The two-dimensional elements derived here include the point, circle, line segment, ellipse, and infinite line, corresponding to polar, elliptical and Cartesian coordinates. Each element is derived for the simplest useful case, an impulse response due to a confined, transient, single-aquifer source. The extension of these elements to include effects due to leaky, unconfined, multi-aquifer, wellbore storage, and inertia is shown for a few simple elements (point and line), with ready extension to other elements. General temporal behavior is achieved using convolution between these impulse and general time functions; convolution allows the spatial and temporal components of an element to be handled independently.Comparisons are made between inverse Laplace transform algorithms; the accelerated Fourier series approach of de Hoog et al. (1982) is found to be the most appropriate for LT-AEM applications. An application and synthetic examples are shown for several illustrative forward and parameter estimation simulations to illustrate LT-AEM capabilities. Extension of LT-AEM to three-dimensional flow and non-linear infiltration are discussed.

The pollution generated by accidental marine oil spills can cause persistent ecological disasters and lead to serious social and economical damages. Numerical simulations are a valuable tool to make proper decisions in emergency situation or to plan response actions beforehand. The main objective of this work was to improve SIMOIL, a computational model developed earlier at URV and capable of predicting the evaporation and spreading of massive oil spills in coastal areas. Specifically, a new coastal current model, based on the resolution of the shallow water equations in generalised coordinates, has been developed and validated and then coupled to SIMOIL. The model was specially designed to describe coastal oceanic flows over topography accounting for Coriolis force, eddy viscosity, seabed friction and to couple with SIMOIL in domain with complex boundaries. The equations have been discretized over generalised domains by means of finite differences of second order accuracy. The code was then implemented in FORTRAN. The code has been validated extensively against numerical and experimental flow studies of the bibliography. Finally, the new complete version of SIMOIL, coupling the shallow water model and the oil slick model, has been applied to the study of two accidental oil spills: • A massive leakage from the Repsol's floating dock in the port of Tarragona • The biggest oil spill ever occurred in the Eastern Mediterranean Sea: the 2006 Lebanon oil spill. In both cases, the new version of SIMOIL, demonstrate more accurate predictions of the behaviour of the oil spill, specially for moderate winds with complex topography. / La contaminación generada por los vertidos accidentales de petróleo puede ser reducida si se actúa y si se toman las decisiones adecuadas a tiempo. Las simulaciones numéricas de vertidos de petróleo permiten predecir la evolución de las manchas de crudo. En este trabajo, el objetivo principal era de mejorar la precisión y el rango de aplicación del código SIMOIL desarrollando e integrando al código un modelo de predicción de corrientes marinas en aguas costeras. Se han derivado las ecuaciones de aguas poco profundas en coordenadas generalizadas. Se han discretizado las ecuaciones y el código se implementó en FORTRAN 90. El modelo así como los métodos numéricos han sido validados con el estudio de flujos experimentales y numéricos de la bibliografía. Finalmente, la nueva versión de SIMOIL se aplicó con éxito a dos casos físicos de vertidos de crudo: • un vertido ficticio desde la monoboya de descarga de Repsol en el puerto de Tarragona • un vertido real, el mas grande ocurrido en el Este del mar Mediterráneo, consecuencia de la guerra en Líbano en julio de 2006. En ambos casos la nueva versión de SIMOIL proporcionó predicciones más precisas, especialmente para vientos moderados y topografías complejas.

Analytic solutions to the static and stationary boundary value field problems relative to an arbitrary configuration of parallel cylinders are obtained by using translational addition theorems for scalar Laplacian polar functions, to express the field due to one cylinder in terms of the polar coordinates of the other cylinders such that the boundary conditions can be imposed at all the cylinder surfaces. The constants of integration in the field expressions of all the cylinders are obtained from a truncated infinite matrix equation. Translational addition theorems are available for scalar cylindrical and spherical wave functions but such theorems are not directly available for the general solution of the Laplace equation in polar coordinates. The purpose of deriving these addition theorems and applying them to field problems involving systems of cylinders is to obtain exact analytic solutions with controllable accuracies, thereby, yielding benchmark solutions to validate other approximate numerical methods.