skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Benchmarking the accuracy of higher-order particle methods in geodynamic models of transient flow
Abstract. Numerical models are a powerful tool for investigating the dynamic processes in the interior of the Earth and other planets, but the reliability and predictive power of these discretized models depends on the numerical method as well as an accurate representation of material properties in space and time. In the specific context of geodynamic models, particle methods have been applied extensively because of their suitability for advection-dominated processes and have been used in applications such as tracking the composition of solid rock and melt in the Earth's mantle, fluids in lithospheric- and crustal-scale models, light elements in the liquid core, and deformation properties like accumulated finite strain or mineral grain size, along with many applications outside the Earth sciences. There have been significant benchmarking efforts to measure the accuracy and convergence behavior of particle methods, but these efforts have largely been limited to instantaneous solutions, or time-dependent models without analytical solutions. As a consequence, there is little understanding about the interplay of particle advection errors and errors introduced in the solution of the underlying transient, nonlinear flow equations. To address these limitations, we present two new dynamic benchmarks for transient Stokes flow with analytical solutions that allow us to quantify the accuracy of various advection methods in nonlinear flow. We use these benchmarks to measure the accuracy of our particle algorithm as implemented in the ASPECT geodynamic modeling software against commonly employed field methods and analytical solutions. In particular, we quantify if an algorithm that is higher-order accurate in time will allow for better overall model accuracy and verify that our algorithm reaches its intended optimal convergence rate. We then document that the observed increased accuracy of higher-order algorithms matters for geodynamic applications with an example of modeling small-scale convection underneath an oceanic plate and show that the predicted place and time of onset of small-scale convection depends significantly on the chosen particle advection method. Descriptions and implementations of our benchmarks are openly available and can be used to verify other advection algorithms. The availability of accurate, scalable, and efficient particle methods as part of the widely used open-source code ASPECT will allow geodynamicists to investigate complex time-dependent geodynamic processes such as elastic deformation, anisotropic fabric development, melt generation and migration, and grain damage.  more » « less
Award ID(s):
1925677 2149126 1550901 1835673
PAR ID:
10534225
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
European Geosciences Publication
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
17
Issue:
10
ISSN:
1991-9603
Page Range / eLocation ID:
4115 to 4134
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Geophysical Journal International)
    SUMMARY Combining finite element methods for the incompressible Stokes equations with particle-in-cell methods is an important technique in computational geodynamics that has been widely applied in mantle convection, lithosphere dynamics and crustal-scale modelling. In these applications, particles are used to transport along properties of the medium such as the temperature, chemical compositions or other material properties; the particle methods are therefore used to reduce the advection equation to an ordinary differential equation for each particle, resulting in a problem that is simpler to solve than the original equation for which stabilization techniques are necessary to avoid oscillations. On the other hand, replacing field-based descriptions by quantities only defined at the locations of particles introduces numerical errors. These errors have previously been investigated, but a complete understanding from both the theoretical and practical sides was so far lacking. In addition, we are not aware of systematic guidance regarding the question of how many particles one needs to choose per mesh cell to achieve a certain accuracy. In this paper we modify two existing instantaneous benchmarks and present two new analytic benchmarks for time-dependent incompressible Stokes flow in order to compare the convergence rate and accuracy of various combinations of finite elements, particle advection and particle interpolation methods. Using these benchmarks, we find that in order to retain the optimal accuracy of the finite element formulation, one needs to use a sufficiently accurate particle interpolation algorithm. Additionally, we observe and explain that for our higher-order finite-element methods it is necessary to increase the number of particles per cell as the mesh resolution increases (i.e. as the grid cell size decreases) to avoid a reduction in convergence order. Our methods and results allow designing new particle-in-cell methods with specific convergence rates, and also provide guidance for the choice of common building blocks and parameters such as the number of particles per cell. In addition, our new time-dependent benchmark provides a simple test that can be used to compare different implementations, algorithms and for the assessment of new numerical methods for particle interpolation and advection. We provide a reference implementation of this benchmark in aspect (the ‘Advanced Solver for Problems in Earth’s ConvecTion’), an open source code for geodynamic modelling. 
    more » « less
  2. ; ; ; ; ; ; ; (, Zenodo)
    We are pleased to announce the release of ASPECT 2.5.0. ASPECT is the Advanced Solver for Problems in Earth's ConvecTion. It uses modern numerical methods such as adaptive mesh refinement, multigrid solvers, and a modular software design to provide a fast, flexible, and extensible mantle convection solver. ASPECT is available from https://aspect.geodynamics.org/ and the release is available from https://geodynamics.org/resources/aspect and https://github.com/geodynamics/aspect/releases/tag/v2.5.0 Among others this release includes the following significant changes: ASPECT now includes version 0.5.0 of the Geodynamic World Builder. (Menno Fraters and other contributors) ASPECT's manual has been converted from LaTeX to Markdown to be hosted as a website at https://aspect-documentation.readthedocs.io. (Chris Mills, Mack Gregory, Timo Heister, Wolfgang Bangerth, Rene Gassmoeller, and many others) New: ASPECT now requires deal.II 9.4 or newer. (Rene Gassmoeller, Timo Heister) ASPECT now supports a DebugRelease build type that creates a debug build and a release build of ASPECT at the same time. It can be enabled by setting the CMake option CMAKE_BUILD_TYPE to DebugRelease or by typing "make debugrelease". (Timo Heister) ASPECT now has a CMake option ASPECT_INSTALL_EXAMPLES that allows building and install all cookbooks and benchmarks. ASPECT now additionally installs the data/ directory. Both changes are helpful for installations that are used for teaching and tutorials. (Rene Gassmoeller) Changed: ASPECT now releases the memory used for storing initial conditions and the Geodynamic World Builder after model initialization unless an owning pointer to these objects is kept. This reduces the memory footprint for models initialized from large data files. (Wolfgang Bangerth) Added: Various helper functions to distinguish phase transitions for different compositions and compositional fields of different types. (Bob Myhill) Added: The 'adiabatic' initial temperature plugin can now use a spatially variable top boundary layer thickness read from a data file or specified as a function in the input file. Additionally, the boundary layer temperature can now also be computed following the plate cooling model instead of the half-space cooling model. (Daniel Douglas, John Naliboff, Juliane Dannberg, Rene Gassmoeller) New: ASPECT now supports tangential velocity boundary conditions with GMG for more geometries, such as 2D and 3D chunks. (Timo Heister, Haoyuan Li, Jiaqi Zhang) New: Phase transitions can now be deactivated outside a given temperature range specified by upper and lower temperature limits for each phase transition. This allows implementing complex phase diagrams with transitions that intersect in pressure-temperature space. (Haoyuan Li) New: There is now a postprocessor that outputs the total volume of the computational domain. This can be helpful for models using mesh deformation. (Anne Glerum) New: Added a particle property 'grain size' that tracks grain size evolution on particles using the 'grain size' material model. (Juliane Dannberg, Rene Gassmoeller) Fixed: Many bugs, see link below for a complete list. (Many authors. Thank you!). A complete list of all changes and their authors can be found at https://aspect.geodynamics.org/doc/doxygen/changes_between_2_84_80_and_2_85_80.html Wolfgang Bangerth, Juliane Dannberg, Menno Fraters, Rene Gassmoeller, Anne Glerum, Timo Heister, Bob Myhill, John Naliboff, and many other contributors. 
    more » « less
  3. ; ; ; ; (, Geoscientific Model Development)
    Abstract. Due to the increasing availability of high-performance computing over the past few decades, numerical models have become an important tool for research in geodynamics.Several generations of mantle convection software have been developed, but due to their differing methods and increasing complexity it is important to evaluate the accuracy of each new model generation to ensure published geodynamic research is reliable and reproducible.Here we explore the accuracy of the open-source, finite-element codes ASPECT and CitcomS as a function of mesh spacing using low to moderate-Rayleigh-number models in steady-state thermal convection.ASPECT (Advanced Solver for Problems in Earth's ConvecTion) is a new-generation mantle convection code that enables modeling global mantle convection with realistic parameters and complicated physical processes using adaptive mesh refinement (Kronbichler et al., 2012; Heister et al., 2017).We compare the ASPECT results with calculations from the finite-element code CitcomS (Zhong et al., 2000; Tan et al., 2006; Zhong et al., 2008), which has a long history of use in the geodynamics community.We find that the globally averaged quantities, i.e., root-mean-square (rms) velocity, mean temperature, and Nusselt number at the top and bottom of the shell, agree to within 1 % (and often much better) for calculations with sufficient mesh resolution.We also show that there is excellent agreement of the time evolution of both the rms velocity and the Nusselt numbers between the two codes for otherwise identical parameters.Based on our results, we are optimistic that similar agreement would be achieved for calculations performed at the convective vigor expected for Earth, Venus, and Mars. 
    more » « less
  4. ; ; ; ; ; (, Geoscientific Model Development)
    Abstract. Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work.For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores. 
    more » « less
  5. ; ; ; (, Geophysical Research Letters)
    Abstract The rheology of crustal mushes is a crucial parameter controlling melt segregation and magma flow. However, the relations between mush dynamics and crystal size and shape distribution remain poorly understood because of the complexity of melt‐crystal and crystal‐crystal interactions. We performed analog experiments to characterize the mechanisms that control pore space reduction associated with repacking. Three suspensions of monodisperse particles with different geometries and aspect ratios (1:1, 2:1, 4:1) in a viscous fluid were tested. Our results show that particle aspect ratios strongly control the melt extraction processes. We identify two competing mechanisms that enable melt extraction at grain scale. The first mechanism leads to continuous deformation and melt extraction and is associated with “diffuse” frictional dissipation between neighboring particles. The second is stochastic, localized, and nearly instantaneous and is associated with the development and destruction of force chains percolating through the granular assembly. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email