It has been a while since the last development update, but in my opinion we had some significant improvements in the last months that might be missed in the many items of the newsletter, so I would like to emphasize them a bit:
As discussed in the online user meetings last year, we now show all available ASPECT intput parameters online at https://aspect.geodynamics.org/doc/parameter_view/parameters.xml (there is a link from aspect.geodynamics.org to this site). The web-view is more interactive and easier to browse than the Chapter in the manual, and should make it simpler to see which options are available. You can click on subsections to expand them, and on parameters to show their documentation. The parameters will be updated periodically in sync with the manual we link to from the website. If you have suggestions for improvements (or know something about web development and want to improve design or functionality) let us know or write a pull request to the
wwwbranch of the ASPECT repository.
We recently merged some pull requests that increase the output resolution of the dynamic topography, geoid, and heat flux visualization output (the one in the vtk, hdf5, or gnuplot files). If you are interested in the visualization output of any of these quantities you should consider updating your ASPECT version. The heat flux option is disabled by default, but can be safely enabled as long as you have no prescribed heat flux boundary and no material flow through the boundary you are interested in. The old output was not wrong, but it was cell-wise constant, which is a poorer approximation to the real values.
Ludovic Jeanniot contributed a nonlinear solver scheme called
no Advection, no Stokesthat does not solve the equations, but just outputs the initial condition. This is very helpful if you are stuck with a model setup where the solver does not converge, and you do not why. Having an output of the initial condition can show you if there is a problem with one of your settings (e.g. if density becomes 0, or your viscosity variation is bigger than expected, or you prescribe inflow on all boundaries).
Jon Robey implemented an alternative method to advect compositional fields. It is called Volume-of-fluid method and is relatively cheap and extremely effective at advecting sharp boundaries between compositions without numerical diffusion. There is currently little documentation, and it is limited to 2D box models, but there are a few test cases in the tests folder that might act as starting point if your problem is suitable for this method.
There is an open pull request that I would like to bring up for a discussion with the broader community, because I would like feedback on it. For a while we had the option to increase the visualization resolution compared to the grid ASPECT computes on. The reasoning for that is that internally ASPECT uses quadratic polynomials to approximate most variables (except pressure), and Paraview only visualizes them as linear functions. This means the actual solution is more accurate than what you see in Paraview. Using this method evaluates the variable at the mid-points of each cell and face and splits all cells into smaller output cells. The result is a more accurate output that resolves smaller features for essentially no additional cost. However, it might be confusing to have a different number of cells in the output than in the computation. I would welcome opinions on this, and if you did not know about the setting, I encourage you to test it, it is in the
subsection Visualizationand the parameter is called
Interpolate output. The pull request is here:
That is it for today, let us know below if you have comments or questions about the features, and thanks for reading.