https://github.com/johnrudge/compaction-code

Code to compute advection and compaction of a partially molten medium under simple shear.
https://github.com/johnrudge/compaction-code

Last synced: 23 days ago
JSON representation

Code to compute advection and compaction of a partially molten medium under simple shear.

• Host: GitHub
• URL: https://github.com/johnrudge/compaction-code
• Owner: johnrudge
• Created: 2022-03-06T12:04:09.000Z (over 2 years ago)
• Default Branch: master
• Last Pushed: 2024-02-05T09:27:01.000Z (9 months ago)
• Last Synced: 2024-02-05T10:52:42.331Z (9 months ago)
• Language: Python
• Size: 282 KB
• Stars: 0
• Watchers: 1
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.rst

Awesome Lists containing this project

README

        
README for compaction-code

Aim:

====

This code computes advection and compaction in a porous medium under

simple shear. There is an option to run with a regular rectangular mesh, 

or with a mesh that has a circular torque-free inclusion.

More details can be found in the paper:

Alisic L., Rudge J.F., Katz R.F., Wells G.N., Rhebergen S.

Compaction around a rigid, circular inclusion in partially molten rock (2014)

J. Geophys. Res. Solid Earth 119:5903-5920

http://dx.doi.org/10.1002/2013JB010906

Running:

========

run using:

> python3 main.py param.cfg

Requires: FEniCS 2019.1 (https://www.fenicsproject.org) and scipy

Contents:

=========

* main.py             : main code

* analysis.py         : library of diagnostics and benchmark quantities

* core.py             : library of core functions to deal with input and output of data, norms, etc

* mesh_gen.py         : code that interfaces with gmsh to create a rectangular mesh with a cylinder

* physics.py          : library of viscosity formulations and initial porosity fields

* param.cfg           : list of run parameters, can be copied to local run directory

* param_uniform_cylinder.cfg : list of run parameters for cylinder models with uniform initial porosity

* param_planewave.cfg : list of run parameters for plane wave benchmark, no cylinder

* param_perlin.cfg    : list of run parameters for initial noise case, no cylinder

* param_random_cylinder.cfg  : list of run parameters for cylinder models with random initial porosity, read in from file

Creating meshes and initial porosity fields:

============================================

* make_mesh.py         : create mesh, parameters provided in param_make_mesh.cfg

* make_initial_porosity.py : create coarse initial porosity field on large mesh, parameters provided in 

  param_make_initial_porosity.cfg

* interpolate_initial_porosity.py : interpolate coarse initial porosity field onto mesh used in numerical

  simulation, parameters provided in param_interpolate_initial_porosity.cfg

* param_make_mesh.cfg  : used by make_mesh.py

* param_make_initial_porosity.cfg : used by make_initial_porosity.py

* param_interpolate_initial_porosity.cfg : used by interpolate_initial_porosity.py

* run_preprocessing.py : example of pre-processing workflow; adjust file paths as necessary.