Towards Accommodating Comprehensive Chemical Reaction Mechanisms in Practical Internal Combustion Engine Simulations
|
|
- Thomas Jacobs
- 6 years ago
- Views:
Transcription
1 Paper # 0701C-0326 Topic: Internal Combustion and gas turbine engines 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 Towards Accommodating Comprehensive Chemical Reaction Mechanisms in Practical Internal Combustion Engine Simulations Mandhapati Raju 1 Mingjie Wang 1 Meizhong Dai 1 Shaoping Quan 1 Peter Kelly Senecal 1 Sibendu Som 2 Matt McNenly 3 Daniel L Flowers 3 1 Convergent Science Inc., Middleton, WI 53562, USA 2 Argonne National Laboratory, Argonne, IL 60439, USA 3 Lawrence Livermore National Laboratory, Livermore, CA 94551, USA The use of detailed chemical mechanisms, although preferable, poses an enormous challenge for practical engine simulations. To accommodate the use of comprehensive detailed mechanisms various acceleration strategies have been implemented into the CONVERGE CFD solver. The different strategies involve (a) dynamic mechanism reduction, (b) adaptive zoning, (c) use of a preconditioned iterative ODE solver, and (d) strategies for reducing the CPU time needed for species transport. The feasibility of using comprehensive detailed reaction mechanisms in practical engine simulations is tested. As a test case, a single-cylinder Diesel engine is simulated using a detailed mechanism for n-heptane developed by Lawrence Livermore National Laboratory (LLNL). The reduction in CPU time for each of the reduction strategies along with the accuracy of the simulated results is reported. 1. Introduction Using detailed chemical reactions mechanisms for engine and gas turbine applications is formidable. Typically mechanisms consisting of less than 50 species or so are used to simulate within a reasonable amount of time. However, the sizes of the detailed chemical reaction mechanism are quite large and impractical for three dimensional simulations [Lu and Law, 2009]. Different strategies are being applied by researchers to reduce the computational time significantly. The different strategies for reducing the simulation time are (i) reducing the number of combustion calls using an adaptive zoning approach [Aceves et al. 2000; Babajimopoulos et al., 2005; Raju et al., 2012a] (ii) reduce the size of system using a pre-reduced mechanism [Raju et al., 2012b] or by using dynamic mechanism reduction [Raju et al., 2013] (iii) reduce the computational time of each combustion call using methods like analytical evaluation of Jacobian, exploiting the sparsity of the Jacobian matrix and using special techniques like preconditioned iterative solver (as will be demonstrated in this paper). When large mechanism is involved, the transport of all the species could also significantly increase the computational time. New strategies have to be developed to reduce the time for transport of all the species. This paper demonstrates the combination of all the above strategies for engine simulations giving the capability to run full detailed chemical reaction mechanisms within reasonable computational times. Some of the strategies like adaptive zoning and dynamic mechanism reduction were already demonstrated elsewhere. The use of preconditioned iterative solver is a novel contribution and hence the advantage of using it is first demonstrated for a HCCI case. This is followed by the results for a diesel engine sector simulation using a detailed n-heptane mechanism from Lawrence Livermore National Laboratory. 2. Acceleration strategies (a) Adaptive zoning Adaptive zoning is used to accelerate the chemistry solver by reducing the number of calls to the chemistry solver. Instead of invoking chemistry at each and every cell, cells of similar physical and thermodynamic properties are grouped into zones. At each discrete time t, the cells are grouped in zones based on the physical and thermodynamic state of the
2 cell. In a two-dimensional zoning strategy, the zoning is done based on two variables, the temperature and the equivalence ratio of the cells. For each zone, the average temperature, pressure and the composition are determined to specify the thermodynamic state of the mixture in that zone. The SAGE chemistry solver is invoked on each zone. The zonal solution is then remapped onto the individual cells. The zoning strategy and remapping strategy are discussed in detail in Raju et al (b) Dynamic Mechanism Reduction Dynamic mechanism reduction involves the reduction of detailed chemical reaction mechanism during run time for each cell. The state of the art Directed relation graph with error propagation (DRGEP) reduction technique [Pepiot and Pitsch, 2008; Niemeyer et al., 2010] is used to generate locally reduced mechanisms for each cell and for each of the CFD time step. Since the local mechanisms are much smaller in size in comparison to the full size detailed mechanism, significant savings in computational time is achieved. Locally reduced mechanisms are created at each cell/zone at each time step. Chemistry solver is invoked in based on the locally reduced mechanism. The production/consumption rate of the unimportant species is taken as zero. There is some overhead for generating the reduced mechanism. The details of the implementation are presented in Raju et al (c) Preconditioned Iterative solver The chemistry solver uses CVODES solver which is part of the public domain Sundials package developed by Lawrence Livermore National Laboratory [Hindmarsh et al., 2005; Serban and Hindmarsh, 2005]. When a mechanism consists of large number of species, the size of the Jacobian becomes very large. The factorization of the Jacobian matrix in a direct solver becomes very expensive and often leads to a high CPU time for solving the system of ODE equations. CVODES solver provides an option to use iterative solver. However, iterative solver fails miserably without a suitable preconditioner. In this work an incomplete LU (ILU) factorization is used as a preconditioner for the original Jacobian matrix. Advantage is taken of the sparsity of the original Jacobian and the preconditioned matrix. The following section describes briefly the methodology of ODE integration using iterative solver in CVODES. The chemistry solver solves the species and the energy equation at each and every cell/zone. The system of ODE equations can be represented as y f t, y, y t0 y0. (1) For stiff systems, the CVODES solver uses the backward differentiation formula (BDF) method, which results in a system of non-linear equations at each time step represented as n n G( yn) y hn n,0 f tn, y an 0 (2) nm 1 H y y nm G y nm (3) H I J, hn n,0 (4) f J y (5) where hn tn tn 1 is step size, a n and n,0 are the coefficients. The assembly of the jacobian matrix J is analytical evaluation of the derivatives which results in a significant computational savings as compared to the conventional evaluation of derivative using numerical finite difference approximation. The results system of linear equations can be solved either using a direct solver or an iterative solver. As mentioned earlier, the use of iterative solver necessitates provision of a good preconditioner. The preconditioned system of linear equations can be written as 1 M H y M 1 G (6) The preconditioner M is approximate H in some sense. This preconditioner is called a left preconditioner as the preconditioning is applied left of the LHS matrix. The preconditioning M is obtained by dropping some of the elements of the original H matrix which are below a certain tolerance level ε. This results in increasing the sparsity of the M matrix. Both the matrix H and the matrix M are stored in sparse format to take advantage of the sparsity. In-house sparse Gaussian elimination routine is used to factorize the matrix M. Being sparser than the original matrix H, Matrix M would take lesser time for factorization and hence is suitable as a good preconditioner. GMRES is used as the iterative solver for solving the system of equations. The value of tolerance level ε is hardcoded as 0.1. Further study is required for determining the proper choice of value for ε. 2
3 (d) Strategies for reducing the time for transport of the species The above strategies (a-c) deal with the reduction of the chemistry computational time. Typically, the chemistry time is the bottleneck for the internal combustion (IC) engine simulations, especially for large size mechanisms. When the chemistry time is significantly reduced by the above strategies, the time for the transport of the species becomes significant. For all the above strategies, the transport of all the species is inevitable. Even though the mechanism is reduced locally, all the species need to be transported. Earlier researchers [Ren et al., 2011] have proposed RCCI type algorithms to reduce the number of species to be transported. In this paper, only a simple strategy is used to reduce the number of transported species. Transport of a species is skipped if the maximum mass fraction of the species in the entire domain is less than a tolerance level (typically 1.e-14). The species mass fractions are renormalized at every time step to avoid losing any mass due to the above approximation. This approximation gives slight reduction in the time required for transport of the species. Ongoing research is still being conducted to add in more robust algorithm to reduce the time for transport of the species. 3. Results and Discussion The strategies (a-c) accelerate the combustion simulations and strategy (d) accelerates the calculation of the transport of the species. The efficiency of strategy (a) has been demonstrated in an earlier study [Raju et al. 2012] and that of strategy (b) has been demonstrated in another study [Raju et al. 2013]. The primary purpose of this study is to first demonstrate the efficiency of the iterative preconditioned solver and then to demonstrate the combination of all the strategies for a real engine simulation. To test the computational efficiency of the iterative preconditioner solver, a HCCI simulation is performed for n- heptane combustion using a full LLNL n-heptane mechanism consisting of 654 species. Two cases are run, one with a dense solver and the other using an iterative solver. Adaptive zoning is used for both the cases as adaptive zoning significantly reduces the CPU time for combustion. There is a negligible difference in the solutions between the two cases (not shown here) as would be expected. Table 1 shows the CPU time comparison for both the cases. For the given case, the combustion time is only a small portion of the total CPU time. Due to the presence of large number of species, the transport time dominates the simulation. A speed up of 13 times is obtained for the combustion simulation using an iterative preconditioned solver. Table 1: Comparison of CPU Times for Dense Solver and Iterative Preconditioned Solver (654 species) Case CPU time for combustion (sec) Total CPU time (sec) Dense solver Iterative preconditioned solver The efficiency of the iterative preconditioned solver using a large size mechanism has been demonstrated using a HCCI case. In a real engine simulation, the time for combustion dominates the total simulation time. To demonstrate the use of a large size mechanism for engine simulation, a diesel sector simulation using detailed n-heptane mechanism is performed. Table 2 shows the engine parameters. For large size mechanisms (> 500 species or so), the use of SAGE (chemistry on each cell) is prohibitively expensive and is not attempted. Adaptive zoning is used by default. Three different simulations are performed (i) only strategy (a), referred as AZ, (ii) strategies (a) and (b), referred as AZ-IP, (iii) strategies (a),(b), (c) and (d), referred as AZ-IP-DMR. For dynamic mechanism reduction, the DRGEP tolerance value is chosen as The chosen target species are n-c7h16, OH and N2. Figure 1 shows the comparison of the average cylinder pressure for cases (ii) and (iii). The results for case (i) and case (ii) are exactly the same and hence case (i) is not shown in the figure. Figure 1 shows the pressure profiles for both case (ii) and case (iii) are close to each other. The peak pressure for case (iii) is slightly less than that for case (ii). The reason for the differences is due to the approximations involved in mechanism reduction on the fly. Table 2 shows the comparison of the computational time for cases (i), (ii) and (iii) on 12 cores. Table 2 shows a huge savings in the computational time using the iterative preconditioned solver (speed up of ~7.5 times). A further speed up of 4 times is obtained by using dynamic mechanism reduction. The overall computational time is reduced from 2 months to 2 days using all the above mentioned strategies. It is to be noted that without the use of adaptive zoning strategy, the simulation is prohibitively expensive. It is projected that it would take around 6 months and so is not attempted here. The use of all 3
4 the acceleration strategies has significantly reduced the total computational time that it is now possible to use full size mechanism with a little loss in accuracy Average Pressure [MPa] AZ-IP AZ-IP-DMR CA Figure 1: Comparison of average cylinder pressure as a function of crank angle for case (ii) and case (iii) Heat release rate [J/m 3 s] AZ-IP AZ-IP-DMR CA Figure 2: Comparison of heat release rate as a function of crank angle for case (ii) and case (iii) 4
5 Table 2: Comparison of the total CPU time for the three different cases Case Set up Total CPU time AZ AZ-IP AZ-IP-DMR Over 60 days 8 days 2 days 4. Conclusions Different strategies for accelerating combustion simulations towards accommodating large size chemical mechanisms for engine simulations have been proposed. The first strategy is to use adaptive zoning to reduce the number of calls to the combustion calculations. The second strategy is to use dynamic mechanism reduction, in which each cell/zone performs combustion simulations on a locally reduced chemical reaction mechanism. The third strategy is to use iterative preconditioned based ODE solver which takes advantage of the sparsity of the jacobian matrix and incomplete LU based preconditioner to accelerate the ODE solver. The fourth strategy is to skip the transport of insignificant species. It is found that for large size mechanism, the use of iterative preconditioned solver itself gives a significant speedup. For a HCCI simulation with 654 species mechanism, it is found that the time for combustion simulation is 13 times faster. With the use of all the above mentioned strategies, it is found that for an engine sector simulation with 654 species, the simulation can be completed within a couple of days, making it possible to use large size mechanisms for practical engine simulations within a reasonable amount of simulation time. Acknowledgements The authors acknowledge the support of Convergent Science Inc. References Aceves, S. M. Flowers, D. L., Westbrook, C., Smith, J. et al., "A Multi-Zone Model for Prediction of HCCI Combustion and Emissions," SAE Technical Paper Babajimopoulos, A., Assanis, D, N., Flowers, D. L., Aceves, S. M., Hessel, R. P, "A fully coupled computational fluid dynamics and multi-zone model with detailed chemical kinetics for the simulation of premixed charge compression ignition engines", International Journal of Engine Research, 6(5), pp , Hindmarsh, A. C., Brown, P. N., Grant, K. E., Lee, S. L., Serban, R., Shumaker, D. E. and Woodward, C. S. "SUNDIALS: Suite of Nonlinear and Differential/Algebraic Equation Solvers," ACM Transactions on Mathematical Software, 31(3), pp , Lu, T., Law, C.K., Toward accommodating realistic fuel chemistry in large-scale computations, Progress in Energy and Combustion Science, 35(2), pp , K.E. Niemeyer, C.-J. Sung, M.P. Raju, Skeletal mechanism generation for surrogate fuels using directed relation graph with error propagation and sensitivity analysis, Combustion and Flame, 157(9), pp , Pepiot-Desjardins, P. and Pitsch, H., An efficient error propagation based reduction method for large chemical kinetics mechanisms, Combustion and Flame, 154(1-2), pp , Raju, M., Wang, M., Dai, M., Piggott, W., Flowers, D. Acceleration of Detailed Chemical Kinetics Using Multi-zone Modeling for CFD in Internal Combustion Engine Simulations, SAE , 2012a. 5
6 Raju, M., Wang, M., Senecal, P.K., Som, S., Longman, D. E., "A Reduced Diesel Surrogate Mechanism For Compression Ignition Engine Applications," Proceedings of the ASME 2012 Internal Combustion Engine Division Fall Technical Conference, ICEF2012, September 23-26, Vancouver, BC, Canada, 2012b. Raju, M., Wang, M., Senecal, P.K., "Dynamic Chemical Mechanism Reduction for Internal Combustion Engine Simulations" SAE , Ren Z., Goldin, G.M., Hiremath, V., Pope, S. B., "Reduced description of reactive flows with tabulation of chemistry," Combustion Theory and Modelling, 15(6), pp , Serban, R., and Hindmarsh, A. C., "CVODES: the Sensitivity-Enabled ODE Solver in SUNDIALS," Proceedings of IDETC/CIE 2005, Sept. 2005, Long Beach, CA. 6
Adaptive Preconditioning Strategies for Integrating Large Kinetic Mechanisms
Paper # 070RK-0377 Topic: Reaction Kinetics 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More informationRESOLVING TURBULENCE- CHEMISTRY INTERACTIONS IN MIXING-CONTROLLED COMBUSTION WITH LES AND DETAILED CHEMISTRY
RESOLVING TURBULENCE- CHEMISTRY INTERACTIONS IN MIXING-CONTROLLED COMBUSTION WITH LES AND DETAILED CHEMISTRY Convergent Science White Paper COPYRIGHT 218 CONVERGENT SCIENCE. All rights reserved. OVERVIEW
More informationDetailed Chemical Kinetics in Multidimensional CFD Using Storage/Retrieval Algorithms
13 th International Multidimensional Engine Modeling User's Group Meeting, Detroit, MI (2 March 23) Detailed Chemical Kinetics in Multidimensional CFD Using Storage/Retrieval Algorithms D.C. Haworth, L.
More informationOptimizing Time Integration of Chemical-Kinetic Networks for Speed and Accuracy
Paper # 070RK-0363 Topic: Reaction Kinetics 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013
More informationThe Use of Dynamic Adaptive Chemistry and Tabulation in Reactive Flow Simulations
Paper # 070RK-0017 Topic: Reaction kinetics 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013
More informationMatrix Reduction Techniques for Ordinary Differential Equations in Chemical Systems
Matrix Reduction Techniques for Ordinary Differential Equations in Chemical Systems Varad Deshmukh University of California, Santa Barbara April 22, 2013 Contents 1 Introduction 3 2 Chemical Models 3 3
More informationDARS Digital Analysis of Reactive Systems
DARS Digital Analysis of Reactive Systems Introduction DARS is a complex chemical reaction analysis system, developed by DigAnaRS. Our latest version, DARS V2.0, was released in September 2008 and new
More informationValidation of a sparse analytical Jacobian chemistry solver for heavyduty diesel engine simulations with comprehensive reaction mechanisms
212-1-1974 Validation of a sparse analytical Jacobian chemistry solver for heavyduty diesel engine simulations with comprehensive reaction mechanisms Federico Perini 1,2, G. Cantore 1, E. Galligani 1,
More informationApplication of Model Fuels to Engine Simulation
Paper #E28 First published in: Topic: Engine 5 th US Combustion Meeting Organized by the Western States Section of the Combustion Institute and Hosted by the University of California at San Diego March
More informationHigh-dimensional, unsupervised cell clustering for computationally efficient engine simulations with detailed combustion chemistry
High-dimensional, unsupervised cell clustering for computationally efficient engine simulations with detailed combustion chemistry Federico Perini 1,a a Dipartimento di Ingegneria Meccanica e Civile, Università
More informationTarget Species Selection for Dynamic Adaptive Chemistry Simulations
University of Connecticut DigitalCommons@UConn Master's Theses University of Connecticut Graduate School 5-9-2014 Target Species Selection for Dynamic Adaptive Chemistry Simulations Nicholas Curtis nicholas.curtis@uconn.edu
More informationA first investigation on using a species reaction mechanism for flame propagation and soot emissions in CFD of SI engines
A first investigation on using a 1000+ species reaction mechanism for flame propagation and soot emissions in CFD of SI engines F.A. Tap *, D. Goryntsev, C. Meijer, A. Starikov Dacolt International BV
More informationSystematic Reduction of Large Chemical Mechanisms
4th Joint Meeting of the U.S. Sections of the Combustion Institute, Philadelphia, PA, 2005 Systematic Reduction of Large Chemical Mechanisms Perrine Pepiot and Heinz Pitsch Department of Mechanical Engineering,
More informationVignolese 905/A, Modena, Italy
1 2 3 4 5 6 7 8 9 A study of direct and Krylov iterative sparse solver techniques to approach linear scaling of the integration of Chemical Kinetics with detailed combustion mechanisms Federico Perini
More informationFaster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs
Faster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs Christopher P. Stone, Ph.D. Computational Science and Engineering, LLC Kyle Niemeyer, Ph.D. Oregon State University 2 Outline
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationPractical Combustion Kinetics with CUDA
Funded by: U.S. Department of Energy Vehicle Technologies Program Program Manager: Gurpreet Singh & Leo Breton Practical Combustion Kinetics with CUDA GPU Technology Conference March 20, 2015 Russell Whitesides
More informationStrategies for Using Detailed Kinetics in Engine Simulations
Strategies for Using Detailed Kinetics in Engine Simulations Ellen Meeks ERC Symposium: Fuels for Future IC Engines June 6-7, 2007 Madison, Wisconsin Outline Role of simulation in design Importance of
More informationIMPLEMENTATION OF REDUCED MECHANISM IN COMPLEX CHEMICALLY REACTING FLOWS JATHAVEDA MAKTAL. Submitted in partial fulfillment of the requirements
IMPLEMENTATION OF REDUCED MECHANISM IN COMPLEX CHEMICALLY REACTING FLOWS by JATHAVEDA MAKTAL Submitted in partial fulfillment of the requirements For the degree of Master of Science in Aerospace Engineering
More informationThe Combination of Detailed Kinetics and CFD in Automotive Applications
The Combination of Detailed Kinetics and CFD in Automotive Applications J. M. Deur and S. Jonnavithula Analysis and Design Application Co., Ltd. Melville, New York E. Meeks Reaction Design San Diego, California
More informationProgress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling
Progress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling Michael McCourt 1,2,Lois Curfman McInnes 1 Hong Zhang 1,Ben Dudson 3,Sean Farley 1,4 Tom Rognlien 5, Maxim Umansky 5 Argonne National
More informationDevelopment of Reduced Chemistry for Efficient Combustion Simulations
University of Connecticut DigitalCommons@UConn Master's Theses University of Connecticut Graduate School 4-13-2016 Development of Reduced Chemistry for Efficient Combustion Simulations Yufeng Liu yufengliu812@gmail.com
More informationarxiv: v1 [physics.chem-ph] 5 Jul 2016
arxiv:1607.05079v1 [physics.chem-ph] 5 Jul 2016 Skeletal mechanism generation for surrogate fuels using directed relation graph with error propagation and sensitivity analysis Abstract Kyle E. Niemeyer
More informationMechanism reduction for multicomponent surrogates: A case study using toluene reference fuels
Mechanism reduction for multicomponent surrogates: A case study using toluene reference fuels Niemeyer, K. E., & Sung, C. J. (2014). Mechanism reduction for multicomponent surrogates: a case study using
More informationCoupling SUNDIALS with OpenFOAM 6 th OpenFOAM Workshop. E. David Huckaby June 14, 2011
Coupling SUNDIALS with OpenFOAM 6 th OpenFOAM Worshop E. David Hucaby June 14, 2011 EDH-OFW2011 Outline Motivation Mathematics Basic Design Eamples Food Web Steady & Unsteady Pitz-Daily 2 Motivation Why
More informationWILLKOMMEN WELCOME VÄLKOMMEN BENVIDO BIENVENIDO VELKOMMEN DOBRO DOSLI KARIBU WELKOM BENVENUTO SELAMAT DATANG BIENVENUE CROESO SOO DHAWOW NAMASTE
WELCOME VÄLKOMMEN BIENVENIDO DOBRO DOSLI WELKOM SELAMAT DATANG CROESO NAMASTE WILLKOMMEN BENVIDO VELKOMMEN KARIBU BENVENUTO BIENVENUE SOO DHAWOW 2016 Convergent Science. All Rights Reserved Combustion
More informationImplementation of an Alternate Chemistry Library into OpenFOAM TM
of an Alternate Chemistry Library into OpenFOAM TM Bernhard F.W. Gschaider 1 Markus Rehm 2 Peter Seifert 2 Bernd Meyer 2 1 ICE Strömungsforschung Leoben, Austria 2 Institute of Energy Process Enineering
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationCurrent progress in DARS model development for CFD
Current progress in DARS model development for CFD Harry Lehtiniemi STAR Global Conference 2012 Netherlands 20 March 2012 Application areas Automotive DICI SI PPC Fuel industry Conventional fuels Natural
More informationA comparison between two different Flamelet reduced order manifolds for non-premixed turbulent flames
8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 A comparison between two different Flamelet
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationEnhancing Scalability of Sparse Direct Methods
Journal of Physics: Conference Series 78 (007) 0 doi:0.088/7-6596/78//0 Enhancing Scalability of Sparse Direct Methods X.S. Li, J. Demmel, L. Grigori, M. Gu, J. Xia 5, S. Jardin 6, C. Sovinec 7, L.-Q.
More informationUQ in Reacting Flows
UQ in Reacting Flows Planetary Entry Simulations High-Temperature Reactive Flow During descent in the atmosphere vehicles experience extreme heating loads The design of the thermal protection system (TPS)
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationEffective Use of Storage/Retrieval-Based Chemistry Acceleration in CFD
2006 International Multidimensional Engine Modeling Users Group Meeting, Detroit, MI, 2 April 2006 Effective Use of Storage/Retrieval-Based Chemistry Acceleration in CFD Y. Wu 1, I. Veljkovic 2 and D.C.
More informationA wide range kinetic modelling study of laminar flame speeds of reference fuels and their mixtures
A wide range kinetic modelling study of laminar flame speeds of reference fuels and their mixtures A. Frassoldati, R. Grana, A. Cuoci, T. Faravelli, E. Ranzi Dipartimento di Chimica, Materiali e Ingegneria
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationTransported PDF Calculations of Combustion in Compression- Ignition Engines
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress Detroit, MI 15 April 2013 Transported PDF Calculations of Combustion in Compression- Ignition Engines V. Raj Mohan
More informationNUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL
NUMERICAL INVESTIGATION OF IGNITION DELAY TIMES IN A PSR OF GASOLINE FUEL F. S. Marra*, L. Acampora**, E. Martelli*** marra@irc.cnr.it *Istituto di Ricerche sulla Combustione CNR, Napoli, ITALY *Università
More informationComparison of a fully coupled and a decoupled solver for the simulation of Fluid Catalytic Cracking Background Fully coupled solver
Comparison of a fully coupled and a decoupled solver for the simulation of Fluid Catalytic Cracking Paper number 341c Edward Baudrez, 1 Juray De Wilde, 2 Geraldine Heynderickx 1, and Guy B. Marin 1 1.
More informationSimulation of soot formation and analysis of the "Advanced soot model" parameters in an internal combustion engine
Simulation of soot formation and analysis of the "Advanced soot model" parameters in an internal combustion engine Marko Ban, dipl. ing.* Power Engineering Department Faculty of Mechanical Engineering
More informationOpenSMOKE: NUMERICAL MODELING OF REACTING SYSTEMS WITH DETAILED KINETIC MECHANISMS
OpenSMOKE: NUMERICAL MODELING OF REACTING SYSTEMS WITH DETAILED KINETIC MECHANISMS A. Cuoci, A. Frassoldati, T. Faravelli, E. Ranzi alberto.cuoci@polimi.it Department of Chemistry, Materials, and Chemical
More informationIgnition delay-time study of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range at high pressure
25 th ICDERS August 2 7, 2015 Leeds, UK Ignition delay-time study of fuel-rich CH 4 /air and CH 4 /additive/air mixtures over a wide temperature range at high pressure Jürgen Herzler, Mustapha Fikri, Oliver
More informationCarbon Science and Technology
ASI RESEARCH ARTICLE Carbon Science and Technology Received:10/03/2016, Accepted:15/04/2016 ------------------------------------------------------------------------------------------------------------------------------
More informationSensitivity Analysis of Differential-Algebraic Equations and Partial Differential Equations
Sensitivity Analysis of Differential-Algebraic Equations and Partial Differential Equations Linda Petzold Department of Computer Science University of California Santa Barbara, California 93106 Shengtai
More informationComparison of ODE Solver for Chemical Kinetics and Reactive CFD Applications
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/269248243 Comparison of ODE Solver for Chemical Kinetics and Reactive CFD Applications CONFERENCE
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationOverview of Reacting Flow
Overview of Reacting Flow Outline Various Applications Overview of available reacting flow models Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Chemical Process Industry
More informationCombustion and Flame
Combustion and Flame 158 (2011) 2113 2127 Contents lists available at ScienceDirect Combustion and Flame journal homepage: www.elsevier.com/locate/combustflame Combined dimension reduction and tabulation
More informationOverview of Turbulent Reacting Flows
Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding
More informationNew sequential combustion technologies for heavy-duty gas turbines
New sequential combustion technologies for heavy-duty gas turbines Conference on Combustion in Switzerland 07.09.2017 ETH Zurich Nicolas Noiray, Oliver Schulz CAPS Lab D-MAVT ETH Nicolas Noiray 07/09/17
More informationRecent advances in sparse linear solver technology for semiconductor device simulation matrices
Recent advances in sparse linear solver technology for semiconductor device simulation matrices (Invited Paper) Olaf Schenk and Michael Hagemann Department of Computer Science University of Basel Basel,
More informationChemical Kinetic Reaction Mechanisms
Chemical Kinetic Reaction Mechanisms Charles Westbrook Lawrence Livermore National Laboratory May 31, 2008 J. Warnatz Memorial Colloquium How do we measure a career? In scientific careers, we look for
More informationAdvantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling
Advantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling Sudeep Bosu TATA Consultancy Services Baskaran Sundaram TATA Consultancy Services, 185 Lloyds Road,Chennai-600086,INDIA
More informationFuel 105 (2013) Contents lists available at SciVerse ScienceDirect. Fuel. journal homepage:
Fuel 105 (2013) 636 644 Contents lists available at SciVerse ScienceDirect Fuel journal homepage: www.elsevier.com/locate/fuel Simulations of a turbulent non-premixed flame using combined dimension reduction
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationThe influence of C ϕ is examined by performing calculations with the values C ϕ =1.2, 1.5, 2.0 and 3.0 for different chemistry mechanisms.
The Influence of Chemical Mechanisms on PDF Calculations of Nonpremixed Piloted Jet Flames Renfeng Cao and Stephen B. Pope Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca,
More informationRobust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations
Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations Rohit Gupta, Martin van Gijzen, Kees Vuik GPU Technology Conference 2012, San Jose CA. GPU Technology Conference 2012,
More informationReduced Chemical Kinetic Mechanisms for JP-8 Combustion
AIAA 2002-0336 Reduced Chemical Kinetic Mechanisms for JP-8 Combustion C. J. Montgomery Reaction Engineering International Salt Lake City, UT S. M. Cannon CFD Research Corporation Huntsville, AL M. A.
More informationModeling and Simulation of Plasma-Assisted Ignition and Combustion
Modeling and Simulation of Plasma-Assisted Ignition and Combustion Vigor Yang and Sharath Nagaraja Georgia Institute of Technology Atlanta, GA AFOSR MURI Fundamental Mechanisms, Predictive Modeling, and
More informationEffects of Hydrogen Addition on NOx Emissions in Hydrogen-Assisted Diesel Combustion
29 International Multidimensional Engine Modeling Users Group Meeting Detroit, MI, 19 April 29 Abstract Effects of Hydrogen Addition on NOx Emissions in Hydrogen-Assisted Diesel Combustion H. Zhang, G.K.
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationDevelopment of chemical reaction kinetics of VOC ozonation
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (2017) 000 000 www.elsevier.com/locate/procedia International High- Performance Built Environment Conference A Sustainable
More informationNumerical Modeling of Laminar, Reactive Flows with Detailed Kinetic Mechanisms
Department of Chemistry, Materials, and Chemical Engineering G. Natta Politecnico di Milano (Italy) A. Cuoci, A. Frassoldati, T. Faravelli and E. Ranzi Numerical Modeling of Laminar, Reactive Flows with
More informationA Systematic Investigation of In Situ Adaptive. Tabulation for Combustion Chemistry. Liuyan Lu, Stephen B. Pope
A Systematic Investigation of In Situ Adaptive Tabulation for Combustion Chemistry Liuyan Lu, Stephen B. Pope FDA 7- October 7 Abstract A systematic up-to-date investigation of different implementations
More informationImprovement in Computational Efficiency for HCCI Engine Modeling by Using Reduced Mechanisms and Parallel Computing
Improvement in Computational Efficiency for HCCI Engine Modeling by Using Reduced Mechanisms and arallel Computing Amr Ali *, Giulio Cazzoli, Song-Charng Kong, Rolf D. Reitz Engine Research Center, University
More informationNumerical modeling of auto-ignition of isolated fuel droplets in microgravity
Department of Chemistry, Materials, and Chemical Engineering Politecnico di Milano (Italy) A. Cuoci, A. Frassoldati, T. Faravelli and E. Ranzi Numerical modeling of auto-ignition of isolated fuel droplets
More informationP. Pepiot-Desjardins a & H. Pitsch a a Department of Mechanical Engineering, Stanford University, CA, Available online: 15 Nov 2008
This article was downloaded by: [Cornell University] On: 21 February 2012, At: 06:42 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationVALIDATION OF THE ASVDADD CONSTRAINT SELECTION ALGORITHM FOR EFFECTIVE RCCE MODELING OF NATURAL GAS IGNITION IN AIR
Proceedings of the ASME 2016 International Mechanical Engineering Congress & Exposition IMECE2016 November 11-17, 2016, Phoenix, Arizona, USA IMECE2016-65323 VALIDATION OF THE ASVDADD CONSTRAINT SELECTION
More informationModeling for Control of HCCI Engines
Modeling for Control of HCCI Engines Gregory M. Shaver J.Christian Gerdes Parag Jain Design Division Stanford, California 9435-421 Email: shaver@stanford.edu Design Division Stanford, California 9435-421
More informationModeling ion and electron profiles in methane-oxygen counterflow diffusion flames
Abstract 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 Modeling ion and electron profiles
More informationAn Efficient Graph Sparsification Approach to Scalable Harmonic Balance (HB) Analysis of Strongly Nonlinear RF Circuits
Design Automation Group An Efficient Graph Sparsification Approach to Scalable Harmonic Balance (HB) Analysis of Strongly Nonlinear RF Circuits Authors : Lengfei Han (Speaker) Xueqian Zhao Dr. Zhuo Feng
More informationA robust multilevel approximate inverse preconditioner for symmetric positive definite matrices
DICEA DEPARTMENT OF CIVIL, ENVIRONMENTAL AND ARCHITECTURAL ENGINEERING PhD SCHOOL CIVIL AND ENVIRONMENTAL ENGINEERING SCIENCES XXX CYCLE A robust multilevel approximate inverse preconditioner for symmetric
More informationA REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL ENTROPY INEQUALITY
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Special Issue/2018, pp. 285 290 A REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL
More informationAdvanced Multi-Physics Modeling & Simulation Efforts for Fast Reactors
Advanced Multi-Physics Modeling & Simulation Efforts for Fast Reactors J.W. Thomas Argonne National Laboratory IAEA Technical Meeting on Priorities in Modeling & Simulation in Fast Reactor Systems April
More informationAssignment on iterative solution methods and preconditioning
Division of Scientific Computing, Department of Information Technology, Uppsala University Numerical Linear Algebra October-November, 2018 Assignment on iterative solution methods and preconditioning 1.
More informationPreconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners
Preconditioning Techniques for Large Linear Systems Part III: General-Purpose Algebraic Preconditioners Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, Georgia, USA
More informationA New Parallel Numerically improved KIVA3V Program for Diesel Combustion Computations.
A New Parallel Numerically improved KIVA3V Program for Diesel Combustion Computations. P.Belardini, C.Bertoli, P.DeMarino, G.Avolio. Istituto Motori National Research Council of Italy. Via Marconi 8 8125
More informationLecture 17: Iterative Methods and Sparse Linear Algebra
Lecture 17: Iterative Methods and Sparse Linear Algebra David Bindel 25 Mar 2014 Logistics HW 3 extended to Wednesday after break HW 4 should come out Monday after break Still need project description
More informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationFAS and Solver Performance
FAS and Solver Performance Matthew Knepley Mathematics and Computer Science Division Argonne National Laboratory Fall AMS Central Section Meeting Chicago, IL Oct 05 06, 2007 M. Knepley (ANL) FAS AMS 07
More informationApproximation of chemical reaction rates in turbulent combustion simulation
Approximation of chemical reaction rates in turbulent combustion simulation Lars Frank Große and Franz Joos * Helmut-Schmidt-University University of the Federal Armed Forces Hamburg - Laboratory of Turbo
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationCombustion Chemistry of a New Biofuel: Butanol
Combustion Chemistry of a New Biofuel: Butanol William H. Green, D.F. Davidson, F. Egolfopoulos, N. Hansen, M. Harper, R.K. Hanson, S. Klippenstein, C.K. Law, C.J. Sung, D.R. Truhlar, & H. Wang Assessing
More informationX. He and M. Neytcheva. Preconditioning the incompressible Navier-Stokes equations with variable viscosity. J. Comput. Math., in press, 2012.
List of papers This thesis is based on the following papers, which are referred to in the text by their Roman numerals. I II III IV V VI M. Neytcheva, M. Do-Quang and X. He. Element-by-element Schur complement
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More informationCombustion and Emission Modeling in CONVERGE with LOGE models
Combustion and Emission Modeling in CONVERGE with LOGE models Corinna Netzer, Harry Lehtiniemi and Fabian Mauss 2015 CONVERGE USER CONFERENCE RICHARD CHILDRESS RACING, WELCOME, NC Outline Objective LOGE
More informationExercises in Combustion Technology
Exercises in Combustion Technology Exercise 4: Turbulent Premixed Flames Turbulent Flow: Task 1: Estimation of Turbulence Quantities Borghi-Peters diagram for premixed combustion Task 2: Derivation of
More informationThe Chemical Kinetics Time Step a detailed lecture. Andrew Conley ACOM Division
The Chemical Kinetics Time Step a detailed lecture Andrew Conley ACOM Division Simulation Time Step Deep convection Shallow convection Stratiform tend (sedimentation, detrain, cloud fraction, microphysics)
More informationFINE-GRAINED PARALLEL INCOMPLETE LU FACTORIZATION
FINE-GRAINED PARALLEL INCOMPLETE LU FACTORIZATION EDMOND CHOW AND AFTAB PATEL Abstract. This paper presents a new fine-grained parallel algorithm for computing an incomplete LU factorization. All nonzeros
More informationComputationally-Efficient Parallel Implementation of Combustion Chemistry in LES/PDF Computations
Fall Technical Meeting of the Eastern States Section of the Combustion Institute Hosted by the University of Connecticut, Storrs, CT October 9-12, 2011 Computationally-Efficient Parallel Implementation
More informationA Study of Grid Resolution and SGS Models for LES under Non-reacting Spray Conditions
ILASS Americas, 25 th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013 A Study of Grid Resolution and SGS Models for LES under Non-reacting Spray Conditions Q. Xue 1*,
More informationParallel on-the-fly adaptive kinetics in direct numerical simulation of turbulent premixed flame
Available online at www.sciencedirect.com Proceedings of the Combustion Institute 36 (017) 05 03 www.elsevier.com/locate/proci Parallel on-the-fly adaptive kinetics in direct numerical simulation of turbulent
More informationNUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER
NUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER Naveen Kumar D 1*, Pradeep R 2 and Bhaktavatsala H R 3 1 Assistant Professor Department of Mechanical Engineering, M S Engineering College,
More informationParallelization of Multilevel Preconditioners Constructed from Inverse-Based ILUs on Shared-Memory Multiprocessors
Parallelization of Multilevel Preconditioners Constructed from Inverse-Based ILUs on Shared-Memory Multiprocessors J.I. Aliaga 1 M. Bollhöfer 2 A.F. Martín 1 E.S. Quintana-Ortí 1 1 Deparment of Computer
More informationStructure, Extinction, and Ignition of Non-Premixed Flames in the Counterflow Configuration
Structure, Extinction, and Ignition of Non-Premixed Flames in the Counterflow Configuration Ryan Gehmlich STAR Global Conference 2013 Orlanda, Florida March 18-20 1 Outline Background Developing Reaction
More informationTOPICAL PROBLEMS OF FLUID MECHANICS 97
TOPICAL PROBLEMS OF FLUID MECHANICS 97 DOI: http://dx.doi.org/10.14311/tpfm.2016.014 DESIGN OF COMBUSTION CHAMBER FOR FLAME FRONT VISUALISATION AND FIRST NUMERICAL SIMULATION J. Kouba, J. Novotný, J. Nožička
More informationDEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
DRAFT Proceedings of ASME IMECE: International Mechanical Engineering Conference & Exposition Chicago, Illinois Nov. 5-10, 2006 IMECE2006-14867 DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
More information