TY - GEN A1 - Oevermann, Michael A1 - Schmidt, Heiko A1 - Kerstein, Alan R. T1 - Linear-Eddy Modeling of Autoignition under Thermal Stratification N2 - The influence of thermal stratification on autoignition at constant volume and high pressure is investigated under turbulent conditions using the one-dimensional Linear-Eddy Model (LEM) and detailed hydrogen/air chemistry. Results are presented for the influence of initial temperature inhomogeneities on the heat release rate and the relative importance of diffusion and chemical reactions. The predicted heat release rates are compared with heat release rates of Chen et al. and Hawkes et al. obtained by two-dimensional Direct Numerical Simulations (DNS). Using the definition of Chen et al. for the displacement speed of the H2 mass fraction tracked at the location of maximum heat release, and a comparison of budget terms, different combustion modes including ignition front propagation and deflagration waves are identified and the results are compared to the DNS data. The LEM approach shows qualitatively and quantitatively reasonable agreement with the DNS data over the whole range of investigated temperature fluctuations. The results presented in this work suggest that LEM is a potential candidate as a sub-model for CFD calculations of HCCI engines. T3 - ZIB-Report - 08-05 KW - HCCI KW - Selbstzündung KW - Turbulenz KW - Linear-Eddy Modell KW - HCCI KW - autoignition KW - thermal explosion KW - turbulence KW - Linear-Eddy Model Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10562 SN - 1438-0064 ER - TY - GEN A1 - Oevermann, Michael A1 - Scharfenberg, Carsten A1 - Klein, Rupert T1 - A sharp interface finite volume method for elliptic equations on Cartesian grids N2 - We present a second order sharp interface finite volume method for the solution of the three-dimensional poisson equation with variable coefficients on Cartesian grids. In particular, we focus on interface problems with discontinuities in the coefficient, the source term, the solution, and the fluxes across the interface. The method uses standard piecewiese trilinear finite elements for normal cells and a double piecewise trilinear ansatz for the solution on cells intersected by the interface resulting always in a compact 27-point stencil. Singularities associated with vanishing partial volumes of intersected grid cells are removed by a two-term asymptotic approach. In contrast to the 2D method presented by two of the authors in [M.~Oevermann, R.~Klein: A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, J.~Comp.~Phys.~219 (2006)] we use a minimization technique to determine the unknown coefficients of the double trilinear ansatz. This simplifies the treatment of the different cut-cell types and avoids additional special operations for degenerated interface topologies. The resulting set of linear equations has been solved with a BiCGSTAB solver preconditioned with an algebraic multigrid. In various testcases -- including large coefficient ratios and non-smooth interfaces -- the method achieves second order of accuracy in the L_inf and L_2 norm. T3 - ZIB-Report - 08-38 KW - elliptische Gleichungen KW - Finite Volumen KW - eingebettete Ränder KW - variable und diskontinuierliche Koeffizienten KW - diskontinuierliche Lösung KW - Elliptic equations KW - finite volume methods KW - embedded interface KW - variable and discontinuous coefficients KW - discontinuous solution Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10900 SN - 1438-0064 ER - TY - GEN A1 - Oevermann, Michael A1 - Klein, Rupert T1 - A cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces N2 - We present a finite volume method for the solution of the two-dimensional Poisson equation $ \nabla\cdot( \beta( {\mbox{\boldmath $x$}}) \nabla u({\mbox{\boldmath $x$}})) = f(\mbox{\boldmath $x$}) $ with variable, discontinuous coefficients and solution discontinuities on irregular domains. The method uses bilinear ansatz functions on Cartesian grids for the solution $u({\mbox{\boldmath $x$})$ resulting in a compact nine-point stencil. The resulting linear problem has been solved with a standard multigrid solver. Singularities associated with vanishing partial volumes of intersected grid cells or the dual bilinear ansatz itself are removed by a two-step asymptotic approach. The method achieves second order of accuracy in the $L^\infty$ and $L^2$ norm. T3 - ZIB-Report - 06-05 KW - Poisson equation KW - finite volume methods KW - embedded interface KW - ariable and discontinuous coefficients KW - discontinuous solution Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8996 ER - TY - GEN A1 - Oevermann, Michael A1 - Klein, Rupert A1 - Berger, Marsha A1 - Goodman, Jonathan T1 - A Projection Method for Two-Phase Incompressible Flow with Surface Tension and Sharp Interface Resolution N2 - We present a fully second order projection method for the simulation of two-phase incompressible flow with surface tension. The Navier-Stokes equations are solved with a projection method on a fixed Cartesian grid. The free interface between the two fluids is tracked with a level set approach. The conditions at the interface for the pressure, the pressure gradient, and the velocity are explicitly incorporated into the scheme leading to a sharp representation of the pressure discontinuity and the interfacial force. The scheme in the presented form does not introduce additional points in the standard finite difference stencils. Computational results are compared with analytic solutions for a static round bubble, damped surface waves, and Rayleigh-Taylor instabilities. T3 - ZIB-Report - 00-17 KW - Two-phase flow KW - immersed interface method KW - incompressible flow KW - projection method KW - level set KW - multi grid Y1 - 2000 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5851 ER - TY - GEN A1 - Schmidt, Heiko A1 - Oevermann, Michael A1 - Bastiaans, Rob J.M. A1 - Kerstein, Alan R. T1 - A Priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds, Extended to Incorporate Strain Effects N2 - In this paper we revisit the a priori turbulent flame speed tabulation (TFST) technique for a given parameter space within the region of flamelet combustion-regimes. It can be used as a subgrid-scale (SGS) model in Large Eddy Simulation (LES). In a first step, stationary laminar flamelets are computed and stored over the progress variable following the ideas of flamelet generated manifolds (FGM). In a second step, the incompressible one-dimensional Navier-Stokes equations supplemented by the equation for the progress variable are solved on a grid that resolves all turbulent scales. Additionally, turbulent transport is implemented via the linear eddy model (LEM). The turbulent flame structures are solved until a statistically stationary mean value of the turbulent flame speed has been reached. The results are stored in a table that could be used by large scale premixed combustion models, e.g. front tracking schemes. First results are compared quantitatively with direct numerical simulations (DNS) taken from literature. Then it is illustrated in one example how the new method could help to fix constants in algebraic models for turbulent flame speeds. Further it is shown how the technique can be extended to incorporate turbulent strain effects. Finally we investigate the effect of the use of detailed and tabulated chemistry under unsteady conditions. T3 - ZIB-Report - 09-09 KW - Turbulente Verbrennung KW - Flammenstrukturen KW - Linear Eddy Model KW - FGM KW - Turbulente Brenngeschwindigkeit KW - Streckungseffekte KW - turbulent premixed combustion KW - flame structures KW - linear eddy model KW - flamelet generated manifolds KW - turbulent flame speed tabulation KW - strain effects Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11214 SN - 1438-0064 ER - TY - GEN A1 - Schmidt, Heiko A1 - Oevermann, Michael A1 - Bastiaans, Rob J.M. A1 - Kerstein, Alan R. T1 - A priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds N2 - In this paper we propose a technique for a priori turbulent flame speed tabulation (TFST) for a given parameter space in standard combustion-regime diagrams. It can be used as a subgrid-scale (SGS) model in Large Eddy Simulation (LES). In a first step, stationary laminar flamelets are computed and stored over the progress variable following the ideas of flamelet generated manifolds (FGM). In a second step, the incompressible one-dimensional Navier-Stokes equations supplemented by the equation for the progress variable are solved on a grid that resolves all turbulent scales. Additionally, turbulent transport is implemented via the linear eddy model (LEM). The turbulent flame structures are solved until a statistically stationary mean value of the turbulent flame speed has been reached. The results are stored in a table that could be used by large scale premixed combustion models, e.g. front tracking schemes. Results are compared to an algebraic model and to direct numerical simulations (DNS). T3 - ZIB-Report - 08-02 KW - Verbrennung KW - Flammenstrukturen KW - stochastisches Turbulenzmodell KW - FGM KW - turbulente Brenngeschwindigkeiten KW - turbulent premixed combustion KW - flame structures KW - linear eddy model KW - flamelet generated manifolds KW - turbulent flame speeds Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10532 SN - 1438-0064 ER -