NUEN 489 GEMIX Project

Python script to run interFoam in parallel. Does blockMesh and decomposePar before running and reconstructPar after running. Has command line options for

• parallel decomposition (how many processors in x, y, z)
• number of volumes (split into x_inlet, x_mixing, y and z)
• turbulence model (k-Epsilon or k-Omega)
• switch for running in 3D

See ./run -h for help.

Contains fvSchemes, fvSolution, turbulenceProperties files corresponding to k-Epsilon or k-Omega. ./run copies them into their correct place in system or constant before running interFoam.

Location of blockMeshDict. Top is the top half of the blockMeshDict that is common to both 3D and 2D runs. 2D and 3D change the boundary conditions to be appropriate for 2 and 3 D geometry. simple is the blockMeshDict for 3D geometry without angled inlets.

./run alters local variables Nx1, Nx2, Ny, Nz which control the number of volumes used, concatenates top and 2D/3D and writes to system/blockMeshDict.

Plots the time step output to make sure the adjustable time step isn't going too small. Run as ./plotTime.py or python3 plotTime

Reads in a CSV file of the axial pressure to calculate the fRe product.

• boundary conditions for omega
• boundary conditions for k
• inlet conditions from data
• top and bottom inlets not mixing
• flow is currently not turbulent (Re ~ 50 for 1 m/s velocity and average viscosity)