An integrated code-base in the existing OpenFoam framework will be developed, such that an "UGent OpenFoam Community" can emerge, taking advantage of code sharing and joint code development. The goal of this community is to enhance collaboration between different research groups, faculties, and scientific domains. This will result in the development of an open source CFD code-base in the OpenFoam framework. The contributed implementations can be read by all members of the community, and in this way, the OpenFoam CFD research will contribute to a more reproducible, peer reviewed way of research. Further, in the open source spirit, issues can be tackled collaboratively so that researchers can contribute their specific expertise to a broad variety of research questions.
A link to the repository will be posted here as soon as a stable version is available (most likely within 2017). For more information on having access to the repository please see the Contact page.
An overview of what is currently available in the repository is listed below:
Combustion
— EDM: Eddy Dissipation Model with a reaction time scale for under-resolved fire dynamics based on [1].
— EDC: Eddy Dissipation Concept based on [2]-[3].
Turbulence
— constantSmagorinsky: The constant Smagorinsky model for compressible flows [4]. Requires the input of the constant cs. The additional coefficients ck and ce (cs^2 = ck*sqrt(ck/ce) [5]) are used to calculate the sub-grid scale kinetic energy k_sgs.
— dynamicSmagorinsky: The dynamic Smagorinsky model for compressible flows [6].
— dynamicSmagorinskyVariablePrt: The dynamic Smagorinsky model for compressible flows [6] with a variable Prandtl number formulation [7].
— oneEqEddyBuoyancy: The one-equation Eddy Viscosity Model for compressible flows with buoyancy production term.
Wall functions
— alphaSgsConvectiveWallFunction: It sets the convective heat flux to be q_conv=\alpha_eff*(T_gas-T_wall). Calculates the convective heat transfer coefficient based on correlations for the Nusselt number. The characteristic length, L, of the configuration is required.
— muSgsConvectiveWallFunction: This boundary condition provides a turbulent dynamic viscosity condition at the walls as mu_sgs=alpha_sgs*Pr_t.
Sovers
— icoFoam_convergence: The icoFoam convergence solver has the same solver properties as the original icoFoam solver: the "core" of the solver is not edited. This solver checks for convergence by means of maximum residuals or maximum averaged squared residuals. The latter is added as a commented code block. If the maximum residual is less then a specified maximum (described in the casedir/constant/convergenceProperties file) on a specific time step, then the runTime-loop will write out the fields at that timestep and end the runtime loop prematurely. This solver is used for obtaining a steady state velocity and pressure profile.
— icoFoam_pingpong: Laminar transient flow OpenFoam solver for solutes subjected by ping pong reaction kinetics. Only half-reactions are implemented. No steady state / reverse reaction. Dictionary solute: make a list of solutes: E, SA, SB, PP, PQ E = enzyme; SA = substrate A; SB = substrate B; PP = product P (from substrate A); PQ = product Q (from substrate B); All the necessary constants are written in this dictionary (diffusion constant, diffusion model, molar mass). Dictionary reaction: reads the necessary reaction kinetic parameters for the ping pong kinetics (catalyst constants, Michaelis Menten constants, core inhibition constants). Member functions are included to enable parameters to be read in the icoFoam_solute_pingpong.C file.
— icoFoam_pingpong_enzymeplacement: This solver is a merge of the enzyme placement utility and the icoFoam ping pong solver.
— PPBB_SS: This solver calculate the scalar transport equations for various solutes coupled with enzyme kinetics. The solver utilises a predefined steady state flow and pressure field.
— pingpong: This solver calculate the scalar transport equations for various solutes coupled with enzyme kinetics. The solver utilises a predefined steady state flow and pressure field.
— pingpong_SS_testing: This solver calculate the scalar transport equations for various solutes coupled with enzyme kinetics. The solver utilises a predefined steady state flow and pressure field. Aims at obtaining steady state concentration in the reactor.
— pingpong_scalarTransport: Scalar transport of the solutes of the pingpong reaction. NO KINETICS.
— pingpong_testing: This solver calculates the scalar transport equations for various solutes coupled with enzyme kinetics. The solver utilises a predefined steady state flow and pressure field.
Utilities
— enzyme_placement_wall (pre-processing): Identifies cells adjacent to the wall in order to set adapt the values for those cells, to mimic a reactor with immobilized enzyme.
Tutorials
— FireFOAM: 1DPyrolysis, compartmentFire, firePlume, flameSpread_corner, flameSpread_parallelPanels, heliumPlume, thermalPlume, tunnelFire.
— icoFoam: A variaty of tutorials based on the solvers mentioned above.
References
[1]: R. McDermott, K. McGrattan, and J. Floyd. A simple reaction time scale for under-resolved fire dynamics. In Fire Safety Science - Proceedings of the 10th International Symposium, pages 809-820, University of Maryland, College Park, Maryland, USA, 2011, p 45.
[2]: B. Panjwani, I. Ertesvag, K.E. Rian and A. Gruber, Subgrid Combustion modeling for large eddy simulation (LES) of turbulent combustion using eddy dissipation concept, V European Conference on Computational Fluid Dynamics (ECCOMAS), Lisbon, Portugal, 14-17 June 2010.
[3]: Z. Chen, J. Wen, B. Xu, S. Dembele, Large eddy simulation of a medium-scale methanol pool fire using the extended eddy dissipation concept, International Journal of Heat and Mass Transfer 70 (2014) 389-408.
[4]: J. Smagorinsky, General Circulation Experiments with the Primitive Equations, Monthly Weather Review 91 (1963) 99-164.
[5]: P.P Sullivan, J.C. McWilliams, C.-H. Moeng, A subgrid-scale model for large-eddy simulation of planetary boundary layer, Bound.-Layer Meteor. 71 (1994) 274-276.
[6]: P. Moin, K. Squires, W. Cabot, and S. Lee. A dynamic subgrid-scale model for compressible turbulence and scalar transport. Phys. Fluids A 3 (1991) 2746-2757.
[7]: Martin, M.P., Piomelli, U. and Candler, G.V., Sub-grid-Scale Models for Compressible Large-Eddy Simulations, Theoret. Comput. Fluid Dynamics 13 (2000) 361-276.