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.本文提供了题为《流体注入诱发干燥多孔介质内空腔膨胀》的研究论文的配套补充数据。本补充数据集包含用于求解该研究所提瞬态模型的Mathematica代码。该代码利用了本问题的准静态与准稳态特性，可将时间作为参数而非变量参与运算。本代码的核心目标为：通过将注入速率分解为与空腔体积变化及渗流过程相关的分量，实现流体诱发空腔膨胀的预测。本附件包含题为《流体注入诱发干燥多孔介质内空腔膨胀》的研究的补充数据，可供下载。该压缩包内包含用于求解不同工况下所提模型的Mathematica代码。使用本数据集时，请遵循以下操作步骤：1. 运行名为"numerical_ODE_solve_s_GT_c"的Mathematica文件，针对本研究提及的一组无量纲参数计算时间"te"；代码将持续运行直至塑性半径"c"达到1，以此求得"te"。2. 运行"numerical_ODE_solve"文件，以获取时间区间小于"te"时的界面参数，包括空腔半径"a"、渗流前锋"s"以及渗流系数"ζ"。3. 借助"numerical_ODE_solve_s_GT_c"文件，运行代码直至塑性半径"c"超过渗流前锋"s"，以此计算时间"tp"；同时针对介于"te"与"tp"之间的时间区间，求解空腔半径"a"、渗流前锋"s"、塑性半径"c"以及渗流系数"ζ"等界面参数。4. 最后运行"numerical_ODE_solve_c_GT_s"文件，以获取时间大于"tp"时的界面参数，包括空腔半径"a"、渗流前锋"s"以及渗流系数"ζ"。各步骤所得结果将用于计算任意选定时间区间内轴对称域内的应力与位移变化，相关计算将依据前述研究论文中详述的方程完成。