MATLAB code used in calculations of 1-D diffusive growth of mushy layers with perfectly and imperfectly conducting boundaries, and analysis of linear stability to convective perturbation

This dataset includes data and processing code from calculations of diffusive mushy layer growth and linear convective stability analysis of this system, and is made available under the GNU General Public License v3:

http://www.gnu.org/licenses/gpl.html

Note that thiscode has only been tested on original author's system, and we accept no responsibility for adverse consequences if used on other systems.

Use at your own risk.

In particular, be wary that MATLAB scripts supplied may include code to modify the system path/write data/clear data from a MATLAB workspace.  Use with caution.

Simulation results produced using code written by Joseph Hitchen and edited by Andrew Wells for use with MATLAB. This code base was used to generate figures contained in the PhD thesis:

Hitchen (2017), “Understanding the growth and convective instability of mushy layers, with application to young sea ice”, PhD Thesis University of Oxford.

and two accompanying papers (in press at the Journal of Fluid Mechanics at the time of initial submission of this repository):

Joseph R. Hitchen & Andrew J. Wells, in press, J. Fluid Mech. doi:10.1017/jfm.2024.808 

Structure of mushy layers grown from perfectly and imperfectly conducting boundaries. Part 1. Diffusive solidification 

Joseph R. Hitchen & Andrew J. Wells, in press, J. Fluid Mech. doi:10.1017/jfm.2024.809

Structure of mushy layers  grown from perfectly and imperfectly conducting boundaries. Part 2. Onset of convection

Contact andrew.wells@physics.ox.ac.uk for further details.