An economical algorithm for the solution of elliptic difference equations independent of user-supplied parameters,

Extra t.p., with thesis statement, inserted.

A posteriori error estimation for discontinuous Galerkin discretizations of H(curl)-elliptic partial differential equations

We develop the a posteriori error estimation of interior penalty discontinuous Galerkin discretizations for H(curl)-elliptic problems that arise in eddy current models. Computable upper and lower bounds on the error measured in terms of a natural (mesh-dependent) energy norm are derived. The proposed a posteriori error estimator is validated by numerical experiments, illustrating its reliability and efficiency for a range of test problems.

Maximum Principle for Elliptic and Parabolic Equations

Nel primo capitolo si riporta il principio del massimo per operatori ellittici. Sarà considerato, in un primo momento, l'operatore di Laplace e, successivamente, gli operatori ellittici del secondo ordine, per i quali si dimostrerà anche il principio del massimo di Hopf. Nel secondo capitolo si affronta il principio del massimo per operatori parabolici e lo si utilizza per dimostrare l'unicità delle soluzioni di problemi ai valori al contorno.

Numerical solution of an elliptic 3-dimensional Cauchy problem by the alternating method and boundary integral equations

We consider the Cauchy problem for the Laplace equation in 3-dimensional doubly-connected domains, that is the reconstruction of a harmonic function from knowledge of the function values and normal derivative on the outer of two closed boundary surfaces. We employ the alternating iterative method, which is a regularizing procedure for the stable determination of the solution. In each iteration step, mixed boundary value problems are solved. The solution to each mixed problem is represented as a sum of two single-layer potentials giving two unknown densities (one for each of the two boundary surfaces) to determine; matching the given boundary data gives a system of boundary integral equations to be solved for the densities. For the discretisation, Weinert's method [24] is employed, which generates a Galerkin-type procedure for the numerical solution via rewriting the boundary integrals over the unit sphere and expanding the densities in terms of spherical harmonics. Numerical results are included as well.

Elliptic curves, group law, and efficient computation

This thesis is about the derivation of the addition law on an arbitrary elliptic curve and efficiently adding points on this elliptic curve using the derived addition law. The outcomes of this research guarantee practical speedups in higher level operations which depend on point additions. In particular, the contributions immediately find applications in cryptology. Mastered by the 19th century mathematicians, the study of the theory of elliptic curves has been active for decades. Elliptic curves over finite fields made their way into public key cryptography in late 1980’s with independent proposals by Miller [Mil86] and Koblitz [Kob87]. Elliptic Curve Cryptography (ECC), following Miller’s and Koblitz’s proposals, employs the group of rational points on an elliptic curve in building discrete logarithm based public key cryptosystems. Starting from late 1990’s, the emergence of the ECC market has boosted the research in computational aspects of elliptic curves. This thesis falls into this same area of research where the main aim is to speed up the additions of rational points on an arbitrary elliptic curve (over a field of large characteristic). The outcomes of this work can be used to speed up applications which are based on elliptic curves, including cryptographic applications in ECC. The aforementioned goals of this thesis are achieved in five main steps. As the first step, this thesis brings together several algebraic tools in order to derive the unique group law of an elliptic curve. This step also includes an investigation of recent computer algebra packages relating to their capabilities. Although the group law is unique, its evaluation can be performed using abundant (in fact infinitely many) formulae. As the second step, this thesis progresses the finding of the best formulae for efficient addition of points. In the third step, the group law is stated explicitly by handling all possible summands. The fourth step presents the algorithms to be used for efficient point additions. In the fifth and final step, optimized software implementations of the proposed algorithms are presented in order to show that theoretical speedups of step four can be practically obtained. In each of the five steps, this thesis focuses on five forms of elliptic curves over finite fields of large characteristic. A list of these forms and their defining equations are given as follows: (a) Short Weierstrass form, y2 = x3 + ax + b, (b) Extended Jacobi quartic form, y2 = dx4 + 2ax2 + 1, (c) Twisted Hessian form, ax3 + y3 + 1 = dxy, (d) Twisted Edwards form, ax2 + y2 = 1 + dx2y2, (e) Twisted Jacobi intersection form, bs2 + c2 = 1, as2 + d2 = 1, These forms are the most promising candidates for efficient computations and thus considered in this work. Nevertheless, the methods employed in this thesis are capable of handling arbitrary elliptic curves. From a high level point of view, the following outcomes are achieved in this thesis. - Related literature results are brought together and further revisited. For most of the cases several missed formulae, algorithms, and efficient point representations are discovered. - Analogies are made among all studied forms. For instance, it is shown that two sets of affine addition formulae are sufficient to cover all possible affine inputs as long as the output is also an affine point in any of these forms. In the literature, many special cases, especially interactions with points at infinity were omitted from discussion. This thesis handles all of the possibilities. - Several new point doubling/addition formulae and algorithms are introduced, which are more efficient than the existing alternatives in the literature. Most notably, the speed of extended Jacobi quartic, twisted Edwards, and Jacobi intersection forms are improved. New unified addition formulae are proposed for short Weierstrass form. New coordinate systems are studied for the first time. - An optimized implementation is developed using a combination of generic x86-64 assembly instructions and the plain C language. The practical advantages of the proposed algorithms are supported by computer experiments. - All formulae, presented in the body of this thesis, are checked for correctness using computer algebra scripts together with details on register allocations.

Adaptive finite element analysis of elliptic problems based on bubble-type local mesh generation

A new mesh adaptivity algorithm that combines a posteriori error estimation with bubble-type local mesh generation (BLMG) strategy for elliptic differential equations is proposed. The size function used in the BLMG is defined on each vertex during the adaptive process based on the obtained error estimator. In order to avoid the excessive coarsening and refining in each iterative step, two factor thresholds are introduced in the size function. The advantages of the BLMG-based adaptive finite element method, compared with other known methods, are given as follows: the refining and coarsening are obtained fluently in the same framework; the local a posteriori error estimation is easy to implement through the adjacency list of the BLMG method; at all levels of refinement, the updated triangles remain very well shaped, even if the mesh size at any particular refinement level varies by several orders of magnitude. Several numerical examples with singularities for the elliptic problems, where the explicit error estimators are used, verify the efficiency of the algorithm. The analysis for the parameters introduced in the size function shows that the algorithm has good flexibility.

Use of Abel integral equations in water wave scattering by two surface-piercing barriers

In this paper the classical problem of water wave scattering by two partially immersed plane vertical barriers submerged in deep water up to the same depth is investigated. This problem has an exact but complicated solution and an approximate solution in the literature of linearised theory of water waves. Using the Havelock expansion for the water wave potential, the problem is reduced here to solving Abel integral equations having exact solutions. Utilising these solutions,two sets of expressions for the reflection and transmission coefficients are obtained in closed forms in terms of computable integrals in contrast to the results given in the literature which,involved six complicated integrals in terms of elliptic functions. The two different expressions for each coefficient produce almost the same numerical results although it has not been possible to prove their equivalence analytically. The reflection coefficient is depicted against the wave number in a number of figures which almost coincide with the figures available in the literature wherein the problem was solved approximately by employing complementary approximations. (C) 2009 Elsevier B.V. All rights reserved.

Free convection about cylinders of elliptic cross section embedded in a porous medium

Numerical solutions are presented for the free convection boundary layers over cylinders of elliptic cross section embedded in a fluid-saturated porous medium. The transformed conservation equations of the nonsimilar boundary layers are solved numerically by an efficient finite-difference method. The theory was applied to a number of cylinders and the results compared very well with published analytical solutions. The results are of use in the design of underground electrical cables, power plant steam, and water distribution lines, among others.

Solutions of cubic equations in quadratic fields

Let K be any quadratic field with O-K its ring of integers. We study the solutions of cubic equations, which represent elliptic curves defined over Q, in quadratic fields and prove some interesting results regarding the solutions by using elementary tools. As an application we consider the Diophantine equation r + s + t = rst = 1 in O-K. This Diophantine equation gives an elliptic curve defined over Q with finite Mordell-Weil group. Using our study of the solutions of cubic equations in quadratic fields we present a simple proof of the fact that except for the ring of integers of Q(i) and Q(root 2), this Diophantine equation is not solvable in the ring of integers of any other quadratic fields, which is already proved in [4].

Strongly magnetized cold degenerate electron gas: Mass-radius relation of the magnetized white dwarf

We consider a relativistic, degenerate electron gas at zero temperature under the influence of a strong, uniform, static magnetic field, neglecting any form of interactions. Since the density of states for the electrons changes due to the presence of the magnetic field (which gives rise to Landau quantization), the corresponding equation of state also gets modified. In order to investigate the effect of very strong magnetic field, we focus only on systems in which a maximum of either one, two, or three Landau level(s) is/are occupied. This is important since, if a very large number of Landau levels are filled, it implies a very low magnetic field strength which yields back Chandrasekhar's celebrated nonmagnetic results. The maximum number of occupied Landau levels is fixed by the correct choice of two parameters, namely, the magnetic field strength and the maximum Fermi energy of the system. We study the equations of state of these one-level, two-level, and three-level systems and compare them by taking three different maximum Fermi energies. We also find the effect of the strong magnetic field on the mass-radius relation of the underlying star composed of the gas stated above. We obtain an exciting result that it is possible to have an electron-degenerate static star, namely, magnetized white dwarfs, with a mass significantly greater than the Chandrasekhar limit in the range 2.3-2.6M(circle dot), provided it has an appropriate magnetic field strength and central density. In fact, recent observations of peculiar type Ia supernovae-SN 2006gz, SN 2007if, SN 2009dc, SN 2003fg-seem to suggest super-Chandrasekhar-mass white dwarfs with masses up to 2.4-2.8M(circle dot) as their most likely progenitors. Interestingly, our results seem to lie within these observational limits.

Elastic fields of interacting elliptic inhomogeneities

This paper presents a method for the calculation of two-dimensional elastic fields in a solid containing any number of inhomogeneities under arbitrary far field loadings. The method called 'pseudo-dislocations method', is illustrated for the solution of interacting elliptic inhomogeneities. It reduces the interacting inhomogeneities problem to a set of linear algebraic equations. Numerical results are presented for a variety of elliptic inhomogeneity arrangements, including the special cases of elliptic holes, cracks and circular inhomogeneities. All these complicated problems can be solved with high accuracy and efficiency.

Magnetohydrodynamics equilibrium of a self-confined elliptic plasma ball

A variational principle is applied to the problem of magnetohydrodynamics (MHD) equilibrium of a self-contained elliptical plasma ball, such as elliptical ball lightning. The principle is appropriate for an approximate solution of partial differential equations with arbitrary boundary shape. The method reduces the partial differential equation to a series of ordinary differential equations and is especially valuable for treating boundaries with nonlinear deformations. The calculations conclude that the pressure distribution and the poloidal current are more uniform in an oblate self-confined plasma ball than that of an elongated plasma ball. The ellipticity of the plasma ball is obviously restricted by its internal pressure, magnetic field, and ambient pressure. Qualitative evidence is presented for the absence of sighting of elongated ball lightning.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Detection of efficient carrier capture in ultrathin InAs/GaAs layers using a degenerate pump-probe technique

By analysing the carrier dynamics based on the rate equations and the change of the refractive index due to the efficient carrier capture, we have calculated the carrier capture process in the InAs/GaAs system detected by a simple degenerate pump-probe technique. The calculated results are found to be in good agreement with the experimental findings. Our results indicate that this simple technique, with the clear advantage of being easy to carry out, can be very useful in studying the carrier dynamics for some specific structures such as InAs ultrathin layers embedded in a GaAs matrix described here.