The difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method. Harbaugh abstract this report presents a finite difference model and its associated. Finitevolume transport on various cubedsphere grids william m. Expressions are given for the first and second derivatives at a grid point in terns of nodal values. I we therefore consider some arbitrary function fx, and suppose we can evaluate it at the uniformly spaced grid points x1,2 3, etc. Mimetic finite difference discretizations on triangular grids. Mimetic finite difference methods for diffusion equations. Finite differences on triangular grids brighi 1998 numerical. Mimetic finite difference methods for diffusion equations on. Finite difference schemes university of manchester. There are some analyses on the multigrid methods for the cell centered finite difference on rectangular meshes 5, 7, 8, where the trivial injection was chosen as the prolongation. An optimized variablegrid finitedifference method for seismic forward modeling chunling wu and jerry m. Comparison of vcycle multigrid method for cellcentered.

A triangular grid finitedifference model for windinduced circulation in shallow lakes. Ganzha, analysis and optimization of inner products for mimetic finite difference methods on triangular grid, mathematics and computers in simulation, 67 2004, pp. Finite differences on triangular grids, numerical methods. Pdf conservative finitedifference methods on general grids. Ifnand g are the number of grid cells and the grid level, respectively, we can write n553223. Abstract we define a finite differences method for triangular grids and we show how to link it to a finite element method. I large grid distortions need to be avoided, and the schemes cannot easily be applied to very complex ow geometry shapes. A secondorder finite difference discretization of 5. In section 3, we propose a new prolongation operator and compare it with other prolongations for smooth and nonsmooth problems. For example, the simplest secondorder finite difference approximation to the gradient, on a square grid, uses only wall neighbors.

Finitedifference solution to the eikonal equation for. A triangular grid finitedifference model for windinduced. It can be used directly for forward modeling on models with complex top surfaces and migration without statics preprocessing. The support operator method first constructs a discrete divergence operator from the divergence theorem and then constructs a discrete gradient operator as the adjoint operator of. In each case, the finitedifference scheme on the triangular grid reduces the anisotropy of. If you have a triangular mesh with data defined at each vertex, you can of course associate a linear function with each element. For the purpose of comparing our results with several different triangular cavity studies with different triangle geometries, a general triangle mapped onto a computational domain is considered. Development of irregulargrid finite difference method ifdm. Triangular grid used for finite difference algorithm. For example, the simplest secondorder finitedifference approximation to the gradient, on a square grid, uses only wall neighbors.

Regular article finitedifference schemes on regular triangular grids. An unstructured grid, finitevolume, threedimensional. The multigrid technique has been applled to triangular. The section 4 dealing with discretization starts b y describing the t yp e of discretization scalar and v ector functions on unstructured triangular grid, con tin ues b y in tro ducing natural and formal scalar pro ducts of grid functions ends with deriv ation of discrete appro. It is certainly possible to construct finite difference operators on square grids and triangular grids in which information from all neighboring cells is used e. As if it were essentially a finite difference problem, namely, instead of the finite element problem that it only appears to be. In section 4, we test various problems and report the eigenvalues, condition numbers, and contraction.

Generation of finite difference formulas on arbitrarily spaced grids by bengt fornberg abstract. Although the approaches used by these pioneers are different, they share one essential characteristic. Finitedifference mesh aim to approximate the values of the continuous function ft, s on a set of discrete points in t, s plane divide the saxis into equally spaced nodes at distance. Introductory finite difference methods for pdes contents contents preface 9 1.

The finitedifference schemes presented on the triangular grid include a secondorder method, a compact fourthorder method, and a modified compact method designed to extend the accurate wave number range of the numerical approximation. What is the difference between finite difference methods. If nh 1, then the number of internal grid points is n 12. To use a finite difference method to approximate the solution to a problem, one must first discretize the problems domain. A new and efficient generalized finite difference method, irregular grid finite difference method ifdm, developed for the purpose of generality and efficiency. Mesh and finite volumes the mesh can consist of di erent types of cells, e. An optimized variablegrid finitedifference method for. 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. An optimized fourthorder staggered grid finite difference fd operator is derived on a mesh with variable grid. The method was called the finite difference method based on variation principle, which was another independent invention of the finite element method. Rectangular plate sides ad and bc, simply supported sides ab and dc cantilever supported sides and plate is. A finite difference algorithm for solution of stationary diffusion equation on unstructured triangular grid has been developed earlier by a support operator method. Finite volume solution of the twodimensional euler.

Despite not being generally used in industrial codes, finite difference schemes. Comparison of vcycle multigrid method for cellcentered finite difference on triangular meshes do y. From this new point of view we then analyze properties of the solution and. Cellcentered finite volume method let tbe a triangular or cartesian grid of. Finitedifference schemes on regular triangular grids. Finite volume solution of the twodimensional euler equations. Finite difference operators on unstructured triangular meshes. Finite difference method for the solution of laplace equation. The method was called the finite difference method based on variation principle.

A modular threedimensional finite difference groundwater flow model by michael g. Simple recursions are derived for calculating the weights in compact finite difference formulas for any order of derivative and to any order of accuracy on onedimensional grids with arbitrary spacing. A study of edge effects in triangular grid taperedslot. Discretization formulas for unstructured grids ntrs nasa. Galerkin finite difference laplacian operators on isolated. The navierstokes equations in general curvilinear coordinates in streamfunction and vorticity.

Finite differences on triangular grids finite differences on triangular grids brighi, bernard. Analysis of rectangular thin plates by using finite. However, to that end, we must look at the problem from a different, or should i rather say a difference perspective. Poststack reversetime migration using a finite difference. Pdf finite difference schemes on triangular cellcentered. A triangular grid nitedi erence model for windinduced circulation in shallow lakes david john mcinerney, hons. Finite difference schemes for the 2d wave equation operating on hexagonal grids and the accompanying. The key to making a finite difference scheme work on an irregular geometry is to have a shape matrix with values that denote points outside, inside, and on the boundary of. Finite differences on triangular grids, numerical methods for. We call the spherical triangular grid generated through this process the raw grid. Jul 20, 2008 the difference method based on triangular grids maintains the simplicity of the difference method and the precision of the finite element method.

Let the square region r, 0 x 1, 0 y 1, be covered by a grid with sides parallel to the coordinate axis and grid spacings such that llx ily h. The cell centered finite difference is a nonconforming example. Thus a finite difference solution basically involves three steps. Mimetic finite differences, finite element methods. A new and efficient generalized finite difference method, irregulargrid finite difference method ifdm, developed for the purpose of generality and efficiency.

An optimized variable grid finite difference method for seismic forward modeling chunling wu and jerry m. Harbaugh abstract this report presents a finitedifference model and its associated. The finite difference equation at the grid point involves five grid points in a fivepoint stencil. Finite difference schemes 201011 5 35 i many problems involve rather more complex expressions than simply derivatives of fitself. These finite difference approximations are algebraic in form. Finite differences on domains with irregular boundaries. In this paper, finite difference method fdm was used to obtain solutions for analysis of thin rectangular flat plates carrying distributed load with the following boundary conditions. Finite difference method for the solution of laplace equation ambar k. Rectangular grid, b skew grid, c triangular grid, and d circular grid. A simple finite difference approach using unstructured. Triangular grid used for finitedifference algorithm. Similar to the finiteelement method in the discretization of a numerical mesh, the grid method is flexible in incorporating surface topography and curved interfaces.

This method can be characterized as a finite volume method having piecewise constant basis functions. A recipe for stability analysis of finitedifference wave. An optimized fourthorder staggeredgrid finitedifference fd operator is derived on a. Finite difference schemes on triangular cellcentered grids. Finite difference schemes on triangular cellcentered. Order of accuracy in a rectangular grid is influenced by. In this paper, with new techniques we give completely zerosum flow numbers for certain classes of triangular grid graphs, namely, regular triangular grids, triangular belts, fans, and wheels.

Multigrid algorithm for cell centered finite difference on. Finite difference schemes 201011 2 35 i finite difference schemes can generally be applied to regularshaped domains using bodytted grids curved grid lines, following domain boundaries. Development of irregulargrid finite difference method. These considerations have led to research on zonal grid ding i, global cartesian gridding, 2 and embedded grids 3 just to name a few. Phillip bording introduction finitedifference 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. A finite element framework for some mimetic finite difference. By careful construction of the dissipative terms, the scheme is designed to be second order accurate in space, provided the grid is smooth, except in the vicinity of shocks, where it behaves as first order accurate. Lee2 1department of mathematics, kaist, taejon 305701, korea 2department of mathematics, seoul national university, seoul 151747, korea received 7. Finite difference methods on rectangular grids are widely used in numerical models of environmental flows. Free printable triangular graph paper created date. Finitevolume transport on various cubedsphere grids. Thesis submitted for the degree of doctor of philosophy in applied mathematics at the university of adelaide faculty of engineering, computer and mathematical sciences school of mathematical sciences. The ifdm is derived in light of greengauss theorem rather than taylorseries expansion theorem that is usually used in any other finite difference method.

Introduction in developing unstructured finite difference equations, jameson if, eq. Norgren 2010, a study of edge effects in triangular grid tapered. On the supraconvergence of elliptic finite difference schemes core. But it is possible to use grids with hanging nodes also, which can be advantageous when doing grid.

To find a numerical solution to equation 1 with finite difference methods, we first need to define a set of grid points in the domaindas follows. Usually the mesh is taken to be conforming in the sense that there are no hanging nodes. In developing unstructured finite difference equations, iameson 1987 eq. The evolution and application of the finite difference. The finite difference schemes presented on the triangular grid include a secondorder method, a compact fourthorder method, and a modified compact method designed to extend the accurate wave number range of the numerical approximation. Generation of finite difference formulas on arbitrarily. Numerical solutions of 2d steady incompressible flow inside a triangular cavity are presented. The center is called the master grid point, where the finite difference equation is used to approximate the pde. Forexample,grids0,1, and2have12,42,and162vertices,respectively. All of the schemes considered are centered and hence nondissipative. Putman a, shianjiann lin b a nasa gsfc software integration and visualization o. The finite difference method relies on discretizing a function on a grid.

It is certainly possible to construct finitedifference operators on square grids and triangular grids in which information from all neighboring cells is used e. Fdem is a blackbox solver that solves by a finite difference method. Regarding grid uniformity, although far superior than the latitudelongitude grid, the cubedsphere is not. I am aware how to do this when dealing with finite difference methods on the standard orthogonal grid, using the standard formulas. S apart, and, the taxis into equally spaced nodes a distance. A recipe for stability analysis of finitedifference wave equation computations laurence r. It turns out that the finite difference equation is a certain nonstandard finite element scheme on triangular grids combined with a special form of quadrature. We present a new numerical modelling algorithm for psvwave propagation in heterogeneous media, which is named the grid method in this paper. But it is an excellent compromise when considering grid uniformity, orthogonality, and ease of applying existing highorder.

1454 588 904 1054 84 1074 80 520 1005 187 139 1185 1458 809 906 180 1115 1132 1044 1376 1232 422 923 288 146 1495 1512 487 797 1179 1388 59 516 883 1437 757 1380