Adaptive mesh refinement

Adaptive mesh refinement
This article is about the use of adaptive meshing in numerical analysis. See Subdivision surface for the use of adaptive techniques in Computer Graphics modelling.

In numerical analysis, adaptive mesh refinement is a method of adaptive meshing. Central to any Eulerian method is the manner in which it discretizes the continuous domain of interest into a grid of many individual elements. This grid may be static, established once and for all at the beginning of the computation, or it may be dynamic, tracking the features of the result as the computation progresses. If the computation has features which one wants to track which are much smaller than the overall scale of the problem, and which move in time, then one must either include many more static grids to cover the region of interest, or adopt a dynamic scheme.[1]

The advantages of a dynamic gridding scheme are:

  1. Increased computational savings over a static grid approach.
  2. Increased storage savings over a static grid approach.
  3. Complete control of grid resolution, compared to the fixed resolution of a static grid approach, or the Lagrangian-based adaptivity of smoothed particle hydrodynamics.

Introduction to adaptive mesh refinement

The image above shows the grid structure of an AMR calculation of a shock impacting an inclined slope. Each of the boxes is a grid; the more boxes it is nested within, the higher the level of refinements. As the image shows, the algorithm uses high resolution grids only at the physical locations and times where they are required.

In a series of papers, Marsha Berger, Joseph Oliger, and Phillip Colella developed an algorithm for dynamic gridding called local adaptive mesh refinement. The algorithm begins with the entire computational domain covered with a coarsely resolved base-level regular Cartesian grid. As the calculation progresses, individual grid cells are tagged for refinement, using a criterion that can either be user-supplied (for example mass per cell remains constant, hence higher density regions are more highly resolved) or based on Richardson extrapolation.

All tagged cells are then refined, meaning that a finer grid is overlaid on the coarse one. After refinement, individual grid patches on a single fixed level of refinement are passed off to an integrator which advances those cells in time. Finally, a correction procedure is implemented to correct the transfer along coarse-fine grid interfaces, to ensure that the amount of any conserved quantity leaving one cell exactly balances the amount entering the bordering cell. If at some point the level of refinement in a cell is greater than required, the high resolution grid may be removed and replaced with a coarser grid.

This allows the user to solve problems that are completely intractable on a uniform grid; for example, astrophysicists have used AMR to model a collapsing giant molecular cloud core down to an effective resolution of 131,072 cells per initial cloud radius, corresponding to a resolution of 1015 cells on a uniform grid.[2]

References

Berger, M. J.; Colella, P. (1989). "Local adaptive mesh refinement for shock hydrodynamics". J. Comput. Phys. (Elsevier) 82: 64–84.

See also


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать курсовую

Look at other dictionaries:

  • Adaptive Meshing — Als Adaptives Meshing bezeichnet man ein spezielles Verfahren zur Diskretisierung eines geometrischen Modells in einer numerischen Simulationsrechnung wie der Finite Elemente Methode oder der Finite Integral Methode. In der herkömmlichen… …   Deutsch Wikipedia

  • Numerical relativity — is one of the branches of general relativity that uses numerical methods and algorithms to solve and analyze problems. To this end, supercomputers are often employed to study black holes, gravitational waves, neutron stars and many other… …   Wikipedia

  • Phillip Colella — Nationality American Fields …   Wikipedia

  • Cactus Framework — Infobox Software name = Cactus Framework developer = [http://www.cactuscode.org/ Cactus Team] latest release version = 4.0 beta 15 latest release date = November 2, 2004 operating system = Cross platform genre = Software framework license = LGPL… …   Wikipedia

  • Сетка (расчетная) — Функция одной переменной Ф, заданная на структурированной сетке {xk} Расчетная (вычислительная) сетка совокупность точек (сеточных уз …   Википедия

  • Практическое применение раскраски графов — Эту статью следует викифицировать. Пожалуйста, оформите её согласно правилам оформления статей. Раскраска графов практически применяется (постановку задачи различиных раскрасок здесь обсуждаться не будет) дл …   Википедия

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • Гидродинамика сглаженных частиц — (англ. Smoothed Particle Hydrodynamics, SPH)  вычислительный метод для симуляции жидкостей и газов. Используется во многих областях исследований, включая астрофизику, баллистику, вулканологию и океанографию. Метод гидродинамики… …   Википедия

  • Richard Klein (astronomer) — Richard Klein is an Adjunct Professor of Astronomy at the University of California, Berkeley and a Scientific Staff Member at the Lawrence Livermore National Laboratory (LLNL). Klein received his bachelor s degree in physics from Rensselaer… …   Wikipedia

  • Smoothed particle hydrodynamics — (SPH) is a computational method used for simulating fluid flows. It has been used in many fields of research, including astrophysics, ballistics, vulcanology and oceanology. It is a mesh free Lagrangian method (where the co ordinates move with… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”