Two-Dimensional Approaches to Sparse Matrix Partitioning

Most of the equations used to describe the behaviour of continua are of the form:u,t +∇ · (Fa − Fv) = S(u) , (9.1) where u,Fa,Fv,S(u) denote the vector of unknowns, advective and diffusive flux tensors, and source-terms respectively. In the case of compressible gases, we haveu = (ρ ; ρvi ; ρe) ,Faj = (ρvj ; ρvivj + pδij ; vj(ρe+ p)) ,Fvj = (0 ; σij ; vlσlj + kT,j) . (9.2)Here ρ, p, e, T, k, vi denote the density, pressure, specific total energy, temperature, conductivity and fluid velocity in direction xi respectively. This set of equations is closed by providing an equation of state, e.g., for a polytropic gas:p = (γ − 1)ρ[e− 1 2 vjvj ] , T = cv[e − 12 vjvj ] , (9.3 a, b)where γ, cv are the ratio of specific heats and the specific heat at constant volume respectively. Furthermore, the relationship between the stress tensor σij and the deformation rate must be supplied. For water and almost all gases, Newton’s hypothesisσij = µ( ∂vi ∂xj+ ∂vj ∂xi) + λ ∂vk ∂xkδij (9.4)complemented with Stokes’ hypothesisλ = −2µ 3(9.5)is an excellent approximation. The compressible Euler equations are obtained by neglecting the viscous fluxes, i.e., setting Fv = 0. The incompressible Euler or Navier-Stokes equations are obtained by assuming that the density is constant and that pressure does not depend on temperature. The Maxwell equations of electromagnetics, the heat conduction equations of solids, and the equations describing elastic solids undergoing small deformation can readily be written in the form given by Equation (9.1).