gdbOF: A debugging tool for OpenFOAM®

Abstract

OpenFOAMR libraries are a great contribution to CFD community and a powerful way to create solvers and other tools. Nevertheless in this creative process a deep knowledge is needed concerning with classes structure, for value storage in geometric fields and also for matrices resulting from equation systems, becoming a hard task for debugging. To help in this process a new tool, called gdbOF, attachable to gdb (GNU Debugger) is presented in this paper. It allows to analyze classes structures at debugging time. This application is implemented by gdb macros, these macros can access to code classes and also to their data in a transparent way, giving the requested information. This tool is tested for different application cases, such as the assemble and storage of matrices in a scalar advectivediffusive problem, non orthogonal correction methods in purely diffusive tests and multiphase solvers based on Volume of Fluid Method. In these tests several type of data are checked, such as: internal and boundary vector and scalar values from solution fields, fluxes in cell faces, boundary patches and boundary conditions. As additional features of this tool data dumping to file and a graphical monitoring of fields are presented. All these capabilities give to gdbOF a wide range of use not only in academic tests but also in real problems.

On a dissolution–diffusion model. Existence, uniqueness, regularity and simulations

“…One shaped as the A-A ′ cut and length 19 cm, the other one shaped as B-B ′ and length 10.37 cm. The agreement between our computations and the measurements is remarkable, and comparable to the computations on the three-dimensional domain reported in [18]; where meshes with almost four million elements were needed.…”

“…1). We recall that it is initially loaded with 1.9 g of Progesterone, and is used by cattle producers to achieve estrous synchronization, which in turn allow a better outcome of artificial insemination; see [18] and references therein.…”

“…The parameter values used for the simulation of the drug release from this device, taken from [18], are stated in Table 1; C 0 denotes the initial concentration of dissolved drug and m 0,p , a 0,p stand for the initial mass and area of each particle of solid drug, respectively, whereas Num is the total number of particles in the device, which we assumed evenly distributed. The diffusion coefficients D ap40 and D ap60 correspond to a liberation media with 40% and 60% ethanol in water, respectively, and k D is the dissolution rate of the particles of solid drug.…”

“…We observe that the computations according to this proposal fit the experimental data very well. Visually, the fit is equal to the one obtained in [18], where three-dimensional computations, on a mesh with 3 816 762 elements were performed. It is interesting to note that adding the contributions of the two bi-dimensional meshes that we used, we arrive at 75 280 elements and 39 031 vertices.…”

An extended mixture model for the simultaneous treatment of small‐scale and large‐scale interfaces

“…From these steps, it is clear that the correct calculation of the flux and its treatment along the solving algorithm plays a central role in the successful implementation of the solver. Because the presented algorithm is not available in the OpenFOAM ® suite, it was programmed starting from VOF solver and using the gdbOF debugging tools .…”