The aim of this thesis is to implement source terms to numerical model for curent sheet in solar atmosphere. The chapters are structured in the way that can present us with basic knowledge of the Sun and processes in its atmosphere. Than we present important equations for numerical solution and for the initial equilibrium of our simulation. One of the last chapters is dedicated to software we use for our numerical simulations called FLASH. In the end of the thesis we present results of our numerical simulations.

Made available in DSpace on 2014-10-09T12:38:52Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:02:08Z (GMT). No. of bitstreams: 1 02806.pdf: 2555475 bytes, checksum: dd5ac6dfe2c65a240c1aa4feafd10dbc (MD5) / Dissertacao (Mestrado) / IPEN/D / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

Made available in DSpace on 2014-10-09T12:38:52Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:02:08Z (GMT). No. of bitstreams: 1 02806.pdf: 2555475 bytes, checksum: dd5ac6dfe2c65a240c1aa4feafd10dbc (MD5) / Dissertacao (Mestrado) / IPEN/D / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

Accurate evaluation of the non-linear wave-wave energy transfer requires a significant proportion of the computation time of ocean wave models. The Discrete Interaction Approximation (DIA) developed within the first version of the WAM model (WAMDI, 1988) is the only algorithm to be used today in operational wave modelling as it is the only way to calculate the wave-wave interactions rapidly enough. In this study, the Two- Scale Approximation (TSA), a potential successor method to the DIA, was successfully implemented in the third generation operational wave model WAVEWATCH IIITM(WW3). Preliminary results (Perrie and Resio, 2009) showing that it offers improved accuracy are confirmed in this study by the modelling of wave evolution under constant winds. Fetch-growth curves and two-dimensional spectra for energy and non-linear wave-wave interactions obtained using the TSA in these conditions show better agreement to more exact computations of non-linear interactions, than the DIA results.

No description available.

Numerical schemes for the partial differential equations used to characterize stiffly forced conservation laws are constructed and analyzed. Partial differential equations of this form are found in many physical applications including modeling gas dynamics, fluid flow, and combustion. Many difficulties arise when trying to approximate solutions to stiffly forced conservation laws numerically. Some of these numerical difficulties are investigated. A new class of numerical schemes is developed to overcome some of these problems. The numerical schemes are constructed using an infinite sequence of conservation laws. Restrictions are given on the schemes that guarantee they maintain a uniform bound and satisfy an entropy condition. For schemes meeting these criteria, a proof is given of convergence to the correct physical solution of the conservation law. Numerical examples are presented to illustrate the theoretical results.

A novel high-order finite volume scheme using flux correction methods in conjunction with structured finite difference schemes is extended to low Mach and incompressible flows on strand grids. Flux correction achieves high-order by explicitly canceling low-order truncation error terms in the finite volume cell. The flux correction method is applied in unstructured layers of the strand grid. The layers are then coupled together using a source term containing the derivatives in the strand direction. Proper source term discretization is verified. Strand-direction derivatives are obtained by using summation-by-parts operators for the first and second derivatives. A preconditioner is used to extend the method to low Mach and incompressible flows. We further extend the method to turbulent flows with the Spalart Allmaras model. We verify high-order accuracy via the method of manufactured solutions, method of exact solutions, and physical problems. Results obtained compare well to analytical solutions, numerical studies, and experimental data. It is foreseen that future application in the Naval field will be possible.

No description available.

No description available.

Ce travail est consacré au développement de schémas numériques pour approcher les solutions de modèles d'écoulement type “depth averaged”. Dans un premier temps, nous détaillons la construction d'approches Volumes Finis pour le système Shallow Water avec termes sources sur maillages non structurés. En se basant sur une reformulation appropriée des équations, nous mettons en place un schéma équilibré et préservant la positivité de la hauteur d'eau, et suggérons des extensions MUSCL adaptées. La méthode est capable de gérer des topographies irrégulières et exhibe de fortes propriétés de stabilité. L'inclusion des termes de friction fait l'objet d'une analyse poussée, aboutissant à l'établissement d'une propriété type “Asymptotic Preserving” à travers l'amélioration d'un autre récent schéma Volumes Finis. La seconde composante de cette étude concerne les méthodes Elements Finis type Galerkin discontinu. Certaines des idées avancées dans le contexte Volumes Finis sont employées pour aborder le système Shallow Water surmaillages triangulaires. Des résultats numériques sont exposés et la méthode se révèle bien adaptée à la description d'une large variété d'écoulements. Partant de ces observations nous proposons finalement d'exploiter ces caractéristiques pour étendre l'approche à une nouvelle famille d'équations type Green-Nadghi. Des validations numériques sont également proposées pour valider le modèle numérique. / This work is devoted to the development of numerical schemes to approximatesolutions of depth averaged flow models. We first detail the construction of Finite Volume approaches for the Shallow Water system with source terms on unstructured meshes. Based on a suitable reformulation of the equations, we implement a well-balanced and positive preserving approach, and suggest adapted MUSCL extensions. The method is shown to handle irregular topography variations and demonstrates strong stabilities properties. The inclusion of friction terms is subject to a thorough analysis, leading to the establishment of some Asymptotic Preserving property through the enhancement of another recent Finite Volume scheme.The second aspect of this study concerns discontinuous Galerkin Finite Elementmethods. Some of the ideas advanced in the Finite Volume context areemployed to broach the Shallow Water system on triangular meshes. Numericalresults are exposed and the method turns out to be well suited to describe a large variety of flows. On these observations we finally propose to exploit its features to extend the approach to a new family of Green-Nadghi equations. Numerical experiments are also proposed to validate this numerical model.