DOI: 10.5061/dryad.pc866t22m
title: Aneurysmal haemodynamics: A three-dimensional fluid-structure interaction approach

description: 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.
description: # Aneurysmal haemodynamics: A three-dimensional fluid-structure
 interaction approach 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
 --- 1. Unzip one simulation folder 2. Open MATLAB 3. Set the folder as the
 MATLAB working directory 4. 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 --- 1. Open one subfolder (Rigid_Wall,
 Low_Pulsatility, or High_Pulsatility) 2. Ensure all .dat files are in the
 same directory 3. 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