Since the 70s of last century, the finite element method has begun to be applied to the shallow water equations. The wave seems to spread out from the center, but very slowly. Pdf explicit finite difference method has been widely used for seismic modeling in heterogeneous media with. So, we will take the semidiscrete equation 110 as our starting point. Its instructive to look at some important examples to see how they arise. So the general answer to learning finite difference methods is to take a class revolving around numerical analysis, numerical methods, or computational physics. The finite difference method allows you to easily investigate the wavefunction dependence upon the total energy. Find materials for this course in the pages linked along the left. Pdf this paper introduces the fourth order compact finite difference method for solving the numerical. The performance and the accuracy of the method are illustrated by solving three test examples of the problem. Pdf finite difference method for solving acoustic wave. This paper introduces the fourth order compact finite difference method for solving the numerical solution of onedimensional wave equations. Keywords compact finite differences wave equation variable coefficients high order.
Convergence of finite difference methods for the wave. Many facts about waves are not modeled by this simple system, including that wave motion in water can depend on the depth of the. Finitedifference approximation of wave equations acoustic waves in 1d to solve the wave equation, we start with the simplemost wave equation. Numerical experiments corroborate our accuracy analysis. Many facts about waves are not modeled by this simple system, including that wave motion in water can depend on the depth of the medium, that waves tend to disperse, and. Pdf splitting algorithms for the highorder compact. Dec 09, 2014 a simple finite difference scheme for the scalar wave equation on a square mesh with inhomogeneous phase velocity. Equation 8 suggests that the finitedifference scheme for the divergence is of the same secondorder form. Finite di erence methods for wave motion github pages. Splitting algorithms for the highorder compact finitedifference schemes in wave equation modeling. This method discretizes the elastic wave equation as a.
A nite di erence method proceeds by replacing the derivatives in the di erential equations by nite di erence approximations. Finite difference methods, clarendon press, oxford. A linearized implicit finite difference method to obtain numerical solution of the onedimensional regularized longwave rlw equation is presented. Finite difference method for solving acoustic wave. A recipe for stability analysis of finitedifference wave. A fast linearized finite difference method for the nonlinear. Finitedifference modelling of twodimensional elastic wave. Zienkiewicz 34, and peraire 22 are among the authors who have worked on this line. The equation is discretized by a finite difference operator satisfying a summation by parts property, and the boundary and grid interface conditions are imposed weakly by the simultaneous approximation term method. In parallel to this, the use of the finite volume method has grown. The following is my matlab code to simulate a 2d wave equation with a gaussian source at center using fdm. Pdf fourth order compact finite difference method for solving. Finite difference method for the diffusiveviscous wave equation haixia zhao, jinghuai gao, and zhangxin chen abstract.
In the finite difference method, we relax the condition that holds at all points in the spacetime domain \. The method of lines for solution of the onedimensional wave. This manual contains a detailed description of the above problems, the. Convergence of finite difference methods for the wave equation in two space dimensions. In this introductory work i will present the finite difference method for hyperbolic equations, focusing on a method which has second order precision both in time and space the socalled staggered leapfrog method and applying it to the case of the 1d and 2d wave equation. Hans petter langtangen 1, 2 1 center for biomedical computing, simula research laboratory 2 department of informatics, university of oslo. Mitra department of aerospace engineering iowa state university introduction laplace equation is a second order partial differential equation pde that appears in many areas of science an engineering, such as electricity, fluid flow, and steady heat conduction. Fourth order finite difference methods for the wave equation with. We first propose a discretization to the multiterm caputo derivative based on the recently established fast l 21. Schrodinger equation but this is more complex to write the mscript and not as versatile as using the finite difference method.
Finite difference method of wave equation physics forums. The 1d scalar wave equation for waves propagating along the x axis. The finite difference method relies on discretizing a function on a grid. Finite di erence method nonlinear ode heat conduction with radiation if we again consider the heat in a metal bar of length l, but this time consider the e ect of radiation as well as conduction, then the steady state equation has the form u xx du4 u4 b gx. Finite difference methods are quite fundamental when it comes to solving differential. Finite difference method for pde using matlab mfile 23. A recipe for stability analysis of finite difference wave equation computations laurence r. The accuracy analysis is in a general framework, here applied to the second order wave equation. Finite difference method for solving acoustic wave equation. Finite difference method applied to 1d convection in this example, we solve the 1d convection equation. Finite difference method for solving differential equations. A brief derivation of the energy and equation of motion of a wave is done before the numerical. Higher order compact finitedifference method for the wave equation a compact finite difference scheme comprises of adjacent point stencils of which differences are taken at the middle node, therefore typically 3, 9 and 27 nodes are used for compact finite difference descretization in one.
We develop the finitedifference algorithm to the acoustic wave equation in. Derivation of the finitedifference equation 23 following the conventions used in figure 21, the width of cells in the row direction, at a given column, j, is designated. Finite difference methods massachusetts institute of. Since both time and space derivatives are of second order, we use centered di erences to approximate them. To use a finite difference method to approximate the solution to a problem, one must first discretize the problems domain. A simple finite difference scheme for the scalar wave equation on a square mesh with inhomogeneous phase velocity.
Convergence of summationbyparts finite difference methods. Study notes on numerical solutions of the wave equation with the. The heart of the finite difference method is the approximation of the second derivative by the. In this introductory work i will present the finite difference method for hyperbolic equations, focusing on a method which has second order precision. A fourth order accurate finite difference scheme for the. Finitedifference numerical methods of partial differential equations. Solving the heat, laplace and wave equations using. Phillip bording introduction finite difference solutions to the wave equation are pervasive in the modeling of seismic wave propagation kelly and marfurt, 1990 and in seismic imaging bording and lines, 1997. Finite di erence methods for wave motion hans petter langtangen 1. Understand what the finite difference method is and how to use it to solve problems. Sep 27, 2016 in this paper, we develop techniques for analyzing this phenomenon, and particularly consider the second order wave equation. Using excel to implement the finite difference method for 2d heat transfer in a mechanical engineering technology course abstract. Bancroft abstract a new method of migration using the finite element method fem and the finite difference method fdm is jointly used in the spatial domain. The mathematical derivation of the computational algorithm is accompanied by python.
Convergence of finite difference methods for the wave equation in two. This gives a large algebraic system of equations to be solved in place of the di erential equation, something that is easily solved on a computer. Many facts about waves are not modeled by this simple system, including that wave motion in water can depend on the depth of the medium, that. Solving wave equations by the finitedifference timedomain method. In this paper, we study a fast and linearized finite difference method to solve the nonlinear timefractional wave equation with multi fractional orders. Advanced finite difference methods for seismic modeling yang liu 1,2 and mrinal k sen 2 1state key laboratory of petroleum resource and prospecting china university of petroleum, beijing, beijing, 102249, china 2the institute for geophysics, john a. W3v1 wave equation week 03 the finitedifference method 1d. Doing physics with matlab 2 introduction we will use the finite difference time domain fdtd method to find solutions of the most fundamental partial differential equation that describes wave motion, the onedimensional scalar wave equation. Study notes on numerical solutions of the wave equation with. The convergence of the method for the problem under. Finite difference method for pde using matlab mfile.
Finitedifference modelling of twodimensional elastic. A fast linearized finite difference method for the. Finite difference fd schemes transform the partial differential equation into a difference equation by discretizing time and space. We analyze two types of summationbyparts finite difference operators for approximating the second derivative with variable.
Finitedifference method for the diffusiveviscous wave equation haixia zhao, jinghuai gao, and zhangxin chen abstract. Pdf fourth order compact finite difference method for. Introductory finite difference methods for pdes the university of. Finite difference, finite element and finite volume methods for the numerical solution of pdes vrushali a. How to learn finitedifference methods for 2d wave equations.
Later, we use this observation to conclude that bordings conjecture for stability of finite difference schemes for the scalar wave equation lines et al. Finitedifference modelling of twodimensional elastic wave propagation in cracked media. The di usiveviscous wave equation plays an important role in seismic exploration and it can be used to explain the frequencydependent re ections observed both in laboratory and eld data. Finite difference method for the solution of laplace equation. Some numerical results and comparisons with the finite difference method presented in 1 are given in. The straightforward finite difference approximation to the second partial derivative is. Sometimes an analytical approach using the laplace equation to describe the problem can be used. This is usually done by dividing the domain into a uniform grid see image to the right.
Lecture notes numerical methods for partial differential. It has been applied to solve a time relay 2d wave equation. The resulting methods are called finite difference methods. Finite difference modelling of the full acoustic wave. Finite difference method for solving acoustic wave equation using locally adjustable timesteps. Finite di erence schemes and the schrodinger equation. Using excel to implement the finite difference method for. Advanced finitedifference methods for seismic modeling. In order to validate our finite difference method, we apply it to a single crack in a homogeneous embedding and compare its scattered wavefield with the scattered field resulting from a direct evaluation of the integral representation given in eq. Download free books at 4 introductory finite difference methods for pdes contents contents preface 9 1. They have tested the new method for two examples from the literature 18. Jackson school of geosciences, the university of texas at austin, 10100. The idea is to replace the derivatives appearing in the differential equation by finite differences that approximate them. A high order compact timespace finite difference scheme for the.
A finite difference solution of the regularized longwave. Comparison of finite difference schemes for the wave. Solving the heat, laplace and wave equations using nite. Common applications of the finite difference method are in computational science and engineering disciplines, such as thermal engineering, fluid mechanics, etc. Chapter 4 the wave equation another classical example of a hyperbolic pde is a wave equation. In this paper, we develop techniques for analyzing this phenomenon, and particularly consider the second order wave equation. Just a little push of my doubts i hope somebody could help me with my confusion of one of our home works. A finite difference approximation was also defined for the wave equation, and the cfl stability condition was shown to be necessary. Discontinuous galerkin finite element method for the wave. A program for the control and stabilization of waves.
Fourth order compact finite difference method for solving one dimensional wave equation. Finite difference method for solving acoustic wave equation using locally adjustable timesteps alexandre j. The wave equation is a secondorder linear hyperbolic pde that describesthe propagation of a variety of waves, such as sound or water waves. Study notes on numerical solutions of the wave equation. Finite element method 2 acoustic wave equation in 1d how do we solve a timedependent problem such as the acoustic wave equation.
Finite difference methods for waves on a string the complete initialboundary value problem input data in the problem. A numerical scheme for solving the secondorder wave equation with given initial conditions and a. Multidimensional heat transfer problems can be approached in a number of ways. Finite difference methods for ordinary and partial. Finite difference method for the solution of laplace equation ambar k. Solve 1d wave equation using finite difference method. Solving heat equation using finite difference method. Stepwave test for the lax method to solve the advection % equation clear. A linearized implicit finite difference method to obtain numerical solution of the onedimensional regularized long wave rlw equation is presented. Pdf finite difference method for solving acoustic wave equation. Finite di erence schemes and the schrodinger equation jonathan king, pawan dhakal june 2, 2014 1 introduction in this paper, we primarily explore numerical solutions to the quantum 1d in nite square well problem. Sep 26, 2018 so the general answer to learning finite difference methods is to take a class revolving around numerical analysis, numerical methods, or computational physics. Finite difference method for solving acoustic wave equation using locally adjustable timesteps article pdf available in procedia computer science 29.
1178 1493 1323 743 974 1469 31 766 248 167 326 1450 777 1348 969 1658 920 1292 1128 640 63 68 22 1418 1371 1141 1065 1160 40 1022 1319 788 962 1371 741