Aneurysmal haemodynamics: A three-dimensional fluid-structure interaction approach
Data files
Dec 02, 2025 version files 289.37 MB
-
FSI_dt_0.0125_puls_0.00075.zip
101.38 MB
-
FSI_dt_0.0125_puls_0.00500.zip
102.05 MB
-
Haemodynamics_dt_0.0125.zip
85.76 MB
-
README.md
5.85 KB
-
TAWSS_OSI.zip
171.52 KB
Abstract
In numerical studies of blood flow in aneurysms, it is essential to consider movements of the arterial wall and their interactions with the fluid. Here, mechanics equations are augmented to an Euler--Lagrange formulation, facilitating the study of blood flow in a pathological geometry during the cardiac cycle. The arterial morphology and pulse are represented by the transformation of the system of equations into a body-fitted approach via generalized curvilinear coordinates. Dynamic three-dimensional governing partial differential equations are discretized with the finite volume method on a collocated grid. The results obtained with this mathematical model reveal that the pulsating wall influences the velocity field, with prominent recirculation zones. In addition, localised lateral pressure gradients are observed within the pathological region. Increased pulsatility causes large augmentations of the Reynolds and Womersley numbers, with divergence from the hydrodynamic case. Analysis of biomedical factors such as the time-averaged wall shear stress and oscillatory shear index on the oscillating boundary provides valuable insight into the shear forces on the arterial wall compared to rigid wall dynamics. These findings underscore the need to model fluid--structure interactions in the context of aneurysm progression and cardiovascular risk assessment and to consider these interactions in patient-based biomedical approaches.
The data is a supplementary part of the paper "Aneurysmal Haemodynamics: A Three-Dimensional Fluid--Structure Interaction Approach", Royal Society Open Science. The data contain the velocity and pressure profiles of the simulated cases in the paper - heamodynamic and Fluid--Structure interaction cases. Additionally, contain the TAWSS and OSI indices data from the same simulations. We include a script for visualizations of these data in each folder. For more details, please read the README file.
Dryad DOI: https://doi.org/10.5061/dryad.pc866t22m
File Formats: ZIP archives, MATLAB (.m) scripts, ASCII (.dat) data files, GIF visualizations
Software Required: MATLAB R2023a or newer
License: CC0 (Public Domain)
OVERVIEW OF THE DATASET
This dataset contains numerical simulation results from a three-dimensional computational study of blood flow in an aneurysmal vessel geometry. The simulations include:
- One rigid-wall haemodynamic case
- Two fluid–structure interaction (FSI) cases with:
- Low pulsatility
- High pulsatility
The dataset supports visualization and post-processing of:
- Velocity fields
- Pressure fields
- Time-Averaged Wall Shear Stress (TAWSS)
- Oscillating Shear Index (OSI)
The data may be reused for solver validation, biomedical flow analysis, rigid vs deformable wall comparison, and educational purposes in computational fluid dynamics (CFD).
All data are fully synthetic and anonymized. No human subject data are included.
MAIN DATASET FOLDERS
The dataset is distributed across four main ZIP archives:
- README.txt
- Haemodynamics_dt_0.0125.zip
- FSI_dt_0.0125_puls_0.00075.zip
- FSI_dt_0.0125_puls_0.00500.zip
- TAWSS_OSI.zip
NAMING CONVENTIONS
- Haemodynamics = Rigid-wall flow
- FSI = Fluid–structure interaction
- dt_0.0125 = Time step size of the simulation (dt = 0.0125)
- puls_0.00075 = Low pulsation of the deformable wall
- puls_0.00500 = High pulsation of the deformable wall
GRID RESOLUTIONS
Each simulation case is provided on two structured grids:
- 40 x 20 x 20
- 50 x 25 x 25
These correspond to the number of computational cells in the axial, radial, and circumferential directions.
CONTENTS OF FSI ZIP FILES
Both FSI_dt_0.0125_puls_0.00075.zip and FSI_dt_0.0125_puls_0.00500.zip contain two folders:
- 40x20x20_dt=0.0125_puls=...
- 50x25x25_dt=0.0125_puls=...
Each grid-resolution folder contains:
GIF FILES:
- 2D_pressure.gif
- 2D_velocity_streamlines.gif
- 3D_pressure.gif
- 3D_u_velocity.gif
MATLAB SCRIPTS:
- contours_2D.m
- contours_3D.m
- main_visuals.m
Additionally, the high pulsatility case also includes:
- mat_to_dat.m
DATA FOLDER:
A folder named "Data" containing:
- p-1_.mat to p-55_.mat (pressure fields)
- u-1_*.mat (x-velocity)
- v-1_*.mat (y-velocity)
- w-1_*.mat (z-velocity)
- x-2.mat to x-55.mat and x-steady.mat
- y-2.mat to y-55.mat and y-steady.mat
- z-2.mat to z-55.mat and z-steady.mat
For the low pulsatility case, pressure files are named:
- p-1_40x20x20_case2.mat
For the high pulsatility case, pressure files are named:
- p-1_40x20x20_case4.mat
CONTENTS OF HAEMODYNAMICS ZIP FILE
Haemodynamics_dt_0.0125.zip contains two folders:
- Results_40x20x20_dt=0.0125
- Results_50x25x25_dt=0.0125
Each folder contains:
- p-1_40x20x20.mat to p-55_40x20x20.mat (pressure fields)
- u-1_40x20x20.mat (x-velocity)
- v-1_40x20x20.mat (y-velocity)
- w-1_40x20x20.mat (z-velocity)
- x, y, z coordinate files
GIF FILES:
- 2D_pressure.gif
- 2D_velocity_streamlines.gif
- 3D_pressure.gif
- 3D_u_velocity.gif
MATLAB SCRIPTS:
- contours_2D.m
- contours_3D.m
- main_visuals.m
Additionally, the high pulsatility case also includes:
- mat_to_dat.m
HOW TO RUN THE HAEMODYNAMIC AND FSI VISUALIZATIONS
-
Unzip one simulation folder
-
Open MATLAB
-
Set the folder as the MATLAB working directory
-
Run the command:
main_visuals
All required .dat files must be present in the same directory as the script.
TAWSS AND OSI VISUALIZATION TOOL
TAWSS_OSI.zip contains three folders:
- High Pulsatility
- Low Pulsatility
- Rigid Wall
Each folder corresponds to one simulation case.
CONTENTS OF EACH TAWSS/OSI FOLDER
Each folder contains:
- tawss_osi_plot.m
- OSI_*.dat
- TAWSS_*.dat
- rho1.dat to rho28.dat
- theta1.dat
- x1.dat
File naming:
High Pulsatility:
- OSI_FSI HP.dat
- TAWSS_FSI HP.dat
Low Pulsatility:
- OSI_FSI LP.dat
- TAWSS_FSI LP.dat
Rigid Wall:
- OSI_Rigid Wall.dat
- TAWSS_Rigid Wall.dat
VARIABLE DEFINITIONS
- TAWSS = Time-Averaged Wall Shear Stress i
- OSI = Oscillating Shear Index
- x = Axial coordinate
- rho = Radial coordinate
- theta = Circumferential coordinate
HOW TO RUN TAWSS AND OSI VISUALIZATIONS
-
Open one subfolder (Rigid_Wall, Low_Pulsatility, or High_Pulsatility)
-
Ensure all .dat files are in the same directory
-
Run the MATLAB script:
tawss_osi_plot.m
SOFTWARE AND HARDWARE REQUIREMENTS
- MATLAB R2023a or newer
- Minimum recommended memory: 16 GB RAM
- Compatible with Windows, Linux, and macOS
DATA PROCESSING LEVELS
Raw Data:
- Full three-dimensional velocity and pressure fields
Processed Data:
- Cross-Sectional velocities and pressures
- TAWSS and OSI fields
Visualization Data:
- Animated GIF files
DATA REUSE AND APPLICATIONS
This dataset may be reused for:
- CFD solver benchmarking
- Fluid–structure interaction validation
- Biomedical flow modeling
- Machine learning datasets for flow prediction
- Graduate-level CFD education
SHARING AND ACCESS INFORMATION
Repository: Dryad Digital Repository
Derived From: Fully numerical simulations
Related Publication: Aneurysmal Haemodynamics: A Three-Dimensional Fluid–Structure Interaction Approach
LEGAL AND ETHICAL STATEMENT
- No human subject or patient-specific data are included
- All data are fully anonymized and synthetic
- Data are released under the Creative Commons CC0 Public Domain dedication
END OF README