Solving One-Dimensional Advection-Dispersion with Reaction Using Some Finite-Difference Methods
|
|
- Frank Simmons
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol., 008, no. 53, 6-68 Solving One-Dimensional Advection-Dispersion with Reaction Using Some Finite-Difference Methods H. Saberi Naafi Department of Mathematics, Faculty of Sciences Guilan University, Iran H. Hainezhad Department of Mathematics, Faculty of Sciences Payame Noor University, Iran Haniye Abstract In this paper we have solved the one-dimensional advection-dispersion equation with reaction using some finite-difference methods, namely, FTCS used by B.Ataie-Ashtiani et al. (996), (we call it FTCS-Ataie) [], Upstream, DuFort-Frankel, (3-3), and Crank-Nicolson methods. The implementation and behavior of the methods have been compared by plotting the figures. It has been shown that the Crank-Nicolson method works more accurately than the other methods. The equation is the same as the one used in the article [Journal of Contaminant Hydrology 3 (996) 49-56]. Keywords: Advection, Dispersion, Finite-difference. Introduction Advection-dispersion phenomena occur in many physical situation including the transfer of heat in fluids, flow through porous media, the spread of contaminants in fluids and in chemical separation processes, and so on. The one-dimensional advection-dispersion has been used and solved by researchers, for example, it has been used to describe thermal pollution in a river system (chaudhry et al. 983) leaching of salts in soil (Bresler, 973), flow in porous media (kumar, 983), the spread of pollutants in streams (Cunge et al. 980, Mc Bride and Rutherford, 984), the dispersion of contaminants in estuaries (Winterwerp, 983), and has been solved by researchers like B.J.Noye [4]
2 6 H. Saberi Naafi and H. Hainezhad and others. For the more general transport equation (e. g, with reaction) it has been studied for example by Ataie-Ashtiani et al. [], Khalid Alhumaizi, (Computers & Chemical Engineering. 8(004) ) and Moldrup et al. (994). The aim of this paper is solving the one-dimensional advection-dispersion with reaction, using: FTCS-Ataie, Upstream, DuFort-Frankel, (3-3) and Crank- Nicolson finite difference methods and study the implementation of the methods by comparison the results.. The Mathematical equation The one-dimensional advection-dispersion equation which describes the transient transport of solutes through a homogeneous soil is: t = D C z u kc () z Including first-order reaction, has been taken from [], where C: is solute concentration, (ML 3 ), t: is time, (T ), z: is soil depth, (L), u: is the pore water velocity, (LT ), D: is the dispersion coefficient, (L T ), k: is the first-order reaction rate coefficient, (T ). Since the solute mass leaving from a segment via convection during the time interval t is less than the solute mass inside at the beginning of the time interval. (see [3]), therefore the mass balance criterion of the system is: t u This condition must come with condition stability to show the stability criteria. 3. The finite-difference methods 3.. The FTCS-Ataie method [] The equation () evaluated at the (,n)th grid point, using the FTCS-Ataie, the equation becomes C n+ =( D t () + u t )Cn +( k t D t () )Cn +(D t () u t )Cn + () Where D = D + t (u kd),u = u ( t)ku, k = k t k Ataie et al, in [] investigated the stability of this method using the matrix method and concluded that to satisfy the stability criteria if t D + k
3 Finite-difference methods The Upstream method [4] In this method the backward-difference form uses for the spatial derivative in the term u, forward-difference form for the time derivative and the centered z -difference form for the space derivative, the equation () becomes C n+ C n t = D Cn + Cn +Cn u Cn Cn kc n () This leads to the explicit finite- difference formula C n+ =(s + c)c n +( s c k t)cn + scn + (3) Where s = D t and c = u t () It is easy to show by matrix method [4] and Gerschgorin s circle theorem [] that, (3) is stable if k t < and s + c + k t <. Since k>0 it is sufficient to see t < D () + u + k Note: The equation (3) satisfies the consistency conditions (see []) therefore by Lax theorem [4] in this case with well-posed initial condition the solution which approximates the equation () is convergent The DuFort-Frankel method [4] Another explicit method for solving equation () is the DuFort and Frankel method and We use the following approximation for reaction term C n = (Cn+ + C n ) Using this method the equation () becomes C n+ C n = D Cn + (Cn+ + C n () )+C n u Cn + Cn k Cn+ + C n Therefore C n+ = c +s +s + k t Cn c s s k t +s + k t Cn + + +s + k t Cn (4) Using von Neumann [4] and Miller theorem [4], (4) is stable if G = ((s + k t) c 4s )Sin β +4s (s + k t) Note: The equation (4) satisfies the consistency conditions and with wellposed initial condition the solution which approximates the equation () is convergent.
4 64 H. Saberi Naafi and H. Hainezhad 3.4. The (3-3) method In this method, we use these approximations C = n+ C n t n+ z n+ C n+ z C n+ t = [ Cn Cn = [ Cn + Cn +Cn () C = n+ +C n + Cn+ C n+ ] + Cn+ + Cn+ +C n+ ] () Therefore the equation () at (, n + )th grid point evaluated as t n+ = D C z n+ u z n+ kc n+ or C n+ C n t = D [Cn + Cn + Cn () u [Cn Cn + Cn+ therefore (c + s)c n+ ( + c +s + k t)c n+ + sc+ n+ + Cn+ + C n+ + C n+ ] () C n+ ] k Cn+ + C n = (c + s)c n ( (c +s + k t))cn scn + (5) Using von Neumann method shows that the equation (5) is always stable. If the resulting system of equation solves by Thomas algorithm, it is always stable The Crank-Nicolson method [4] This is an implicit method which uses the same approximations as (3-3) method and z n+ = [Cn + Cn Therefore equation () in this case becomes + Cn+ + C n+ ] (c +s)c n+ +(4+4s +k t)c n+ (s c)c+ n+ =(c +s)c n +(4 4s k t)cn +(s c)cn + (6) The von Neumann method shows that the equation (6) is always stable. If the resulting system of equation solves by Thomas algorithm it is stable if c s.
5 Finite-difference methods Numerical tests To show the implementation and behavior of each method and compare the results obtained from the methods, the analytic solution of the equation () will be used, which is [5] C(z, t) = C 0 [exp((u v)z D z vt )erfc( )+exp( (Dt) (u + v)z D z + vt )erfc( )] (Dt) Where v =(u +4kD) the conditions are C =0, t =0,z >0 C =0, t 0,z C = C 0, t 0,z =0 We assume C 0 = mgl, D =.0cm, u =0.5cmh, k =0., 0.h, =.0cm and =.0h The error function is Error = C Num C Ana where C Num and C Ana respectively are numerical and analytical solution. The results are shown in figures -4.The values of k, t and z used at each test are also shown. 5.Summary and Conclusions The finite difference methods applied to the equation () are i) The FTCS-Ataie, ii) The Upstream, iii)the DuFort-Frankel, iv) The (3-3), v) TheCrank-Nicolson. The results obtained from the methods are plotted in the Figures -4 and describe the following items: a)in the DuFort-Frankel method although the centered-difference formula has been used but as Figures and 3 show the error increases fast, because of the coefficients in (4). b) The FTCS-Ataie works better than Upstream method and (3-3) method, because the approximations in FTCS-Ataie are more accurate, see Figures and 3. c) The Crank-Nicolson method works more accurate than FTCS-Ataie, see Figures and 4. Therefore we suggest using the Crank-Nicolson method for solving the equation ().
6 66 H. Saberi Naafi and H. Hainezhad References [] B.Ataie-Ashtiani, D.A. Lokhington, R.E. Volker, Numerical correction for finite-difference solution of the advection-dispersion equation with reaction, J.Contaminant Hydrology, 996. [] B.N. Datta, Numerical Linear Algebra and Aplications, Brooks/Cole Publishing company, 995. [3] W.Kinzelbach, Groundwater Modeling: An Introduction with Sample Pograms in Basic, Elsevier Amsterdam,986,-65 [4] B.J.Noye, Finite difference Method For Partial Differential Equation, Applied Mathematics Leacture Notes, Department of Applied Mathematics, The university of Adelaide, Australia., 993. [4] M.Th.Van Genuchten and W.J.Alves, Analytical solution of the one-dimention convective-dispersive solute transport equation. U.S. Dep. Agric.,Tech. Bull. No. 66, 98. Received: March 9, 008
7 Finite-difference methods 67 Figure : Shows that errors of the five methods, with k=0., z=0 and t=0()30. The minimum error belongs to the FTCS-Ataie method and Crank- Nicolson method. Figure : shows that errors of Crank-Nicolson and FTCS-Ataie methods with k=0., t=00 and z=0()0. The Crank-Nicolson has the minimum error.
8 68 H. Saberi Naafi and H. Hainezhad Figure 3: Shows that errors of the five methods, with k=0. z=0 and t=0()5. The minimum error belongs to the FTCS-Ataie method and Crank- Nicolson method. Figure 4: shows that errors of Crank-Nicolson and FTCS-Ataie methods with k=0., t=00 and z=0()0. The Crank-Nicolson has the minimum error.
Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
Australian Journal of Basic and Applied Sciences, 5(6): 536-543, 0 ISSN 99-878 Numerical Solution of One-dimensional Advection-diffusion Equation Using Simultaneously Temporal and Spatial Weighted Parameters
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationA semi-lagrangian scheme for advection-diffusion equation
EPiC Series in Engineering Volume 3, 2018, Pages 162 172 HIC 2018. 13th International Conference on Hydroinformatics Engineering A semi-lagrangian scheme for advection-diffusion equation Ersin Bahar 1,
More informationPollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationNUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 130 137 NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
More information1906 Dilip kumar Jaiswal et al./ Elixir Pollution 31 (2011) Available online at (Elixir International Journal)
196 Dilip kumar Jaiswal et al./ Eliir Pollution 31 (211) 196-191 ARTICLE INF O Article history: Received: 21 December 21; Received in revised form: 16 January 211; Accepted: 1 February 211; Keywords Advection,
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationQUESTIONS FOR THE TESTS-ECI 146
QUESTIONS FOR THE TESTS-ECI 46 IMPORTANT: Please note that the questions in the exams WILL BE DIFFERENT than these. (These questions are valid for the midterm and for the final exam) Remember: A) both
More informationIOSR Journal of Mathematics (IOSRJM) ISSN: Volume 2, Issue 1 (July-Aug 2012), PP
IOSR Journal of Mathematics (IOSRJM) ISSN: 78-578 Volume, Issue (July-Aug ), PP - Analytical Solutions of One-Dimensional Temporally Dependent Advection-Diffusion Equation along Longitudinal Semi-Infinite
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationA note on benchmarking of numerical models for density dependent flow in porous media
Advances in Water Resources 29 (2006) 1918 1923 www.elsevier.com/locate/advwatres A note on benchmarking of numerical models for density dependent flow in porous media B. Ataie-Ashtiani *, M.M. Aghayi
More informationSummer Intern Project/Model at Stroud Water Research Center. Xiaojuan (Cathy) Yu
Summer Intern Project/Model at Stroud Water Research Center Xiaojuan (Cathy) Yu If we knew what we were doing, it wouldn t be called Research. A. Einstein Outline Introducing OTIS Mathematical model in
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationRADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA
GEOPHYSICS RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA C. BUCUR 1, M. OLTEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, Pitesti, Romania, crina.bucur@scn.ro 2 Academy of Scientists Bucharest,
More information1 of 7 2/8/2016 4:17 PM
1 of 7 2/8/2016 4:17 PM Continuous issue-4 January - March 2015 Abstract The similarity solution of Diffusion equation using a technique of infinitesimal transformations of groups The similarity solution
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationOn the advective-diffusive transport in porous media in the presence of time-dependent velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 3, L355, doi:.29/24gl9646, 24 On the advective-diffusive transport in porous media in the presence of time-dependent velocities A. P. S. Selvadurai Department of Civil
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More information1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationComputer Aided Design of Thermal Systems (ME648)
Computer Aided Design of Thermal Systems (ME648) PG/Open Elective Credits: 3-0-0-9 Updated Syallabus: Introduction. Basic Considerations in Design. Modelling of Thermal Systems. Numerical Modelling and
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationTEMPORALLY DEPENDENT DISPERSION THROUGH SEMI-INFINITE HOMOGENEOUS POROUS MEDIA: AN ANALYTICAL SOLUTION
IJRRAS 6 () February www.arpapress.com/volumes/vol6issue/ijrras_6 5.pdf TEMPORALLY EPENENT ISPERSION THROUGH SEMI-INFINITE HOMOGENEOUS POROUS MEIA: AN ANALYTICAL SOLUTION R.R.Yadav, ilip Kumar Jaiswal,
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationCHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH
SCA006-49 /6 CHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH Mostafa FOURAR LEMTA Ecole des Mines de Nancy, Parc de Saurupt, 54 04 Nancy, France
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationCoupling TRIGRS and TOPMODEL in shallow landslide prediction. 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 2016/10/13
Coupling TRIGRS and TOPMODEL in shallow landslide prediction 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 016/10/13 Outline Introduction Literature review Methodology Pre-result Future work Introduction 3 Motivation
More informationStabilitätsanalyse des Runge-Kutta- Zeitintegrationsschemas für das konvektionserlaubende Modell COSMO-DE
Stabilitätsanalyse des Runge-Kutta- Zeitintegrationsschemas für das konvektionserlaubende Modell COSMO-DE DACH-Tagung 20.-24.09.200, Bonn Michael Baldauf GB FE, Deutscher Wetterdienst, Offenbach M. Baldauf
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
More informationConditional stability of Larkin methods with non-uniform grids
Theoret. Appl. Mech., Vol.37, No., pp.139-159, Belgrade 010 Conditional stability of Larkin methods with non-uniform grids Kazuhiro Fukuyo Abstract Stability analysis based on the von Neumann method showed
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationRT3D Rate-Limited Sorption Reaction
GMS TUTORIALS RT3D Rate-Limited Sorption Reaction This tutorial illustrates the steps involved in using GMS and RT3D to model sorption reactions under mass-transfer limited conditions. The flow model used
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale Ekkehard Holzbecher, Sandra Oehlmann October 10 th, 2012 Excerpt from the Proceedings of the 2012 COMSOL Conference in Milan Heat & Mass Transfer
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationAdvanced numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
More informationA breakthrough curve analysis of unstable density-driven flow and transport in homogeneous porous media
WATER RESOURCES RESEARCH, VOL. 40,, doi:10.1029/2003wr002668, 2004 A breakthrough curve analysis of unstable density-driven flow and transport in homogeneous porous media M. Wood, C. T. Simmons, and J.
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More informationThermohaline Flow and Reactive Solute Transport in Fractured Rock
Thermohaline Flow and Reactive Solute Transport in Fractured Rock Thomas GRAF Institute of Fluid Mechanics in Civil Engineering Alberta Geological Survey, Edmonton 2008-2009 Basin Hydrogeologist Sasdf
More informationDeutscher Wetterdienst
Stability Analysis of the Runge-Kutta Time Integration Schemes for the Convection resolving Model COSMO-DE (LMK) COSMO User Seminar, Langen 03.+04. March 2008 Michael Baldauf Deutscher Wetterdienst, Offenbach,
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More information1D Heat equation and a finite-difference solver
1D Heat equation and a finite-difference solver Guillaume Riflet MARETEC IST 1 The advection-diffusion equation The original concept applied to a property within a control volume from which is derived
More informationCOMPLETED RICHARDSON EXTRAPOLATION IN SPACE AND TIME
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING, Vol. 13, 573±582 (1997) COMPLETED RICHARDSON EXTRAPOLATION IN SPACE AND TIME SHANE A. RICHARDS Department of Applied Mathematics, The University of Adelaide,
More informationNATURAL ZEOLITE AS A PERMEABLE REACTIVE BARRIER
NATURAL ZEOLITE AS A PERMEABLE REACTIVE BARRIER PREDICTION OF LEAD CONCENTRATION PROFILE THROUGH ZEOLITE BARRIER N. Vukojević Medvidović, J. Perić, M. Trgo, M. Ugrina, I. Nuić University of Split, Faculty
More informationA three-point time discretization technique for parabolic partial differential equations
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech. 2011; 35:406 418 Published online 23 March 2010 in Wiley Online Library (wileyonlinelibrary.com)..902
More informationInternational Research Journal of Engineering and Technology (IRJET) e-issn: Volume: 04 Issue: 09 Sep p-issn:
International Research Journal of Engineering and Technology (IRJET) e-issn: 395-56 Volume: 4 Issue: 9 Sep -7 www.irjet.net p-issn: 395-7 An analytical approach for one-dimensional advection-diffusion
More informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More informationJournal of Applied Mathematics and Computation (JAMC), 2018, 2(3), 67-83
Journal of Applied Mathematics and Computation (JAMC) 8 (3) 67-83 http://www.hillpublisher.org/journal/jamc ISSN Online:576-645 ISSN Print:576-653 Analysis of two-dimensional solute transport through heterogeneous
More informationAnalytical approach predicting water bidirectional transfers: application to micro and furrow irrigation
Advances in Fluid Mechanics VI 633 Analytical approach predicting water bidirectional transfers: application to micro and furrow irrigation D. Crevoisier Irrigation Research Unit, Cemagref Montpellier,
More informationThe generation and mobility of colloids in soils Production et mobilité des colloï des dans les sols
Scientific registration n o : 1112 Symposium n o : 4 Presentation : poster The generation and mobility of colloids in soils Production et mobilité des colloï des dans les sols WATERS Angela, CHITTLEBOROUGH
More informationO.R. Jimoh, M.Tech. Department of Mathematics/Statistics, Federal University of Technology, PMB 65, Minna, Niger State, Nigeria.
Comparative Analysis of a Non-Reactive Contaminant Flow Problem for Constant Initial Concentration in Two Dimensions by Homotopy-Perturbation and Variational Iteration Methods OR Jimoh, MTech Department
More informationThe distortion observed in the bottom channel of Figure 1 can be predicted from the full transport equation, C t + u C. y D C. z, (1) x D C.
1 8. Shear Dispersion. The transport models and concentration field solutions developed in previous sections assume that currents are spatially uniform, i.e. u f(,y,). However, spatial gradients of velocity,
More informationSunday September 28th. Time. 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn. CP: Contributed Presentation 25 mn
Sunday September 28th 06:00 pm 09:00 pm Registration (CAAS) IP: Invited Presentation (IS or SP) 55 mn CP: Contributed Presentation 25 mn Monday September 29th 07:30 am 08:30 am Registration 08:30 am 09:00
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationSimulation of seawater intrusion in coastal aquifers: Some typical responses
Sādhanā, Vol. 26, Part 4, August 2001, pp. 317 352. Printed in India Simulation of seawater intrusion in coastal aquifers: Some typical responses 1. Introduction AMLAN DAS 1 and BITHIN DATTA 2 1 Department
More informationDissipation and Dispersion
Consider the problem with periodic boundary conditions Dissipation and Dispersion u t = au x 0 < x < 1, t > 0 u 0 = sin 40 πx u(0, t) = u(1, t) t > 0 If a > 0 then the wave is moving to the left and if
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationABSTRACT INTRODUCTION
Transport of contaminants from non-aqueous phase liquid pool dissolution in subsurface formations C.V. Chrysikopoulos Department of Civil Engineering, University of California, ABSTRACT The transient contaminant
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationAccuracy of Muskingum-Cunge flood routing
Alternative Hydraulics Paper 3, 02.03.2011 Accuracy of Musingum-Cunge flood routing Institute of Hydraulic and Water Resources Engineering Vienna University of Technology Karlsplatz 13/222, 1040 Vienna,
More informationChemical Hydrogeology
Physical hydrogeology: study of movement and occurrence of groundwater Chemical hydrogeology: study of chemical constituents in groundwater Chemical Hydrogeology Relevant courses General geochemistry [Donahoe]
More informationFrozen Ground Containment Barrier
Frozen Ground Containment Barrier GEO-SLOPE International Ltd. www.geo-slope.com 1200, 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: +1 403 269 2002 Fax: +1 888 463 2239 Introduction Frozen soil barriers
More informationABSTRACT GOVERNING EQUATIONS
A three dimensional finite element model for fluid flow and transport in confined or unconfined aquifer F. Jacob*, J.M. Crolef, P. Lesaint*, J. Mania* "Laboratoire de calcul scientifique, ^Laboratoire
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationAn Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD
An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface
More informationModeling the Fate and Transport of Engineered Nanoparticles in Porous Media
Modeling the Fate and Transport of Engineered Nanoparticles in Porous Media Nadim Copty Institute of Environmental Sciences Bogazici University Istanbul, Turkey NM2PorousMedia, Dubrovnik, Croatia, Sept.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationEnvironmental Hydraulics of Open Channel Flows
Environmental Hydraulics of Open Channel Flows Hubert Chanson ME, ENSHM Grenoble, INSTN, PhD (Cant), DEng (Qld) Eur Ing, MIEAust, MIAHR 13th Arthur Ippen awardee (IAHR) Reader in Environmental Fluid Mechanics
More informationNumerical Solution of Burger s_fisher Equation in One - Dimensional Using Finite Differences Methods
Fluid Mechanics 2018; 41): 20-26 http://www.sciencepublishinggroup.com/j/fm doi: 10.11648/j.fm.20180401.13 ISSN: 2575-1808 Print); ISSN: 2575-1816 Online) Numerical Solution of Burger s_fisher Equation
More information