High Resolution Schemes for Conservation Laws With Source Terms.

Author

Martínez i Gavara, Anna

Director

Donat Beneito, Rosa M.

Date of defense

2008-10-24

ISBN

9788437074283

Legal Deposit

V-1287-2009



Department/Institute

Universitat de València. Departament de Matemàtica Aplicada

Abstract

This memoir is devoted to the study of the numerical treatment of<br/>source terms in hyperbolic conservation laws and systems. In particular,<br/>we study two types of situations that are particularly delicate from<br/>the point of view of their numerical approximation: The case of balance<br/>laws, with the shallow water system as the main example, and the case of<br/>hyperbolic equations with stiff source terms.<br/>In this work, we concentrate on the theoretical foundations of highresolution<br/>total variation diminishing (TVD) schemes for homogeneous<br/>scalar conservation laws, firmly established. We analyze the properties<br/>of a second order, flux-limited version of the Lax-Wendroff scheme which<br/>avoids oscillations around discontinuities, while preserving steady states.<br/>When applied to homogeneous conservation laws, TVD schemes prevent<br/>an increase in the total variation of the numerical solution, hence guaranteeing<br/>the absence of numerically generated oscillations. They are successfully<br/>implemented in the form of flux-limiters or slope limiters for<br/>scalar conservation laws and systems. Our technique is based on a flux<br/>limiting procedure applied only to those terms related to the physical<br/>flow derivative/Jacobian. We also extend the technique developed by Chiavassa<br/>and Donat to hyperbolic conservation laws with source terms and<br/>apply the multilevel technique to the shallow water system.<br/>With respect to the numerical treatment of stiff source terms, we take<br/>the simple model problem considered by LeVeque and Yee. We study<br/>the properties of the numerical solution obtained with different numerical<br/>techniques. We are able to identify the delay factor, which is responsible<br/>for the anomalous speed of propagation of the numerical solution<br/>on coarse grids. The delay is due to the introduction of non equilibrium values through numerical dissipation, and can only be controlled<br/>by adequately reducing the spatial resolution of the simulation.<br/>Explicit schemes suffer from the same numerical pathology, even after reducing<br/>the time step so that the stability requirements imposed by the<br/>fastest scales are satisfied. We study the behavior of Implicit-Explicit<br/>(IMEX) numerical techniques, as a tool to obtain high resolution simulations<br/>that incorporate the stiff source term in an implicit, systematic,<br/>manner.

Subjects

51 - Mathematics

Knowledge Area

Facultat de Matemàtiques

Documents

martinez.pdf

4.150Mb

 

Rights

ADVERTIMENT. L'accés als continguts d'aquesta tesi doctoral i la seva utilització ha de respectar els drets de la persona autora. Pot ser utilitzada per a consulta o estudi personal, així com en activitats o materials d'investigació i docència en els termes establerts a l'art. 32 del Text Refós de la Llei de Propietat Intel·lectual (RDL 1/1996). Per altres utilitzacions es requereix l'autorització prèvia i expressa de la persona autora. En qualsevol cas, en la utilització dels seus continguts caldrà indicar de forma clara el nom i cognoms de la persona autora i el títol de la tesi doctoral. No s'autoritza la seva reproducció o altres formes d'explotació efectuades amb finalitats de lucre ni la seva comunicació pública des d'un lloc aliè al servei TDX. Tampoc s'autoritza la presentació del seu contingut en una finestra o marc aliè a TDX (framing). Aquesta reserva de drets afecta tant als continguts de la tesi com als seus resums i índexs.

This item appears in the following Collection(s)