Adaptive mesh refinement: Difference between revisions
m citations corrected and completed |
|||
| Line 8: | Line 8: | ||
[[File:amrgridimg.jpg|framed|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.]] |
[[File:amrgridimg.jpg|framed|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 [[scientific paper|papers]], [[Marsha Berger]], Joseph Oliger, and [[Phillip Colella]] developed an [[algorithm]] for dynamic gridding called ''local adaptive mesh refinement'' |
In a series of [[scientific paper|papers]], [[Marsha Berger]], Joseph Oliger, and [[Phillip Colella]] developed an [[algorithm]] for dynamic gridding called ''local adaptive mesh refinement'' <ref>{{cite journal |last1=Berger |first1=Marsha J. |last2=Oliger |first2=Joseph |date=1984 |title=Adaptive mesh refinement for hyperbolic partial differential equations |url=https://apps.dtic.mil/sti/pdfs/ADA130162.pdf |journal=Journal of Computational Physics |volume=53 |issue=3 |pages=484--512 |doi=10.1016/0021-9991(84)90073-1,access-date=2021-07-22}} |
||
</ref><ref>{{cite journal |last1=Berger |first1=Marsha J. |last2=Colella |first2=Philipp |date=1989 |title=Local adaptive mesh refinement for shock hydrodynamics |url=https://crd.lbl.gov/assets/pubs_presos/AMCS/ANAG/A113.pdf |journal= Journal of Computational Physics|volume=82 |issue=1 |pages=64--84 |doi=10.1016/0021-9991(89)90035-1}}</ref>. The algorithm begins with the entire computational [[Domain of a function|domain]] covered with a coarsely resolved base-level regular [[Cartesian coordinate system|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. |
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. |
||
| Line 14: | Line 15: | ||
This allows the user to solve problems that are completely intractable on a [[regular grid|uniform grid]]; for example, [[astrophysics|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 10<sup>15</sup> cells on a uniform grid.<ref>{{cite journal |last1=Klein |first1=Richard |year=1999 |title=Star formation with 3-D adaptive mesh refinement: the collapse and fragmentation of molecular clouds |journal=Journal of Computational and Applied Mathematics |volume=109 |issue=1–2 |pages=123–152 |doi=10.1016/S0377-0427(99)00156-9 |doi-access=free }}</ref> |
This allows the user to solve problems that are completely intractable on a [[regular grid|uniform grid]]; for example, [[astrophysics|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 10<sup>15</sup> cells on a uniform grid.<ref>{{cite journal |last1=Klein |first1=Richard |year=1999 |title=Star formation with 3-D adaptive mesh refinement: the collapse and fragmentation of molecular clouds |journal=Journal of Computational and Applied Mathematics |volume=109 |issue=1–2 |pages=123–152 |doi=10.1016/S0377-0427(99)00156-9 |doi-access=free }}</ref> |
||
Advanced mesh refinement has been introduced via functionals.<ref>{{cite book |last1=Huang| first1=Weizhang| last2=Russell|first2=Robert D.| title=Adaptive Moving Mesh |
Advanced mesh refinement has been introduced via functionals.<ref>{{cite book |last1=Huang| first1=Weizhang| last2=Russell|first2=Robert D.| title=Adaptive Moving Mesh Method |date=2010 |publisher=Springer |isbn=978-1-4419-7916-2}}</ref> Functionals allow the ability to generate grids and provide mesh adaptation. Some advanced functionals include the Winslow and modified Liao functionals.<ref>{{cite journal |last1=Khattri | first1=Sanjay Kumar |title=Grid generation and adaptation by functionals |date=2007 |url=https://www.scielo.br/j/cam/a/KQfnf3bLSnzRcdyJCmCkTzL/?format=pdf&lang=en |journal=Computational & Applied Mathematics |volume=26 |issue=2 |pages=235--249 |access-date=2021-07-22}}</ref> |
||
==Applications of adaptive mesh refinement== |
==Applications of adaptive mesh refinement== |
||
| Line 27: | Line 28: | ||
#Complete control of grid resolution, compared to the fixed resolution of a static grid approach, or the Lagrangian-based adaptivity of [[smoothed particle hydrodynamics]]. |
#Complete control of grid resolution, compared to the fixed resolution of a static grid approach, or the Lagrangian-based adaptivity of [[smoothed particle hydrodynamics]]. |
||
#Compared to pre-tuned static meshes, the adaptive approach requires less detailed a priori knowledge on the evolution of the solution. |
#Compared to pre-tuned static meshes, the adaptive approach requires less detailed a priori knowledge on the evolution of the solution. |
||
#The computational costs inherit properties of the physical system.<ref> |
#The computational costs inherit properties of the physical system.<ref>{{cite journal |last1=Popinet |first1=Stéphane |date=2015 |title=A quadtree-adaptive multigrid solver for the Serre–Green–Naghdi equations |url=https://www.archives-ouvertes.fr/hal-01163101/document |journal=Journal of Computational Physics |volume=302 |issue= |pages=336-358 |doi=10.1016/j.jcp.2015.09.009 |access-date=2021-07-22}}</ref> |
||
==References== |
==References== |
||
* Berger, M. J.; Colella, P. (1989). "Local adaptive mesh refinement for shock hydrodynamics". J. Comput. Phys. (Elsevier) 82: 64–84. |
|||
{{reflist}} |
{{reflist}} |
||
Revision as of 10:46, 22 July 2021
In numerical analysis, adaptive mesh refinement (AMR) is a method of adapting the accuracy of a solution within certain sensitive or turbulent regions of simulation, dynamically and during the time the solution is being calculated. When solutions are calculated numerically, they are often limited to pre-determined quantified grids as in the Cartesian plane which constitute the computational grid, or 'mesh'. Many problems in numerical analysis, however, do not require a uniform precision in the numerical grids used for graph plotting or computational simulation, and would be better suited if specific areas of graphs which needed precision could be refined in quantification only in the regions requiring the added precision. Adaptive mesh refinement provides such a dynamic programming environment for adapting the precision of the numerical computation based on the requirements of a computation problem in specific areas of multi-dimensional graphs which need precision while leaving the other regions of the multi-dimensional graphs at lower levels of precision and resolution.
This dynamic technique of adapting computation precision to specific requirements has been accredited to Marsha Berger, Joseph Oliger, and Phillip Colella who developed an algorithm for dynamic gridding called local adaptive mesh refinement. The use of AMR has since then proved of broad use and has been used in studying turbulence problems in hydrodynamics as well as in the study of large scale structures in astrophysics as in the Bolshoi Cosmological Simulation.
Development of adaptive mesh refinement
In a series of papers, Marsha Berger, Joseph Oliger, and Phillip Colella developed an algorithm for dynamic gridding called local adaptive mesh refinement [1][2]. 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.[3]
Advanced mesh refinement has been introduced via functionals.[4] Functionals allow the ability to generate grids and provide mesh adaptation. Some advanced functionals include the Winslow and modified Liao functionals.[5]
Applications of adaptive mesh refinement
When calculating a solution to the shallow water equations, the solution (water height) might only be calculated for points every few feet apart—and one would assume that in between those points the height varies smoothly. The limiting factor to the resolution of the solution is thus the grid spacing: there will be no features of the numerical solution on scales smaller than the grid-spacing. Adaptive mesh refinement (AMR) changes the spacing of grid points, to change how accurately the solution is known in that region. In the shallow water example, the grid might in general be spaced every few feet—but it could be adaptively refined to have grid points every few inches in places where there are large waves.
If the region in which higher resolution is desired remains localized over the course of the computation, then static mesh refinement can be used - in which the grid is more finely spaced in some regions than others, but maintains its shape over time.
The advantages of a dynamic gridding scheme are:
- Increased computational savings over a static grid approach.
- Increased storage savings over a static grid approach.
- Complete control of grid resolution, compared to the fixed resolution of a static grid approach, or the Lagrangian-based adaptivity of smoothed particle hydrodynamics.
- Compared to pre-tuned static meshes, the adaptive approach requires less detailed a priori knowledge on the evolution of the solution.
- The computational costs inherit properties of the physical system.[6]
References
- ^ Berger, Marsha J.; Oliger, Joseph (1984). "Adaptive mesh refinement for hyperbolic partial differential equations" (PDF). Journal of Computational Physics. 53 (3): 484--512. doi:10.1016/0021-9991(84)90073-1,access-date=2021-07-22.
- ^ Berger, Marsha J.; Colella, Philipp (1989). "Local adaptive mesh refinement for shock hydrodynamics" (PDF). Journal of Computational Physics. 82 (1): 64--84. doi:10.1016/0021-9991(89)90035-1.
- ^ Klein, Richard (1999). "Star formation with 3-D adaptive mesh refinement: the collapse and fragmentation of molecular clouds". Journal of Computational and Applied Mathematics. 109 (1–2): 123–152. doi:10.1016/S0377-0427(99)00156-9.
- ^ Huang, Weizhang; Russell, Robert D. (2010). Adaptive Moving Mesh Method. Springer. ISBN 978-1-4419-7916-2.
- ^ Khattri, Sanjay Kumar (2007). "Grid generation and adaptation by functionals". Computational & Applied Mathematics. 26 (2): 235--249. Retrieved 2021-07-22.
- ^ Popinet, Stéphane (2015). "A quadtree-adaptive multigrid solver for the Serre–Green–Naghdi equations". Journal of Computational Physics. 302: 336–358. doi:10.1016/j.jcp.2015.09.009. Retrieved 2021-07-22.