From 93e9cfe84468c9c68bae4cafd294bfea7066c420 Mon Sep 17 00:00:00 2001 From: David Seus <david.seus@ians.uni-stuttgart.de> Date: Tue, 13 Oct 2020 13:51:13 +0200 Subject: [PATCH] correct resolution explanation --- Usecases/README.md | 17 +++++++++-------- 1 file changed, 9 insertions(+), 8 deletions(-) diff --git a/Usecases/README.md b/Usecases/README.md index c74d53b..92dce42 100644 --- a/Usecases/README.md +++ b/Usecases/README.md @@ -110,29 +110,30 @@ number_of_timesteps = 20 This block sets mesh parameters, number of timesteps to be simulated, timestep size as well as whether or not a mesh study is being performed. In detail: -- `python mesh_study = True/False`: Toggle mesh study mode. +- `mesh_study = True/False`: Toggle mesh study mode. This changes the output slightly. Have a look at examples in the `mesh_study` folders. If the parameters is set to `True` the word mesh_study is appended to the output directory and all mesh sizes are output into the same directory. If set to `False`, data for each mesh size gets its own directory. -- `resolutions`: is a dictionary containing pairs `mesh_resolution: error_tol` - the key of each pair is the `mesh_resolution` parameter for `mshr` and the - value is `error_tol` is the error tolerance that is used to stop the iteration. The error criterion is the L^2-norm of subsequent iterates. Once all subsequent errors on all subdomains are lower than `error_tol`, the +- `resolutions`: is a dictionary containing pairs `mesh_resolution: error_tol`. + The key of each pair, `mesh_resolution`, is the `mesh_resolution` parameter for `mshr` determing the mesh size h and the + value `error_tol` is the error tolerance that is used to stop the iteration. The error criterion is the L^2-norm of subsequent iterates. Once all subsequent errors on all subdomains are lower than `error_tol`, the LDD solver stops and the calculation of the time step is considerd finished. Doubling the `mesh_resolution` parameter should result in halfing the grid width. This is dependent on the geometry and is not always true for the first - few values given in the list above. + few values given in the list above. This is a `mshr` issue and cannot be + remidied without changing the mesh generation to another mesh generator. If more than one entry is given in `resolutions`, all pairs are being calculated in parallel. This is usefull if one is interested in the results - of the same usecase but calcuated for different mesh sizes. - Notably, if `mesh_study` is set you should have most of the above commented out. + of the same usecase but calcuated for different mesh sizes h. + Notably, if `mesh_study` is `True`you should have most of the above commented out, and data for a mesh study is performed. Notably space-time errornorms are put out. ''' **!Warning!** Notably the last to pairs take a long time to calculate and memory shortage might be an issue depending on your machine. ''' -- `python starttimes = {0: 0.0}` starttimes gives a list of starttimes to run the simulation from. +- `starttimes = {0: 0.0}` starttimes gives a list of starttimes to run the simulation from. The list is looped over and a simulation is run with t_0 as initial time for each element t_0 in starttimes - timestep_size = 0.001 -- GitLab