No description available.

Waveletteori er en forholdsvis ny matematisk disciplin, som har vakt stor interesse indenfor b°ade teoretisk og anvendt matematik i løbet af det seneste °arti. De altafgørende egenskaber ved wavelets er at de kan analysere forskellige dele af en funktion p°a forskellige skalatrin, samt at de kan repræsentere polynomier nøjagtigt op til en given grad. Dette fører til, at funktioner med hurtige oscillationer eller singulariteter indenfor lokaliserede omr°ader kan approksimeres godt med en linearkombination af forholdsvis f°a wavelets. Til sammenligning skal man medtage mange led i en Fourierrække for at opn°a en god tilnærmelse til den slags funktioner. Disse egenskaber ved wavelets har med held været anvendt indenfor signalbehandling. Denne afhandling omhandler wavelets rolle indenfor scientific computing og den best°ar af tre dele: Del I giver en gennemgang af teorien for ortogonale, kompakt støttede wavelets med udgangspunkt i multiskala analyse. S°adanne wavelets er særligt attraktive, fordi de giver anledning til en stabil og særdeles effektiv algoritme, kaldet den hurtige wavelet transformation (FWT). Vi giver estimater for approksimationsegenskaberne af wavelets og demonstrerer, hvordan og hvorfor FWT-algoritmen kan bruges som første led i en effektiv billedkomprimerings metode. Del II omhandler forskellige implementeringer af FWT algoritmen p°a vektorcomputere og parallelle datamater. Vi udvikler en effektiv og skalerbar parallel FWT algoritme og angiver en model for dens ydeevne. Del III omfatter et studium af mulighederne for at bruge wavelets særlige egenskaber til at løse partielle differentialligninger numerisk. Flere forskellige tilgange identificeres og to af dem beskrives detaljeret. De udviklede algoritmer anvendes p°a den ikke-lineære Schr¨odinger ligning og Burgers ligning. Numeriske undersøgelser viser, at algoritmerne kan være effektive under forudsætning af at problemerne er store, at løsningerne er stærkt lokaliserede og at de forskellige numeriske metode-parametre kan vælges p°a passende vis afhængigt af det p°agældende problem.

When evaluating integration schemes used in molecular dynamics (MD) simulations, energy conservation is often cited as the primary criterion by which the integrators should be com- pared. As a result variable stepsize Runge-Kutta methods are often ruled out of consideration due to their characteristic energy drift. We have shown that by appropriately modifying the stepsize selection strategy in a variable stepsize RK method it is possible for the MD practitioner to obtain substantial control over the energy drift during the course of a simulation. This ability has been previously unreported in the literature, and we present numerical examples to illustrate that it can be achieved without sacrificing computational efficiency under currently obtainable timescales.

When evaluating integration schemes used in molecular dynamics (MD) simulations, energy conservation is often cited as the primary criterion by which the integrators should be com- pared. As a result variable stepsize Runge-Kutta methods are often ruled out of consideration due to their characteristic energy drift. We have shown that by appropriately modifying the stepsize selection strategy in a variable stepsize RK method it is possible for the MD practitioner to obtain substantial control over the energy drift during the course of a simulation. This ability has been previously unreported in the literature, and we present numerical examples to illustrate that it can be achieved without sacrificing computational efficiency under currently obtainable timescales.

In this thesis study, a fast implicit iteration scheme called Alternating Cell Directions Imp licit method is combined with Approximate Factorization scheme. This application aims to offer a mathematically well defined version of the Alternating Cell Directions Implicit Method and increase the accuracy of the iteration scheme that is being used for the numerical solutions of the partial differential equations. The iteration scheme presented here is tested using unsteady diffusion equation, Laplace equation and advection-diffusion equation. The accuracy, convergence character and the stability character of the scheme compared with suitable iteration schemes for structured and unstructured quadrilateral grids. Besides, it is shown that the proposed scheme is applicable to triangular and hybrid polygonal grids. A transonic full potential solver is generated using the current scheme. The flow around a 2-D cylinder is solved for subcritical and supercritical cases. Axi-symmetric flow around cylinder is selected as a benchmark problem since the potential flow around bodies with a blunt leading edge is a more challenging problem than slender bodies. Besides, it is shown that, the method is naturally appropriate for parallelization using shared memory approach without using domain decomposition applications. The parallelization that is performed here is partially line, partially point parallelization. T he performance of the application is presented for a 3-D unsteady diffusion problem using Cartesian cells and 2-D unsteady diffusion problem using both structured and unstructured quadrilateral cells.

Have you ever wondered if you are better than the average bowler? If so, there are a variety of ways to compute the average score of a bowling game, including methods that account for a bowler’s skill level. In this thesis, we discuss several different ways to generate bowling scores randomly. For each distribution, we give results for the expected value and standard deviation of each frame's score, the expected value of the game’s final score, and the correlation coefficient between the score of the first and second roll of a single frame. Furthermore, we shall generalize the results in each distribution for an frame game on pins. Additionally, we shall generalize the number of possible games when bowling frames on pins. Then, we shall derive the frequency distribution of each frame’s scores and the arithmetic mean for frames on pins. Finally, to summarize the variety of distributions, we shall make tables that display the results obtained from each distribution used to model a particular bowler’s score. We evaluate the special case when bowling 10 frames on 10 pins, which represents a standard bowling game.

No description available.

No description available.

Nuestro trabajo se ha concentrado en la simulación y el análisis de modelos no-lineales aplicados a las curvas de luz ópticas de estrellas variables de largo periodo. Cuando las estrellas de masa pequeña o intermedia (de 1 a 11 masas solares) llegan a las fases finales de su evolución, tales como la Rama Asintótica de las (Super-)Gigantes, presentan oscilaciones regulares e irregulares de larga duración (i.e. estrellas Miras). Las pulsaciones regulares e irregulares facilitan los episodios de pérdida de masa y hacen que estas estrellas jueguen un papel crucial en la evolución química de las galaxias. Los modelos hidrodinámicos detallados existentes en la literatura proporcionan resultados muy completos, pero a menudo son difíciles de interpretar. Por lo tanto se ha creado la necesidad de entender, a un nivel más básico, los procesos que llevan a pulsaciones irregulares o comportamientos en contradicción con las observaciones. En este sentido, hemos considerado la estrella variable como constituida por un nucleo compacto y una envoltura cuyo movimiento forzado está dado por las ondas de presión originadas en el interior de la estrella. En el presente enfoque, aproximamos la envoltura mediante una sola capa y por lo tanto, los modelos basados en esta aproximación llevan el nombre de modelos de una capa. Entre ellos,los que hemos estudiado a lo largo de los últimos tres años dan cuenta de dichas irregularidades a través de un sistema de ecuaciones diferenciales que incluye, en las versiones más sencillas, la ecuación del movimiento y la conservación de masa, a las cuales hemos añadido a lo largo del estudio, la ecuación de la variación y transporte de energia en la envoltura estelar. En la presente tesis presentamos los resultados obtenidos, asi como una discusión sobre el estado de la modelización de estrellas variables en el contexto de los modelos de una capa. En el primer capítulo hemos introducido los tipos de estrellas variables y los mecanismos que llevan a su variabilidad, mientras que en el segundo capítulo hemos hecho un breve repaso de los modelos de una capa existentes en la literatura. Los modelos mencionados han sido seleccionados por su relación con los modelos que proponemos. En el primer estudio realizado que describimos en el tercer capítulo, hemos utilizado como punto de partida los resultados obtenidos por Icke et al., 1992, A&A, 258, 341, sobre un oscilador forzado que contiene la mínima dinámica necesaria para describir las oscilaciones estelares: la ecuación del movimiento y la conservación de masa. Aunque ideado para la variabilidad de estrellas de masa pequeña, nosotros hemos extendido el modelo para describir también las estrellas variables más masivas, en la fase de la Rama Asintótica de las Super-Gigantes. Hemos llevado a cabo un estudio paramétrico concienzudo a fin de investigar los tipos de com portamiento proporcionados por el modelo, identificando las bifurcaciones que producen dichos comportamientos y los rangos de los parámetros asociados a ellas. Para ello, hemos integrado en nuestro análisis los métodos característicos de la teoria de las bifurcaciones y del análisis tiempo-frecuencia, con el objetivo de determinar los posibles escenarios de transición al caos. Desde un punto de vista matemático, el análisis ha supuesto el estudio del mapa de Poincaré asociado a nuestro sistema, que, como se ha mencionado, está caracterizado por una perturbación periódica. Los resultados incluyen una serie de bifurcaciones locales y globales, entre las cuales la más importante es una triplicación. Entre las consecuencias de esta bifurcación mencionamos la adquisición por parte del mapa de Poincaré de la propiedad de nontwist que conlleva unas características del mapa de Poincaré típicas de los mapas nontwist (e.g. reconexión, meandros). Debido a la particular forma de la perturbación, la dinámica del sistema se diferencia de la del mapa cubico de Hénon, que se considera el prototipo de los mapas nontwist, y hacen del sistema investigado un ejemplo de dinámica prevista teóricamente, pero para la cual no se conocía ningún ejemplo. Desde un punto de vista astrofísico, la comparación con los resultados obtenidos por Icke et al. (1992) nos ha llevado a concluir que las estrellas variables más masivas presentan pulsaciones más irregulares que las estrella de masa pequeña, en acuerdo con las observaciones. El comportamiento irregular lleva a una pronunciada pérdida de masa, resultado comprobado por los datos observacionales. En el cuarto capítulo hemos extendido el modelo descrito anteriormente, añadiendole la ecuación del transporte de energia que permite una mejor comparación con las observaciones. Los resultados proporcionados por el modelo ampliado, entre los cuales destacamos las series temporales asociadas a la fluctuación de la luminosidad estelar, presentan sorprendentes similitudes con algunas de las más estudiadas y peculiares estrellas variables de largo periodo (las estrellas Miras). Adicionalmente, la dinámica encontrada conduce a series temporales en forma de pulsos energéticos a intervalos de tiempo del orden de mil años, que se pueden relacionar con las periodicidades encontradas en las capas circumestelares que rodean a ciertas nebulosas planetarias. En el quinto capítulo de la presente tesis hemos introducido el acoplamiento entre la convección y la pulsación estelar, proceso que se considera en la literatura como indispensable para una correcta modelización de la evolución de estas estrellas. Para ello hemos utilizado el modelo convectivo de una capa introducido en Stellingwerf, R.F., 1986, ApJ, 303, 119. Los resultados de dicho artículo nos han llamado la atención por algunas discrepancias relacionadas con el análisis de la morfologia de las curvas de luz y velocidad para casos en los cuales no existían ciclos límite. Por consiguiente, nuestra investigación se ha concentrado en el estudio parametrico del sistema con la intención de identificar las condiciones necesarias para la existencia de ciclos límite. Una vez identificadas la regiones de ciclos límite en el espacio parametrico, la morfologia de las curvas de luz asociadas ha revelado la existencia de una banda de inestabilidad semejante a la banda de inestabilidad de las variables Cefeidas. Adicionalmente, hemos ampliado el modelo considerando una forma más realista para los factores geometricos que describen el estado evolutivo de la estrella. Los resultados de esta ampliación indican una progresión de las curvas de luz - o para ser más exacto, de sus amplitudes y periodos - que recuerda la llamada Progresión Hertzsprung de las Cefeidas de tipo Bump. Hemos identificado también los tipos de contrapartidas observacionales que son susceptibles de ajustarse al modelo estudiado. Aunque existen muchos trabajos en la literatura basados en el modelo de Stellingwerf et al. (1986), consideramos que todavía muchas facetas e ideas quedan por aclarar y desarrollar.Para concluir, hemos comentado las implicaciones de nuestros resultados y estudios en el contexto de los modelos no lineales de una capa evidenciando sus capacidades y sus límites.

No description available.