This block sets mesh parameters, number of timesteps to be simulated, timestep
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.
size as well as whether or not a mesh study is being performed.
In detail:
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`
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
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
to the output directory and all mesh sizes are output into the same
directory.
directory.
If set to `False`, data for each mesh size gets its own directory.
If set to `False`, data for each mesh size gets its own directory.
-`resolutions`: is a dictionary containing pairs `mesh_resolution: error_tol`
-`resolutions`: is a dictionary containing pairs `mesh_resolution: error_tol`.
the key of each pair is the `mesh_resolution` parameter for `mshr` and the
The key of each pair, `mesh_resolution`, is the `mesh_resolution` parameter for `mshr`determing the mesh size h 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
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.
LDD solver stops and the calculation of the time step is considerd finished.
Doubling the `mesh_resolution` parameter should result in halfing the grid
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
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
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.
of the same usecase but calcuated for different mesh sizes h.
Notably, if `mesh_study` is set you should have most of the above commented out.
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!**
**!Warning!**
Notably the last to pairs take a long time to calculate and memory shortage
Notably the last to pairs take a long time to calculate and memory shortage
might be an issue depending on your machine.
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
The list is looped over and a simulation is run with t_0 as initial time
