chaovite/crack_pipe: crack_pipe_1.0

This repository contains the Matlab code that couples a multi-section pipes with multiple rectangular fluid filled cracks. The subroutines for SBP operators and imex time integrators are contributed by Ossian O'Reilly. Currently, the 1d pipe model can be used to couple with multiple square cracks (assuming Poiseuille's flow across the width direction). However, the 2D radial symmetric pipe is only allow to be coupled to a crack at the bottom while the crack itself can be 3D tabular crack. This model captures acoustic gravity waves in a 2D radially symmetric conduit and crack waves along 3D crack with proper treatment of viscous boundary layers in both the conduit and the crack. This code was used in publication: Liang, C., L. Karlstrom, and E. M. Dunham, 2019. Magma oscillations in a conduit-reservoir system, applications to very long period (VLP) seismicity at basaltic volcanoes-Part I: Theory, Journal of Geophysical Research, in press