In this thesis, we extend the spectral collocation methods(SCM) (i.e., pseudo-spectral method) in Quarteroni and Valli [27] for the semilinear, parameter-dependentproblems(PDP) in the square with the Dirichlet boundary condition. The optimal error bounds are derived in this thesis for both H1 and L2 norms. For the solutions sufficiently smooth, the very high convergence rates can be obtained. The algorithms of the SCM are simple and easy to carry out. Only a few of basis functions are needed so that not only can the high accuracy of the PDP solutions be achieved, but also a great deal of CPU time may be saved. Moreover, for PDP the stability analysis of SCM is also made, to have the same growth rates of condition number as those for Poisson¡¦s equation. Numerical experiments are carried out to verify the theoretical analysis made.

Yes / In this paper, the Chebychev spectral collocation method is applied for the thermal analysis of convective-radiative straight fins with the temperature-dependent thermal conductivity. The developed heat transfer model was used to analyse the thermal performance, establish the optimum thermal design parameters, and also, investigate the effects of thermo-geometric parameters and thermal conductivity (nonlinear) parameters on the thermal performance of the fin. The results of this study reveal that the rate of heat transfer from the fin increases as the convective, radioactive, and magnetic parameters increase. This study establishes good agreement between the obtained results using Chebychev spectral collocation method and the results obtained using Runge-Kutta method along with shooting, homotopy perturbation, and adomian decomposition methods.

In this paper, we use the newly developed method of particular solution (MPS) and one-stage method of fundamental solution (MFS-MPS) for solving partial differential equation (PDE). In the 1-D Poisson equation, we prove the solution of MFS-MPS is converge to Spectral Collocation Method using Polynomial, and show that the numerical solution similar to those of using the method of particular solution (MPS), Kansa's method, and Spectral Collocation Method using Polynomial (SCMP). In 2-D, we also test these results for the Poisson equation and find the error behaviors.

Many types of radial basis functions, such as multiquadrics, contain a free parameter called shape factor, which controls the flatness of RBFs. In the 1-D problems, Fornberg et al. [2] proved that with simple conditions on the increasingly flat radial basis function, the solutions converge to the Lagrange interpolating. In this report, we study and extend it to the 1-D Poisson equation RBFs direct solver, and observed that the interpolants converge to the Spectral Collocation Method using Polynomial. In 2-D, however, Fornberg et al. [2] observed that limit of interpolants fails to exist in cases of highly regular grid layouts. We also test this in the PDEs solver and found the error behavior is different from interpolating problem.

Dans le cadre de l'étude des écoulements ouverts, deux types de dynamiques coexistent. Les écoulements de type oscillateur qui présentent une fréquence propre d'oscillation indépendante des perturbations extérieures (dynamique intrinsèque), ainsi que les écoulements de type amplificateur sélectif de bruit comme les écoulements de jets ou de couches limites décollées, caractérisés par une plus large gamme de fréquences dépendantes essentiellement de bruit extérieur (dynamique extrinsèque). Les études de couches limites décollées en régime incompressible ont montré un lien entre le phénomène auto-entretenu de basse fréquence qui apparaît et l'interaction non normale des modes globaux instables existants pour ce type de configuration. L'objectif de ce travail consiste à étendre cette interprétation lorsque l'écoulement est en régime subsonique. Dans ce but, un travail d'adaptation des conditions aux limites non-réfléchissantes aux problèmes de stabilité globale a été réalisé. Une méthode de zone absorbante de type Perfectly Matched Layer a été implémentée dans un code de simulation numérique utilisant des méthodes de collocation spectrale. Une méthode de décomposition de domaine adaptée aux calculs des solutions stationnaires ainsi qu'aux problèmes de stabilité globale a également été utilisée pour permettre la validation des conditions aux limites implémentées sur un cas d'écoulement rayonnant de cavité ouverte. Enfin, les études de stabilité d'un écoulement de couche limite décollée derrière une géométrie de type bosse ont été réalisées. L'étude des instabilités bidimensionnelles, responsables du phénomène basse fréquence (flapping), et réalisées en régime subsonique montre que le mécanisme observé en régime incompressible est aussi observé en régime subsonique. La stabilité de cet écoulement vis-à-vis de perturbations tri-dimensionnelles, et plus particulièrement les instabilités centrifuges ont aussi été étudiées en fonction du nombre de Mach. / In open flows context, there are generally two types of dynamic : oscillators, such as cylinder flow, exhibit a well defined frequency insensitive to external perturbations (intrinsic dynamics) and noise amplifiers, such as boundary layers, jets or in some cases the separated flows, which are characterized by wider spectrum bands that depend essentially on the external noise (dynamic extrinsic). Previous studies have shown that separated flows are subject to self-induced oscillations of low frequency in incompressible regime. These studies have revealed links between the interaction of non-normal modes and low oscillations in an incompressible boundary-layer separation and it will be to establish the validity of this interpretation in a compressible regime. In this regard, non-reflecting boundary conditions have been developed to solve the eigenvalue problem formed by linearised Navier-Stokes equations. An absorbing region known as Perfectly Matched Layer has been implemented in order to damp acoustic perturbations which are generated when the compressibility of the flow is considered. A multi-domain approach using spectral collocation discretisation has also been developed in order to study the influence of this absorbing region on the stability analysis of an open cavity flow which is known to generate acoustic perturbations. Finally, we focused on separated boundary layer induced by a bump geometry in order to understand what are the effects of compressibility on the bidimensional low frequency phenomenon and also on transverse instabilities which are known to be unstable for a lots of separated flows.

The computation of Schrödinger equations in the semiclassical regime presents several enduring challenges due to the presence of the small semiclassical parameter. Standard approaches for solving these equations commence with spatial discretisation followed by exponentiation of the discretised Hamiltonian via exponential splittings. In this thesis we follow an alternative strategy${-}$we develop a new technique, called the symmetric Zassenhaus splitting procedure, which involves directly splitting the exponential of the undiscretised Hamiltonian. This technique allows us to design methods that are highly efficient in the semiclassical regime. Our analysis takes place in the Lie algebra generated by multiplicative operators and polynomials of the differential operator. This Lie algebra is completely characterised by Jordan polynomials in the differential operator, which constitute naturally symmetrised differential operators. Combined with the $\mathbb{Z}_2$-graded structure of this Lie algebra, the symmetry results in skew-Hermiticity of the exponents for Zassenhaus-style splittings, resulting in unitary evolution and numerical stability. The properties of commutator simplification and height reduction in these Lie algebras result in a highly effective form of $\textit{asymptotic splitting:} $exponential splittings where consecutive terms are scaled by increasing powers of the small semiclassical parameter. This leads to high accuracy methods whose costs grow quadratically with higher orders of accuracy. Time-dependent potentials are tackled by developing commutator-free Magnus expansions in our Lie algebra, which are subsequently split using the Zassenhaus algorithm. We present two approaches for developing arbitrarily high-order Magnus--Zassenhaus schemes${-}$one where the integrals are discretised using Gauss--Legendre quadrature at the outset and another where integrals are preserved throughout. These schemes feature high accuracy, allow large time steps, and the quadratic growth of their costs is found to be superior to traditional approaches such as Magnus--Lanczos methods and Yoshida splittings based on traditional Magnus expansions that feature nested commutators of matrices. An analysis of these operatorial splittings and expansions is carried out by characterising the highly oscillatory behaviour of the solution.