Aperture and flow patterns in synthetic fractures with varying JRC for "On the anisotropy of the aperture of rough fractures"

This set contains aperture and fluid flow simulation data (pressure, pressure gradient, and velocity distributions) used for upscaling equivalent aperture tensors for eight samples of synthetic fractures with varying joint roughness coefficients (JRC) and varying contact areas.
<br>
The aperture data was generated through fractional Brownian motion by combining multiple steps of simplex noise. For that purpose, we used the Python library <strong>noise </strong>(https://github.com/caseman/noise). The fluid flow simulations were carried out with the academic, multi-institutional numerical framework <strong>CSMP++</strong> (https://www.geoenergynumerics.com/software/open-csmp), developed in C++.
<br>
The mesh files in this dataset are in VTK format and may be read using the Python library <strong>meshio </strong>(https://github.com/nschloe/meshio). The files are named <strong>sample-</strong><em><strong>{sample_name}</strong></em><strong>_JRC-</strong><em><strong>{JRC_value}</strong></em><strong>_CA-</strong><em><strong>{contact_area_value}</strong></em><strong>_</strong><em><strong>{property}</strong></em><strong>.vtk</strong>. For each sample, namely F1-F8, and each contact area value (between 0.0 and 0.4), there are 7 VTK files. These are as follows:

<strong>sample-</strong><em><strong>{sample_name}</strong></em><strong>_JRC-</strong><em><strong>{JRC_value}</strong></em><strong>_CA-</strong><em><strong>{contact_area_value}_</strong></em><strong>aperture.vtk</strong>: these files contain the aperture distribution, that is, the opening of the fractures, for each sample and each contact area value. The aperture values are stored in the center of the mesh cells and are constant for each cell.
<strong>sample-</strong><em><strong>{sample_name}</strong></em><strong>_JRC-</strong><em><strong>{JRC_value}</strong></em><strong>_CA-</strong><em><strong>{contact_area_value}_</strong></em><strong>fluid-pressure-{</strong><em><strong>#n</strong></em><strong>}.vtk</strong>: these files contain the pressure distribution for problem number <em>#n</em> (1 or 2), for each sample and each contact area value. The pressure values are stored in the nodes of the mesh.
<strong>sample-</strong><em><strong>{sample_name}</strong></em><strong>_JRC-</strong><em><strong>{JRC_value}</strong></em><strong>_CA-</strong><em><strong>{contact_area_value}_</strong></em><strong>fluid-pressure-gradient-{</strong><em><strong>#n</strong></em><strong>}.vtk</strong>: these files contain the pressure gradient distribution for problem number <em>#n</em> (1 or 2), for each sample and each contact area value. The pressure gradient values are derived from the nodal pressure values. Like aperture, they are stored in the center of the mesh cells and are constant for each cell.
<strong>sample-</strong><em><strong>{sample_name}</strong></em><strong>_JRC-</strong><em><strong>{JRC_value}</strong></em><strong>_CA-</strong><em><strong>{contact_area_value}_</strong></em><strong>velocity-{</strong><em><strong>#n</strong></em><strong>}.vtk</strong>: these files contain the fluid velocity for problem number <em>#n</em> (1 or 2), for each sample and each contact area value. The velocity values are derived from the pressure gradient and the aperture values using the cubic law: velocity = (aperture^3/12*viscosity)*pressure_gradient. Like the aperture and pressure gradient, the velocity values are stored in the center of the mesh cells and are constant for each cell.

<br>
If you use this dataset in your work, include this citation:
Ferreira, Carlos A. S.; Nick, Hamidreza M. (2023). Aperture and flow patterns in synthetic fractures with varying JRC for "On the anisotropy of the aperture of rough fractures". Technical University of Denmark. Dataset. https://doi.org/10.11583/DTU.23586945.

本数据集包含用于开展8组不同节理粗糙度系数（Joint Roughness Coefficient, JRC）及接触面积的人工裂隙样本的等效孔径张量尺度升尺度研究所需的孔径与流体渗流模拟数据，涵盖压力、压力梯度及速度分布。

孔径数据通过分形布朗运动结合多步单纯形噪声生成，本次研究使用Python库`noise`（https://github.com/caseman/noise）完成该部分生成工作。流体渗流模拟则依托多机构联合开发的学术级C++数值框架**CSMP++**（https://www.geoenergynumerics.com/software/open-csmp）开展。

本数据集内的网格文件均采用VTK格式，可通过Python库`meshio`（https://github.com/nschloe/meshio）读取。文件命名规则为：`sample-{sample_name}_JRC-{JRC_value}_CA-{contact_area_value}_{property}.vtk`。针对F1至F8共8组样本，以及0.0至0.4区间内的所有接触面积取值，每组均包含7个VTK文件，具体如下：
1. `sample-{sample_name}_JRC-{JRC_value}_CA-{contact_area_value}_aperture.vtk`：该类文件存储对应样本与接触面积取值下的孔径分布（即裂隙开度）数据。孔径数值存储于网格单元中心，且每个单元内的孔径值恒定不变。
2. `sample-{sample_name}_JRC-{JRC_value}_CA-{contact_area_value}_fluid-pressure-{#n}.vtk`：该类文件存储对应样本、接触面积取值下第#n（1或2）个渗流算例的压力分布数据。压力数值存储于网格节点处。
3. `sample-{sample_name}_JRC-{JRC_value}_CA-{contact_area_value}_fluid-pressure-gradient-{#n}.vtk`：该类文件存储对应样本、接触面积取值下第#n（1或2）个渗流算例的压力梯度分布数据。压力梯度数值由节点压力值推导得到，与孔径数据一致，存储于网格单元中心且每个单元内数值恒定。
4. `sample-{sample_name}_JRC-{JRC_value}_CA-{contact_area_value}_velocity-{#n}.vtk`：该类文件存储对应样本、接触面积取值下第#n（1或2）个渗流算例的流体速度数据。速度数值由压力梯度与孔径数值通过立方定律推导得到，公式为：速度 = (孔径³/(12×黏度))×压力梯度。与孔径及压力梯度数据一致，速度数值存储于网格单元中心且每个单元内数值恒定。

若您在研究中使用本数据集，请引用如下文献：
Ferreira, Carlos A. S.; Nick, Hamidreza M. (2023). 变JRC人工裂隙的孔径与渗流特征——基于"粗糙裂隙孔径各向异性"研究. 丹麦技术大学. 数据集. https://doi.org/10.11583/DTU.23586945.