wyllduck / aerospike Goto Github PK

OpenFoam (v9) - Aerospike CFD

How to Run It?

IMPORTANT: By default the UNV file has already been imported to OpenFOAM, thus you can skip to step 3)

0. ./Allclean
1. UNV file to OpenFoam ideasUnvToFoam meshes/aerospike_refined_2.unv
2. Go to polyMesh/boundaries and change the types!

• walls -> type: wall
• front -> type: empty
• back -> type: empty
• symmetry -> type: symmetry

3. Field initialization to help the solver setFields
4. Decompose in 4 parts (can do more or less depending on your setup) DecomposePar
5. Run the solver in PARALLEL ("> log" saves the output to a file named log) mpirun -np 4 rhoCentralFoam -parallel > log
6. ReconstructPar
7. paraFoam

OTPIONAL: You can check the residuals (in post-processing or on real-time) using the bash files ./plot_k_omega.sh ./plot_U.sh

Parameters simulation

• 1 atm - 30 bar chamber

Co 0.15 endTime 0.1800

• 0.6 atm - 30 bar chamber

Co 0.15 endTime 0.2050

NOTE: The rest of the parameters are keept the same.

Get Results Simulation + 3D model

(public access) https://1drv.ms/u/s!AhNN3_xNz6xdg-NUF0unrIyeq6NvCg?e=idaQrV

If you want to see the CFD results, download the .zip files and paste and replace all the files inside into this git repository. Then run paraFoam.

SOURCES:

Very useful introduction https://maplekeylabs.com/understanding-the-openfoam-case-structure/

KOmegaSST https://www.openfoam.com/documentation/guides/latest/doc/guide-turbulence-ras-k-omega-sst.html#sec-turbulence-ras-k-omega-sst-initialisation

KOmegeSST Setup Boundary Conditions (IMPORTANT) https://curiosityfluids.com/2019/02/11/basic-komega-sst-boundary-conditions/

nut: kinematic turbulence viscosity

omega: specific dissipation rate

k: turbulent kinetic energy

alphat: __________

Turbulence Kinematic Energy - Equation Explained https://en.wikipedia.org/wiki/Turbulence_kinetic_energy

Decompose Video https://www.youtube.com/watch?v=bHMdh_l45M0&t=1s