In collaboration with Miguel Angel Palomera-Pérez
Yguazú is a parallel, 3D, adaptive grid code developed to solve gasdynamic equations1 based on Van Leer's 'flux vector splitting' algorithm. It uses a fast, semi-implicit method for solving a system of rate equations for atomic/ionic and molecular species [RNV00]. The application was originally developed on a shared memory platform using Fortran and OpenMP as its programming environment. It integrates a 3D Euler (gasdynamic) system of equations and a series of rate equations. The system of equations is written as [RNV00]:
In this equation, the variables U, F, G, H and S are expressed in vector terms as [RNV00]:
In all these equations, the vector (u,v,w) represents the three components of flow velocity in the directions (x,y,z), ρ is density, P is pressure, T is temperature, and n1...nR are the number densities of the atomic/ionic and molecular species with masses m1...mR and charges z1...zR. The specific heat at constant volume (Cv) has a complex behavior, which for the purposes of this application is approximated ...