Coupled solution of the species conservation equations using unstructured finite-volume method

2009 ◽  
Vol 64 (4) ◽  
pp. 409-442 ◽  
Author(s):  
Ankan Kumar ◽  
Sandip Mazumder
2008 ◽  
Author(s):  
Ankan Kumar ◽  
Sandip Mazumder

Many reacting flow applications mandate coupled solution of the species conservation equations. A low-memory coupled solver was developed to solve the species transport equations on an unstructured mesh with implicit spatial as well as species-to-species coupling. First, the computational domain was decomposed into sub-domains comprised of geometrically contiguous cells—a process termed internal domain decomposition (IDD). This was done using the binary spatial partitioning (BSP) algorithm. Following this step, for each sub-domain, the discretized equations were developed using the finite-volume method, written in block implicit form, and solved using an iterative solver based on Krylov sub-space iterations, i.e., the Generalized Minimum Residual (GMRES) solver. Overall (outer) iterations were then performed to treat explicitness at sub-domain interfaces and non-linearities in the governing equations. The solver is demonstrated for a laminar ethane-air flame calculation with five species and a single reaction step, and for a catalytic methane-air combustion case with 19 species and 22 reaction steps. It was found that the best performance is manifested for sub-domain size of about 1000 cells, the exact number depending on the problem at hand. The overall gain in computational efficiency was found to be a factor of 2–5 over the block Gauss-Seidel procedure.


2015 ◽  
Vol 31 (3) ◽  
pp. 937-950 ◽  
Author(s):  
Perry L. Johnson ◽  
Jared M. Pent ◽  
Hrvoje Jasak ◽  
J. Enrique Portillo

Author(s):  
Dong Jin Kang ◽  
Sang Soo Bae ◽  
Jae Won Kim

A Navier-Stokes simulation of the MIT flapping foil experiment is presented. The MIT experiment was designed to provide a good quality database for unsteady boundary layer flows. The unsteady boundary layer around a hydrofoil was generated by flapping two airfoils upstream of the hydrofoil. Present Navier-Stokes simulation is carried out on the entire experimental domain, including the flapping airfoils as well as the downstream fixed hydrofoil. Present Navier-Stokes code uses an unstructured finite volume method based on the SIMPLE algorithm. It uses QUICK scheme for the convective terms and the second order Euler backward differencing for time derivatives to keep second order accuracy spatially and temporally. All other spatial derivatives are approximated by using central difference scheme. All comparisons of present time averaged and unsteady solutions with the corresponding experimental data are satisfactory: all unsteady solutions are compared in terms of time mean and first harmonic. The first harmonic of the velocity shows a peak inside the boundary layer along the surfaces of the hydrofoil and has a local minimum near the edge of the boundary layer. The local minimum becomes manifest as the boundary layer grows. The unsteadiness in the free stream is transferred inside the boundary layer when an unsteady vortex impinges on the surface. The entrained unsteadiness travels with a local velocity slower than that in the free stream. This causes phase lag of the first harmonic between the free stream and the boundary layer and local minimum of the first harmonic near the edge of the boundary layer.


Sign in / Sign up

Export Citation Format

Share Document