The complexity of numerical ocean circulation models requires careful checking with a variety of test problems. The purpose of this paper is to develop a test problem involving Rossby and gravity waves in a two-layer fluid in a channel. The goal is to compute very accurate solutions to this test problem. These solutions can then be used as a part of the checking process for numerical ocean circulation models. Here, Chebychev pseudo-spectral methods are used to solve the governing equations with a high degree of accuracy. Chebychev pseudo-spectral methods can be described in the following way: For a given function, find the polynomial interpolant at a particular non-uniform grid. The derivative of this polynomial serves as an approximation to the derivative of the original function. This approximation can then be inserted to differential equations to solve for approximate solutions. Here, the governing equations reduce to an eigenvalue problem with eigenvectors and eigenvalues corresponding to the spatial dependences of modal solutions and the frequencies of those solutions, respectively. The results of this method are checked in two ways. First, the solutions using the Chebychev pseudo-spectral methods are analyzed and are found to exhibit the properties known to belong to physical Rossby and gravity waves. Second, in the special case where the two-layer model degenerates to a one-layer system, some analytic solutions are known. When the numerical solutions are compared to the analytic solutions, they show an exponential rate of convergence. The conclusion is that the solutions computed using the Chebychev pseudo-spectral methods are highly accurate and could be used as a test problem to partially check numerical ocean circulation models. / Graduation date: 2012

Many known extra-solar giant planets lie close to their host stars. Around 60 have their semi-major axes smaller than 0.05 AU. In contrast to planets further out, the vast majority of these close-in planets have low eccentricity orbits. This suggests that their orbits have been circularized likely due to tidal dissipation inside the planets. These exoplanets share with our own Jupiter at least one trait in common: when they are subject to periodic tidal forcing, they behave like a lossy spring, with a tidal ``quality factor'', Q, of order 10^5. This parameter is the ratio between the energy in the tide and the energy dissipated per period. To explain this, a possible solution is resonantly forced internal oscillation. If the frequency of the tidal forcing happens to land on that of an internal eigenmode, this mode can be resonantly excited to a very large amplitude. The damping of such a mode inside the planet may explain the observed Q value. The only normal modes that fall in the frequency range of the tidal forcing (~ few days) are inertial modes, modes restored by the Coriolis force. We present a new numerical technique to solve for inertial modes in a convective, rotating sphere. This technique combines the use of an ellipsoidal coordinate system with a pseudo-spectral method to solve the partial differential equation that governs the inertial oscillations. We show that, this technique produces highly accurate solutions when the density profile is smooth. In particular, the lines of nodes are roughly parallel to the ellipsoidal coordinate axes. In particular, using these accurate solutions, we estimate the resultant tidal dissipation for giant planets, and find that turbulent dissipation of inertial modes in planets with smooth density profiles do not give rise to dissipation as strong as the one observed. We also study inertial modes in density profiles that exhibit discontinuities, as some recent models of Jupiter show. We found that, in this case, our method could not produce convergent solutions for the inertial modes. Additionally, we propose a way to observe inertial modes inside Saturn indirectly, by observing waves in its rings that may be excited by inertial modes inside Saturn.

Many known extra-solar giant planets lie close to their host stars. Around 60 have their semi-major axes smaller than 0.05 AU. In contrast to planets further out, the vast majority of these close-in planets have low eccentricity orbits. This suggests that their orbits have been circularized likely due to tidal dissipation inside the planets. These exoplanets share with our own Jupiter at least one trait in common: when they are subject to periodic tidal forcing, they behave like a lossy spring, with a tidal ``quality factor'', Q, of order 10^5. This parameter is the ratio between the energy in the tide and the energy dissipated per period. To explain this, a possible solution is resonantly forced internal oscillation. If the frequency of the tidal forcing happens to land on that of an internal eigenmode, this mode can be resonantly excited to a very large amplitude. The damping of such a mode inside the planet may explain the observed Q value. The only normal modes that fall in the frequency range of the tidal forcing (~ few days) are inertial modes, modes restored by the Coriolis force. We present a new numerical technique to solve for inertial modes in a convective, rotating sphere. This technique combines the use of an ellipsoidal coordinate system with a pseudo-spectral method to solve the partial differential equation that governs the inertial oscillations. We show that, this technique produces highly accurate solutions when the density profile is smooth. In particular, the lines of nodes are roughly parallel to the ellipsoidal coordinate axes. In particular, using these accurate solutions, we estimate the resultant tidal dissipation for giant planets, and find that turbulent dissipation of inertial modes in planets with smooth density profiles do not give rise to dissipation as strong as the one observed. We also study inertial modes in density profiles that exhibit discontinuities, as some recent models of Jupiter show. We found that, in this case, our method could not produce convergent solutions for the inertial modes. Additionally, we propose a way to observe inertial modes inside Saturn indirectly, by observing waves in its rings that may be excited by inertial modes inside Saturn.

This paper is devoted to studying the problem of constructing Ds-optimal design for d-th degree polynomial regression with analytic weight function on the interval [m-a,m+a],m,a in R. It is demonstrated that the structure of the optimal design depends on d, a and weight function only, as a close to 0. Moreover, the Taylor polynomials of the scaled versions of the optimal support points and weights can be computed via a recursive formula.

The development of wireless communications increases the challenges on antenna performance to improve the capability of the whole system. New fabrication technologies are emerging that not only can improve the performance of components but also provide more options for materials and geometries. One of the advanced technologies, referred to as deep X-ray lithography (XRL), can improve the performance of RF components while providing interesting opportunities for fabrication. Since this fabrication technology enables the objects of high aspect ratio (tall) structure with high accuracy, it offers RF/microwave components some unique advantages, such as higher coupling energy and compacted size. The research presented in that thesis investigates the properties of deep XRL fabricated tall microstrip transmission line and describes some important features such as characteristic impedance, attenuation, and electromagnetic field distribution. Furthermore, since most of traditional feeding structure cannot supply enough coupling energy to excite the low permittivity DRA element (εr≤10), three novel feeding schemes composed by tall microstrip line on exciting dielectric resonator antennas (DRA) with low permittivity are proposed and analyzed in this research. Both simulation and experimental measured results exhibit excellent performance. Additionally, a new simulation approach to realize Dolph-Chebyshev linear series-fed DRA arrays by using the advantages of tall microstrip line feeding structure is proposed. By using a novel T shape feeding scheme, the array exhibits wide band operation due to the low permittivity (εr=5) DRA elements and good radiation pattern due to the novel feeding structure. The tall metal transmission line feed structure and the polymer-based DRA elements could be fabricated in a common process by the deep XRL technology. This thesis firstly illustrates properties and knowledge for both DRA element and the tall transmission line. Then the three novel feeding schemes by using the tall transmission line on exciting the low permittivity DRA are proposed and one of the feeding structures, side coupling feeding, is analyzed through the simulation and experiments. Finally, the T shape feeding structure is applied into low permittivity linear DRA array design work. A novel method on designing the Dolph-Chebyshev array is proposed making the design work more efficient.

The study of nano-electronic devices is fundamental to the advancement of the semiconductor industry. As electronic devices become increasingly smaller, they will eventually move into a regime where the classical nature of the electrons no longer applies. As the quantum nature of the electrons becomes increasingly important, classical or semiclassical theories and methods will no longer serve their purpose. For example, the simplest non-classical effect that will occur is the tunnelling of electrons through the potential barriers that form wires and transistors. This results in an increase in noise and a reduction in the device?s ability to function correctly. Other quantum effects include coulomb blockade, resonant tunnelling, interference and diffraction, coulomb drag, resonant blockade and the list goes on. This thesis develops both a theoretical model and computational method to allow nanoelectronic devices to be studied in detail. Through the use of computer code and an appropriate model description, potential problems and new novel devices may be identified and studied. The model is as accurate to the physical realisation of the devices as possible to allow direct comparison with experimental outcomes. Using simple geometric shapes of varying potential heights, simple devices are readily accessible: quantum wires; quantum transistors; resonant cavities; and coupled quantum wires. Such devices will form the building blocks of future complex devices and thus need to be fully understood. Results obtained studying the connection of a quantum wire with its surroundings demonstrate non-intuitive behaviour and the importance of device geometry to electrical characteristics. The application of magnetic fields to various nano-devices produced a range of interesting phenomenon with promising novel applications. The magnetic field can be used to alter the phase of the electron, modifying the interaction between the electronic potential and the transport electrons. This thesis studies in detail the Aharonov-Bohm oscillation and impurity characterisation in quantum wires. By studying various devices considerable information can be added to the knowledge base of nano-electronic devices and provide a basis to further research. The computational algorithms developed in this thesis are highly accurate, numerically efficient and unconditionally stable, which can also be used to study many other physical phenomena in the quantum world. As an example, the computational algorithms were applied to positron-hydrogen scattering with the results indicating positronium formation.

In this work, the fast evaluation and reconstruction of multivariate trigonometric polynomials with frequencies supported on arbitrary index sets of finite cardinality is considered, where rank-1 lattices are used as spatial discretizations. The approximation of multivariate smooth periodic functions by trigonometric polynomials is studied, based on a one-dimensional FFT applied to function samples. The smoothness of the functions is characterized via the decay of their Fourier coefficients, and various estimates for sampling errors are shown, complemented by numerical tests for up to 25 dimensions. In addition, the special case of perturbed rank-1 lattice nodes is considered, and a fast Taylor expansion based approximation method is developed. One main contribution is the transfer of the methods to the non-periodic case. Multivariate algebraic polynomials in Chebyshev form are used as ansatz functions and rank-1 Chebyshev lattices as spatial discretizations. This strategy allows for using fast algorithms based on a one-dimensional DCT. The smoothness of a function can be characterized via the decay of its Chebyshev coefficients. From this point of view, estimates for sampling errors are shown as well as numerical tests for up to 25 dimensions. A further main contribution is the development of a high-dimensional sparse FFT method based on rank-1 lattice sampling, which allows for determining unknown frequency locations belonging to the approximately largest Fourier or Chebyshev coefficients of a function. / In dieser Arbeit wird die schnelle Auswertung und Rekonstruktion multivariater trigonometrischer Polynome mit Frequenzen aus beliebigen Indexmengen endlicher Kardinalität betrachtet, wobei Rang-1-Gitter (rank-1 lattices) als Diskretisierung im Ortsbereich verwendet werden. Die Approximation multivariater glatter periodischer Funktionen durch trigonometrische Polynome wird untersucht, wobei Approximanten mittels einer eindimensionalen FFT (schnellen Fourier-Transformation) angewandt auf Funktionswerte ermittelt werden. Die Glattheit von Funktionen wird durch den Abfall ihrer Fourier-Koeffizienten charakterisiert und mehrere Abschätzungen für den Abtastfehler werden gezeigt, ergänzt durch numerische Tests für bis zu 25 Raumdimensionen. Zusätzlich wird der Spezialfall gestörter Rang-1-Gitter-Knoten betrachtet, und es wird eine schnelle Approximationsmethode basierend auf Taylorentwicklung vorgestellt. Ein wichtiger Beitrag dieser Arbeit ist die Übertragung der Methoden vom periodischen auf den nicht-periodischen Fall. Multivariate algebraische Polynome in Chebyshev-Form werden als Ansatzfunktionen verwendet und sogenannte Rang-1-Chebyshev-Gitter als Diskretisierungen im Ortsbereich. Diese Strategie ermöglicht die Verwendung schneller Algorithmen basierend auf einer eindimensionalen DCT (diskreten Kosinustransformation). Die Glattheit von Funktionen kann durch den Abfall ihrer Chebyshev-Koeffizienten charakterisiert werden. Unter diesem Gesichtspunkt werden Abschätzungen für Abtastfehler gezeigt sowie numerische Tests für bis zu 25 Raumdimensionen. Ein weiterer wichtiger Beitrag ist die Entwicklung einer Methode zur Berechnung einer hochdimensionalen dünnbesetzten FFT basierend auf Abtastwerten an Rang-1-Gittern, wobei diese Methode die Bestimmung unbekannter Frequenzen ermöglicht, welche zu den näherungsweise größten Fourier- oder Chebyshev-Koeffizienten einer Funktion gehören.

No description available.

No description available.

Queueing theory has often been applied to study communication and service queueing systems such as call centers, hospital emergency departments and ride-sharing platforms. Unfortunately, it is complicated to analyze queueing systems. That is largely because the arrival and service processes that mainly determine a queueing system are uncertain and must be represented as stochastic processes that are difficult to analyze. In response, service providers might be able to partially capture the main characteristics of systems given partial data information and limited domain knowledge. An effective engineering response is to develop tractable approximations to approximate queueing characteristics of interest that depend on critical partial information. In this thesis, we contribute to developing high-quality approximations by studying tight bounds for the transient and the steady-state mean waiting time given partial information. We focus on single-server queues and multi-server queues with the unlimited waiting room, the first-come-first-served service discipline, and independent sequences of independent and identically distributed sequences of interarrival times and service times. We assume some partial information is known, e.g., the first two moments of inter-arrival and service time distributions. For the single-server GI/GI/1 model, we first study the tight upper bounds for the mean and higher moments of the steady-state waiting time given the first two moments of the inter-arrival time and service-time distributions. We apply the theory of Tchebycheff systems to obtain sufficient conditions for classical two-point distributions to yield the extreme values. For the tight upper bound of the transient mean waiting time, we formulate the problem as a non-convex non-linear program, derive the gradient of the transient mean waiting time over distributions with finite support, and apply classical non-linear programming theory to characterize stationary points. We then develop and apply a stochastic variant of the conditional gradient algorithm to find a stationary point for any given service-time distribution. We also establish necessary conditions and sufficient conditions for stationary points to be three-point distributions or special two-point distributions. Our studies indicate that the tight upper bound for the steady-state mean waiting time is attained asymptotically by two-point distributions as the upper mass point of the service-time distribution increases and the probability decreases, while one mass of the inter-arrival time distribution is fixed at 0. We then develop effective numerical and simulation algorithms to compute the tight upper bound. The algorithms are aided by reductions of the special queues with extremal inter-arrival time and extremal service-time distributions to D/GI/1 and GI/D/1 models. Combining these reductions yields an overall representation in terms of a D/RS(D)/1 discrete-time model involving a geometric random sum of deterministic random variables, where the two deterministic random variables have different values, so that the extremal waiting times need not have a lattice distribution. We finally evaluate the tight upper bound to show that it offers a significant improvement over established bounds. In order to understand queueing performance given only partial information, we propose determining intervals of likely performance measures given that limited information. We illustrate this approach for the steady-state waiting time distribution in the GI/GI/K queue given the first two moments of the inter-arrival time and service time distributions plus additional information about these underlying distributions, including support bounds, higher moments, and Laplace transform values. As a theoretical basis, we apply the theory of Tchebycheff systems to determine extremal models (yielding tight upper and lower bounds) on the asymptotic decay rate of the steady-state waiting-time tail probability, as in the Kingman-Lundberg bound and large deviations asymptotics. We then can use these extremal models to indicate likely intervals of other performance measures. We illustrate by constructing such intervals of likely mean waiting times. Without extra information, the extremal models involve two-point distributions, which yield a wide range for the mean. Adding constraints on the third moment and a transform value produces three-point extremal distributions, which significantly reduce the range, yielding practical levels of accuracy.