Boundary conditions in computational fluid dynamics wikipedia. Because of the staggered grid configuration, we only need to impose surface boundary conditions for the second and fourthorder moments. Boundary conditions when solving the navierstokes equation and continuity equation, appropriate initial conditions and boundary conditions need to be applied. Simulating seismic wave propagation in 3d elastic media using staggeredgrid finite differences by robert w. Finally, we write gridpoint contributions of the new staggered grid stencil for fourthorder derivative approximations that are valid for the whole medium, including pml absorbing boundary conditions.
The point to note is, that a change within the setting even in this trivial example changes a signi. Furthermore these conditions effectively introduce a dirichlet condition for the pressure. However, staggered grid allows for very natural and accurate formulation of several crucial partial differential equations such as stokes and continuity equations with finite differences. Using staggered grids with characteristic boundary conditions when solving compressible reactive navierstokes equations. A spongelike perfectly matched layer pml absorbing boundary condition. Mar 02, 2016 a staggered grid is a setting for the spatial discretisation, in which the variables are not defined at the same position as in case of the collocated grid. These tests showed that in led driven cavity, for low reynolds numbers the results of staggered and collocated grid methods are the same. Very efficient on cartesian grids are also spectral schemes with an approximate solution being defined as linear combination of basis functions that are nonzero in the whole computational domain. Pdf mixedgrid and staggeredgrid finitedifference methods for. Using staggered grids with characteristic boundary conditions when solving compressible reactive navierstokes equations eric albin, yves dangelo, luc vervisch to cite this version. Among them, the simplest method requires a large enough compu. Details regarding the boundary condition inputs you must supply and the internal treatment at boundaries are provided.
Ive read some literature on it and i still dont quite understand if the pressure boundary conditions are treated differently. Staggeredgrid finitedifference acoustic modeling with. A compact and fast matlab code solving the incompressible. Boundary treatment for fourthorder staggered mesh discretizations. Our boundary conditions on all four sides are when solving the pressure poisson equation on the 4 boundaries is the procedure the same as if we were using the staggered grid. Pdf a rotated staggered grid finitedifference with the absorbing.
Staggered grid combines several types of nodal points located in different geometrical positions and looks rather complex and nonlogical. Improvement in the fault boundary conditions for a staggered. The staggered grid finitedifference method is a powerful tool in seismology and is commonly used to study earthquake source dynamics. We note that the staggeredgrid approach results in. We note that the staggered grid approach results in a stencil that incorporates grid nodes see fig. The solution of viscous incompressible jet flows using non. Implementation of freesurface condition for finite. As i am aware, we should use periodic boundary conditions on all the four boundaries. Solving the navierstokes methods for the navierstokes. Pdf a rotated staggered grid finitedifference with the. Lecture 6 boundary conditions applied computational fluid. Sep 30, 2016 thereafter, it can be seen that a staggered grid and the associated physical law based finite volume method fvm are presented to clearly demonstrate one of the first remedy to probably the biggest challenge in the cfd development, a pressure. I present an accurate, stable, and simple implementation of the elastic free surface for staggered grid modeling schemes.
Using staggered grids with characteristic boundary. In this scheme, a fictitious grid is introduced to describe the discontinuous velocity at the fracture interface and a rotated staggered grid is used to accurately indicate the location of the fracture. In the example here, a noslip boundary condition is applied at the solid wall. Diagonally staggered grid for elastodynamic analysis using. It deals with the initial and boundary conditions encountered in the various cfd. Staggered grid methods in cartesian coordinates typically locate the pressure at the control volume centers i, j, k and the velocity components at the surface centers see fig. Using staggered grids with characteristic boundary conditions when. For all boundary conditions, we locate the boundary such that the boundarynormal displacements u n, corresponding outofplane shear stresses. Lectures in computational fluid dynamics of incompressible. Boundary conditions across discontinuities are not. A comparison of staggered and nonstaggered grid navier. I boundary conditions involve only one eld usually e i completely dominating in commercial codes i based on decoupling of electric and magnetic elds. A rotated staggered grid finitedifference with the absorbing boundary condition of a perfectly matched layer article pdf available in chinese science bulletin 5119. In the past many nonstaggered grid navierstokes solvers have been found to su.
Finally, we write gridpoint contributions of the new staggeredgrid stencil for fourthorder derivative approximations that are valid for the whole medium, including pml absorbing boundary conditions. Further if the c2 scheeme is used, it is not really possible to set the derivative at the boundary, if we calculate from stagnonstag not sure whether this is a bug, or the way it should be. A staggered grid is a setting for the spatial discretisation, in which the variables are not defined at the same position as in case of the collocated grid. A standard approach is to prescribe homogeneous neumann boundary conditions for p wherever noslip boundary conditions are prescribed for the velocity. The question, which boundary conditions are appropriate for. An explicit staggeredgrid method for numerical simulation. Fully staggered grid continuity, xmomentum, and ymomentum equa tions satisfied on three different sets of grid points 2 5. The new fd scheme satisfies nonwelded contact boundary conditions, unlike traditional approaches. The continuity equation is then discretized by intergating it at control volumes that coincide with each. In the staggered grid finitedifference method stress and particle velocity components are calculated at different grid points, and a faulting problem is a mixed boundary problem, therefore different implementations of fault boundary conditions have been proposed. The velocity at nodes on the wall is simply the wall velocity the pressure boundary is more complex computational fluid dynamics. Acoustic boundary condition treatments 9, 10 are then applied in a natural manner, despite the staggered character of the grid. Nonstaggered grids store all the unknowns at the same locations, as shown in figure 5 where both nonstaggered and staggered grid con. Simulating seismic wave propagation in 3d elastic media.
Solution to twodimensional incompressible navierstokes. When constructing a staggered grid, it is common to implement boundary conditions by adding an extra node across the physical boundary. Freesurface boundary conditions for elastic staggered. Simulating seismic wave propagation in 3d elastic media using staggered grid finite differences by robert w.
I want to solve two dimensional navier stokes equations on a staggered grid for the case of taylorgreen vortex. Lecture 6 boundary conditions applied computational. In principle, the proposed effective boundary saving will reduce 4 nt 4 n 2 bytes for regular grid finite difference, compared with the method of dussaud et al. The information in this chapter is divided into the following sections. Compressible subsonic flow on a staggered grid ubc computer. With the conventional staggered grids,1,2 as shown in fig.
May 08, 2010 staggered grid combines several types of nodal points located in different geometrical positions and looks rather complex and nonlogical. Staggered grid finitedifference acoustic modeling with the timedomain atmospheric acoustic propagation suite tdaaps neill p. Mixedgrid and staggeredgrid finitedifference methods. This grid topology and procedure can be seen as an hybrid colocated staggered strategy 26. Since i work on the fvm for the incompressible navierstokes equations using staggered gri. Spectral method solution of the stokes equations on. In section 4 we develop standardized compatibility rules for boundary conditions at junctions where two or more. Comparison of staggered and collocated grids for solving navierstokes equations m. Hidden beauty and elegancy of staggered grid become obvious after. The method is based on the ansatz that the elastic hookes tensor on the free surface can be taken as similar to a transversely isotropic medium. Semi staggered grid continuity and momentum equations satisfied on two different sets of grid points lo 6. Rigid boundary conditions for staggeredgrid modelling crewes. Values needed for periodic boundaries are taken from the other side of the domain when needed.
Freesurface boundary conditions for elastic staggeredgrid. Acoustic waveguide eigenmode solver based on a staggered. The nodes just outside the inlet of the system are used to assign the inlet conditions and the physical boundaries can coincide with the scalar control volume boundaries. We show that the accuracy of the fourthorder staggered grid stencil is slightly better than the mixed grid stencil when averaging of the mass acceleration term is applied to the staggered grid. Simulating seismic wave propagation in 3d elastic media using. Sbpsat finite difference discretization of acoustic wave equations.
Because of its simple geometry and boundary conditions, the lid. Write a function to create the finitedifference approximation of the 2nd derivative operator matrix for a staggered grid. Then exactly the same solution as above has not changed the number of extrema. There are many advantages with weakly enforced boundary conditions. Nonstaggered grid, overdetermined system formed by increasing the num. Improvement in the fault boundary conditions for a. The storage requirement of staggered grid based effective boundary saving is about 2 n. A range of models and numerical methods is available for solution of the initial boundary value problem ibvp of gas flow through a pipeline. Boundary conditions will be treated in more detail in this lecture. The qualitative and quantitative results confirmed that the effective thickness eth of the first soil layer be. Acoustic waveguide eigenmode solver based on a staggeredgrid. The most simple and most common boundary conditions used in computer graphics simulations are solid wall boundaries and free surface boundaries. Staggeredgrid finitedifference acoustic modeling with the.
Many researchers have proposed special treatments for outlet boundary conditions owing to lack of information at the outlet. A simple and efficient outflow boundary condition for the. In the staggered grid finitedifference method stress and particle velocity components are calculated at different grid points, and a faulting problem is a mixed boundary problem, therefore different implementations of fault boundary conditions have been. Computational fluid dynamics of incompressible flow.
The turbulent fluxes of momentum, heat, and moisture can either be imposed or computed using a bulk aerodynamic formula. Lectures in computational fluid dynamics of incompressible flow. Mixedgrid and staggeredgrid finitedifference methods for. Boundary conditions this chapter describes the boundary condition options available influent. Lines of nodes is put on the line j, and the values on the virtual line are set antisymmetric to those on the line j1, i. The question, which boundary conditions are appropriate for the poisson equation for the pressure p, is complicated. Boundary conditions in computational fluid dynamics. An explicit staggeredgrid method for numerical simulation of. Periodic boundary conditions for solving navier stokes. In section 4 the dispersion properties of the mixed grid and fourthorder staggered grid stencils are compared through a classical harmonic analysis. Acoustic boundary condition treatments, 14 are then applied in a natural manner, despite the staggered character of the grid.
Staggeredgrid finitedifference acoustic modeling with the timedomain atmospheric acoustic propagation suite tdaaps neill p. I am facing some trouble regarding implementation of this bc as i could not find enough literature about it. Rigid boundary conditions for staggeredgrid modelling. Graves abstract this article provides an overview of the application of the staggered grid finitedifference technique to model wave propagation problems in 3d elastic media. Nevertheless for high reynolds numbers, the results of staggered grid demonstrate better agreement. The most important advantage is the straightforward path to stability for highorder methods in summationbyparts form see discussion below. Note that very recent independent work by topalian et al. A staggered grid is used to store the variables where the pressure is stored at the cell. A compact and fast matlab code solving the incompressible navierstokes equations on rectangular domains mit18086 navierstokes. Especially for the yee scheme and the staggered grid, the formulation of boundary conditions and discontinuous material properties become difficult.
In this paper we propose new boundary conditions for the fourth order. Dec, 2017 for all boundary conditions, we locate the boundary such that the boundarynormal displacements u n, corresponding outofplane shear stresses. The boundary fitted coordinate method developed by thompson et al. Interpolation of velocities on staggered grid in pic.
207 916 1064 1067 291 1481 1058 170 922 810 59 264 1532 1081 983 745 1062 1207 86 1021 1080 612 215 1333 745 862 1283 32 1433 924 350 1491 688 1290 796 1072 388 56 1200 1441