tailieunhanh - shallow liquid simulation using matlab 2001 neumann

AMATH 581 Homework 2 Shallow Liquid Simulation Erik Neumann 610 N. 65th St., Seattle, WA 98103 erikn@ November 19, 2001 Abstract A model of shallow fluid behavior is evaluated using a variety of numerical solving techniques. The model is defined by a pair of partial differential equations which have two dimensions in space and one dimension of time. The equations concern the vorticity ω and the stream function ψ which are related to the velocity field of the fluid. The equations are first discretized in time and space. The time behavior is evaluated using a Runge-Kutta ordinary differential equation solver. The spatial. | AMATH 581 Homework 2 Shallow Liquid Simulation Erik Neumann 610 N. 65th St. Seattle WA 98103 erikn@ November 19 2001 Abstract A model of shallow fluid behavior is evaluated using a variety of numerical solving techniques. The model is defined by a pair of partial differential equations which have two dimensions in space and one dimension of time. The equations concern the vorticity w and the stream function b which are related to the velocity field of the fluid. The equations are first discretized in time and space. The time behavior is evaluated using a Runge-Kutta ordinary differential equation solver. The spatial behavior is solved using either Fast Fourier Transform Gaussian Elimination LU Decomposition or iterative solvers. The performance of these techniques is compared in regards to execution time and accuracy. Contents 1 Introduction and Overview 2 2 Theoretical Background 4 Solving for b - Matrix Method. 5 Solving for b - FFT Method. 6 Discretize the Advection-Diffusion Equation. 6 3 Algorithm Implementation and Development 7 Construction of Matrix A. 8 Construction of Matrix B. 9 Pinning the Value of b 1 1 . 10 1 Comparing Solvers. 11 An FFT problem. 11 4 Computational Results 11 Results for various initial conditions . 11 Running times . 17 Accuracy of solvers. 17 Symmetry of Solution. 18 Time Resolution Needed. 18 Mesh Drift Instability. 18 5 Summary and Conclusions 21 A MATLAB functions used 22 B MATLAB code 23 . 23 . 23 . 23 . 23 . 24 C Calculations 34 1 Introduction and Overview We consider the governing equations associated with shallow fluid modeling. The intended application is the flow of the earth s atmosphere or ocean circulation. The model assumes a 2-dimensional flow with not much movement up or down. Another assumption is that the fluid is shallow ie. that the vertical dimension is much smaller than the horizontal .