The pseudospectral method is a family of numerical methods for the solution of differential equations based on the expansion of basis functions defined on a set of grid points. In this thesis, the relationship between the distribution of grid points and the accuracy and convergence of the solution is emphasized. The polynomial and sinc pseudospectral methods are extensively studied along with many applications to quantum and statistical mechanics involving the Fokker-Planck and Schroedinger equations. The grid points used in the polynomial methods coincide with the points of quadrature, which are defined by a set of polynomials orthogonal with respect to a weight function. The choice of the weight function plays an important role in the convergence of the solution. It is observed that rapid convergence is usually achieved when the weight function is chosen to be the square of the ground-state eigenfunction of the problem. The sinc method usually provides a slow convergence as the grid points are uniformly distributed regardless of the behaviour of the solution. For both polynomial and sinc methods, the convergence rate can be improved by redistributing the grid points to more appropriate positions through a transformation of coordinates. The transformation method discussed in this thesis preserves the orthogonality of the basis functions and provides simple expressions for the construction of discretized matrix operators. The convergence rate can be improved by several times in the evaluation of loosely bound eigenstates with an exponential or hyperbolic sine transformation. The transformation can be defined explicitly or implicitly. An explicit transformation is based on a predefined mapping function, while an implicit transformation is constructed by an appropriate set of grid points determined by the behaviour of the solution. The methodologies of these transformations are discussed with some applications to 1D and 2D problems. The implicit transformation is also used as a moving mesh method for the time-dependent Smoluchowski equation when a function with localized behaviour is used as the initial condition. / Science, Faculty of / Mathematics, Department of / Graduate

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

In this thesis we present a numerical algorithm to solve the singular Inverse Sturm-Liouville problems with symmetric potential functions. The singularity, which comes from the unbounded domain of the problem, is treated by considering the limiting case of the associated problem on the symmetric finite interval. In contrast to regular problems which are considered on a finite interval the singular inverse problem has an ill-conditioned structure despite of the limiting treatment. We use the regularization techniques to overcome the ill-posedness difficulty. Moreover, since the problem is nonlinear the iterative solution procedures are needed. Direct computation of the eigenvalues in iterative solution is handled via psoudespectral methods. The numerical examples of the considered problem are given to illustrate the accuracy and convergence behaviour.

In this dissertation the solutions of the dynamics and real-time optimal control of magnetic attitude control and formation flying systems are presented. In magnetic attitude control, magnetic actuators for the time-optimal rest-to-rest maneuver with a pseudospectral algorithm are examined. The time-optimal magnetic control is bang-bang and the optimal slew time is about 232.7 seconds. The start time occurs when the maneuver is symmetric about the maximum field strength. For real-time computations, all the tested samples converge to optimal solutions or feasible solutions. We find the average computation time is about 0.45 seconds with the warm start and 19 seconds with the cold start, which is a great potential for real-time computations. Three-axis magnetic attitude stabilization is achieved by using a pseudospectral control law via the receding horizon control for satellites in eccentric low Earth orbits. The solutions from the pseudospectral control law are in excellent agreement with those obtained from the Riccati equation, but the computation speed improves by one order of magnitude. Numerical solutions show state responses quickly tend to the region where the attitude motion is in the steady state. Approximate models are often used for the study of relative motion of formation flying satellites. A modeling error index is introduced for evaluating and comparing the accuracy of various theories of the relative motion of satellites in order to determine the effect of modeling errors on the various theories. The numerical results show the sequence of the index from high to low should be Hill's equation, non- J2, small eccentricity, Gim-Alfriend state transition matrix index, with the unit sphere approach and the Yan-Alfriend nonlinear method having the lowest index and equivalent performance. A higher order state transition matrix is developed using unit sphere approach in the mean elements space. Based on the state transition matrix analytical control laws for formation flying maintenance and reconfiguration are proposed using low-thrust and impulsive scheme. The control laws are easily derived with high accuracy. Numerical solutions show the control law works well in real-time computations.

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

Many problems in fluid mechanics and material sciences deal with liquid-vapour flows. In these flows, the ideal gas assumption is not accurate and the van der Waals equation of state is usually used. This equation of state is non-convex and causes the solution domain to have two hyperbolic regions separated by an elliptic region. Therefore, the governing equations of these flows have a mixed elliptic-hyperbolic nature. Numerical oscillations usually appear with standard finite-difference space discretization schemes, and they persist when the order of accuracy of the semi-discrete scheme is increased. In this study, we propose to use a Chebyshev pseudospectral method for solving the governing equations. A comparison of the results of this method with very high-order (up to tenth-order accurate) finite difference schemes is presented, which shows that the proposed method leads to a lower level of numerical oscillations than other high-order finite difference schemes, and also does not exhibit fast-traveling packages of short waves which are usually observed in high-order finite difference methods. The proposed method can thus successfully capture various complex regimes of waves and phase transitions in both elliptic and hyperbolic regimes

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.

The physical processes in lakes remain only partially understood despite successful data collection from a variety of sources spanning several decades. Although numerical models are already frequently employed to simulate the physics of lakes, especially in the context of water quality management, improved methods are necessary to better capture the wide array of dynamically important physical processes, spanning length scales from ~ 10 km (basin-scale oscillations) - 1 m (short internal waves). In this thesis, high-order numerical methods are explored for specialized model equations of lakes, so that their use can be taken into consideration in the next generation of more sophisticated models that will better capture important small scale features than their present day counterparts. The full three-dimensional incompressible density-stratified Navier-Stokes equations remain too computationally expensive to be solved for situations that involve both complicated geometries and require resolution of features at length-scales spanning four orders of magnitude. The main source of computational expense lay with the requirement of having to solve a three-dimensional Poisson equation for pressure at every time-step. Simplified model equations are thus the only way that numerical lake modelling can be carried out at present time, and progress can be made by seeking intelligent parameterizations as a means of capturing more physics within the framework of such simplified equation sets. In this thesis, we employ the long-accepted practice of sub-dividing the lake into vertical layers of different constant densities as an approximation to continuous vertical stratification. We build on this approach by including weakly non-hydrostatic dispersive correction terms in the model equations in order to parameterize the effects of small vertical accelerations that are often disregarded by operational models. Favouring the inclusion of weakly non-hydrostatic effects over the more popular hydrostatic approximation allows these models to capture the emergence of small-scale internal wave phenomena, such as internal solitary waves and undular bores, that are missed by purely hydrostatic models. The Fourier and Chebyshev pseudospectral methods are employed for these weakly non-hydrostatic layered models in simple idealized lake geometries, e.g., doubly periodic domains, periodic channels, and annular domains, for a set of test problems relevant to lake dynamics since they offer excellent resolution characteristics at minimal memory costs. This feature makes them an excellent benchmark to compare other methods against. The Discontinuous Galerkin Finite Element Method (DG-FEM) is then explored as a mid- to high-order method that can be used in arbitrary lake geometries. The DG-FEM can be interpreted as a domain-decomposition extension of a polynomial pseudospectral method and shares many of the same attractive features, such as fast convergence rates and the ability to resolve small-scale features with a relatively low number of grid points when compared to a low-order method. The DG-FEM is further complemented by certain desirable attributes it shares with the finite volume method, such as the freedom to specify upwind-biased numerical flux functions for advection-dominated flows, the flexibility to deal with complicated geometries, and the notion that each element (or cell) can be regarded as a control volume for conserved fluid quantities. Practical implementation details of the numerical methods used in this thesis are discussed, and the various modelling and methodology choices that have been made in the course of this work are justified as the difficulties that these choices address are revealed to the reader. Theoretical calculations are intermittently carried out throughout the thesis to help improve intuition in situations where numerical methods alone fall short of giving complete explanations of the physical processes under consideration. The utility of the DG-FEM method beyond purely hyperbolic systems is also a recurring theme in this thesis. The DG-FEM method is applied to dispersive shallow water type systems as well as incompressible flow situations. Furthermore, it is employed for eigenvalue problems where orthogonal bases must be constructed from the eigenspaces of elliptic operators. The technique is applied to the problem calculating the free modes of oscillation in rotating basins with irregular geometries where the corresponding linear operator is not self-adjoint.