Something Old and Something New: An Assessment of Two Numerical Solution Techniques for a Parabolic Partial Differential Equation Arising in Turbulent Particle Transport

Presentation Type

Oral and/or Visual Presentation

Presenter Major(s)

Mathematics

Mentor Information

James McNair

Department

Annis Water Resource Institute (AWRI)

Location

Kirkhof Center 2266

Start Date

11-4-2012 12:30 PM

Keywords

Life Science, Mathematical Science

Abstract

Turbulent particle transport is a biologically important physical process in stream ecosystems. One of the models employed in studying this process is called the Local Exchange Model (LEM), which is a second-order parabolic partial differential equation with boundary conditions at the water surface and stream bed. A widely used numerical method known as the Crank-Nicolson scheme was used to obtain numerical solutions of the LEM. Though this method is commonly claimed to be unconditionally stable, it produced spurious oscillations for many parameter values when applied to the LEM. These oscillations were caused by the boundary conditions, which are not considered in the usual proof of stability for Crank-Nicolson. A review of the recent literature uncovered a new numerical method called exponential fitting, which was specifically designed to avoid spurious oscillations. This method successfully eliminated the oscillation problem when applied to the LEM.

This document is currently not available here.

Share

COinS
 
Apr 11th, 12:30 PM

Something Old and Something New: An Assessment of Two Numerical Solution Techniques for a Parabolic Partial Differential Equation Arising in Turbulent Particle Transport

Kirkhof Center 2266

Turbulent particle transport is a biologically important physical process in stream ecosystems. One of the models employed in studying this process is called the Local Exchange Model (LEM), which is a second-order parabolic partial differential equation with boundary conditions at the water surface and stream bed. A widely used numerical method known as the Crank-Nicolson scheme was used to obtain numerical solutions of the LEM. Though this method is commonly claimed to be unconditionally stable, it produced spurious oscillations for many parameter values when applied to the LEM. These oscillations were caused by the boundary conditions, which are not considered in the usual proof of stability for Crank-Nicolson. A review of the recent literature uncovered a new numerical method called exponential fitting, which was specifically designed to avoid spurious oscillations. This method successfully eliminated the oscillation problem when applied to the LEM.