Fluid injection - induced cavity expansion in dry porous medium

Supplementary data accompanying the research paper titled "Fluid Injection-Induced Cavity Expansion in Dry Porous Medium" is provided herein. This supplementary dataset includes the Mathematica code employed for solving the proposed transient model. This code capitalizes on the quasi-static and quasi-stationary characteristics of the problem, enabling the treatment of time as a parameter rather than a variable. The primary objective of this code is to facilitate predictions of fluid-induced cavity expansion through the partitioning of the injection rate into components related to changes in cavity volume and infiltration.

The file attached contains the supplementary data for the research study titled "Fluid Injection-Induced Cavity Expansion in Dry Porous Medium" and is available for download. Inside this zip file, you will find the Mathematica code used to solve the proposed model for different cases. To make use of this dataset, follow these steps:
1. Execute the Mathematica file named "numerical_ODE_solve_s_GT_c" to calculate the time "te" for a set of dimensionless parameters referenced in the study. The code computes "te" by running until the plastic radius "c" reaches a value of 1.
2. Run the file "numerical_ODE_solve" to obtain interface parameters such as the cavity radius "a," permeation front "s," and permeation coefficient "ζ" for time intervals less than "te".
3. Utilize the file "numerical_ODE_solve_s_GT_c" to calculate "tp" by running the code until the plastic radius "c" overtakes the permeation front "s." Execute the same file to derive interface parameters like the cavity radius "a," permeation front "s," plastic radius "c," and permeation coefficient "ζ" for time intervals greater than "te" but less than "tp".
4. Finally, run the file "numerical_ODE_solve_c_GT_s" to obtain interface parameters such as the cavity radius "a," permeation front "s," and permeation coefficient "ζ" for times greater than "tp".

The results derived from each of these procedural steps will be utilized to compute variations in stress and displacement within the axisymmetric domain at any selected time interval. These computations will be performed using the equations detailed in the aforementioned research paper.