ParOSol TU Wien

Introduction

ParOSol stands for PARallel Octree SOLver.

It is a fully-parallel micro-FE analysis code based on Octree to solve linear elasticity problems. Because it works on equal sized voxel geometry, it solves the problem direct on 3D images, which can be obtained by CT scans.

ParOSol TU Wien is a fork of ParOSol. The following changes have been made:

1. The code now works with a consistent system of units (m and Pa or mm and MPa).
2. The Poisson ratio that is read from the input HDF5 data is used, not the hard coded one.
3. The ParOSol output is formatted like ParFE output (but only to the level needed by the postprocessing tools used at ILSB TU Wien.)

Note that the data format used by the original ParOSol is not compatible with this version! You can use the .h5 as well, but they need to be patched first.

The changes in code have been tested well at the ILSB TU Wien, but are still provided without warranty of any kind. Use at your own discretion.

Licence

Refer to ParOSol licensing.

Building

ParOSol was tested on Debian Bookworm, Ubuntu 22.04 and Ubuntu 24.04.

To build on Debian based rystems run the following commands:

apt update
apt install -y git build-essential cmake libopenmpi-dev libhdf5-openmpi-dev libeigen3-dev
git clone https://github.com/reox/parosol-tu-wien.git
cd parosol-tu-wien
mkdir build
cd build
cmake .. -DCMAKE_BUILD_TYPE=Release  # Note: Building "Release" brings a huge speedup!
make
cd ..

Then you can run the example meshes:

tar xvf mesh.tar.bz2
mpirun -np 8 build/parosol mesh/h5/sphere.h5

(Adjust -np to the number of processors you have)

And look at the results with paraview. NOTE: Select XDMF Reader!

python3 tools/createxmf.py mesh/h5/sphere.h5
paraview mesh/h5/sphere.xmf