In March 2021, CFD Direct replaced fvOptions with new fvModels and fvConstraints classes. fvModels provided a more complete framework for implementation of complex models which could be optionally selected during a simulation. New fvModels were implemented for heat transfer and sources, which are easier to configure than a general fvOption. Pressure limiting was implemented as an fvConstraint named limitPressure, consistent with limitVelocity and limitTemperature. With additional fvModels such as cloud and surfaceFilm, solvers such as reactingFoam replicate the functionality of several existing solvers which are deprecated.
In 2020, CFD Direct were tasked to develop multicomponent diffusion modelling in OpenFOAM with funding from the Process Engineering Consortium. The Consortium need models, e.g. Maxwell-Stefan, that represent complex diffusion of species in fluid mixtures, affecting both specie concentration and energy. As part of the work, we created a new ThermophysicalTransportModels library with a new interfaces for: fluxes of heat, q() and divq(); and fluxes of mass, j() and divj(). As a consequence, heat transfer became consistent across applications, boundary conditions and data processing. This critical code redesign is typical of our strategy for sustainable development of OpenFOAM, supported by dedicated maintenance funding from OpenFOAM supporters.
In July 2020, CFD Direct released a framework for surface interpolation in interface capturing for multiphase flows in OpenFOAM. All interpolation schemes, used for flux calculation in advection, operate with the MULES-based solvers, e.g. interFoam. The framework is easy to use, extensible, and robust due to bounded, conservative numerics. The new framework comes with a new family of interpolation schemes based on piecewise-linear interface calculation (PLIC). The PLIC schemes provide a promising alternative to standard interface compression, performing better for meshes with regions of cell refinement. All schemes are conservative, bounded and robust and can run efficiently with large time steps.
In May 2017, we implemented water wave modelling into the public, development line of OpenFOAM (OpenFOAM-dev), to simulate flow conditions experienced in marine and offshore engineering, e.g. on a ship in open water. The initial implementation includes the first order Airy wave and the second order Stokes wave (Stokes II), a second-order accurate utility to apply waves for case initialisation to the entire domain. A verticalDamping source (fvOption) provides wave damping towards outlet regions and an interfaceHeight function object processes interface elevation at specified locations.
In April 2017, we added barycentric tracking into the public, development line of OpenFOAM (OpenFOAM-dev), as a complete replacement to the tracking algorithm that existed in OpenFOAM for over 10 years. Barycentric tracking works on any decomposed tet mesh, irrespective of mesh quality, including poor quality, flat and inverted tetrahedra. The new development was funded from the €100 k for OpenFOAM maintenance, raised through the OpenFOAM Foundation in 2017, in which “particles and tracking” was identified as an area in OpenFOAM “requiring significant code refactoring and/or rewriting”.
OpenFOAM is free software, meaning users have the freedom to run, copy, distribute, study, change and improve the software. Users benefit from these freedoms, which account for much of OpenFOAM’s popularity. The OpenFOAM Foundation is the copyright holder of OpenFOAM, which it licenses exclusively under the GPL. It maintains a strong legal position to enforce the licence and preserve its freedoms, by being the single owner of OpenFOAM. This requires contributors to the project to assign copyright in their OpenFOAM contributions to it, through its Contributor Agreement. Organisations with a serious commitment to free software are signing the Agreement, including CFD Direct, blueCAPE, VTT Technical Research of Finland Ltd and Intel.