New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

Total Page:16

File Type:pdf, Size:1020Kb

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics SCIENCE CHINA Physics, Mechanics & Astronomy • Research Paper • February 2010 Vol. 53 No.2: 237−243 Special Topics: Computational Explosion Mechanics doi: 10.1007/s11433-009-0266-z New numerical algorithms in SUPER CE/SE and their applications in explosion mechanics WANG Gang1,2, WANG JingTao1 & LIU KaiXin1,3* 1 LTCS and College of Engineering, Peking University, Beijing 100871, China; 2 LHD, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China; 3 Center for Applied Physics and Technology, Peking University, Beijing 100871, China Received June 28, 2009; accepted November 11, 2009 The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched al- gorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algo- rithms presented in this paper may be applied to the numerical research of explosion mechanics. new algorithm, SUPER CE/SE, numerical simulation, explosion mechanics The research of explosion mechanics, usually involving jing Institute of Applied Physics and Computational many disciplines such as compressible flows, incompressi- Mathematics) and MMIC [2] by BIT (Beijing Institute of ble flows, elastic-plastic flows, chemically reacting flows Technology). The present software packages of computa- and two-phase flows, contributes significantly to the field of tional explosion mechanics mainly adopt well-developed industry, aeronautics, and military engineering. Numerical numerical algorithms and models, for example, VOF [3] simulation of explosion mechanics needs not only accurate (Volume of Fluid) method for interface capturing and tra- numerical schemes and interface capturing methods, but jectory method [4] for particle tracing. also a perfect combination of the selected numerical algo- This paper has studied the new numerical algorithms in rithms and the physical and chemical models. SUPER CE/SE and their applications in explosion mechan- The imported commercial software, such as LS-DYNA ics. SUPER CE/SE uses an improved CE/SE (space-time and AUTODYN, currently plays a dominant role in the ap- Conservation Element and Solution Element) method [5,6] plication of explosion mechanics in China. Further devel- to discrete governing equations, a local hybrid particle level opment based on the imported commercial software is quite set method to capture the interface, and a two-fluid model to difficult, because the kernel modules of the software are treat two-phase flows. SUPER CE/SE is also able to adopt always held as commercial secrets. The software of compu- common chemical reaction models, such as two-step, de- tational explosion mechanics in China has undergone a tailed and Sichel’s two-step chemical reaction models. With step-by-step development. The well-established software the above algorithms and models, problems of shock wave packages include METH [1] developed by BIAPCM (Bei- reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in gas-droplet system are investigated. Comparisons and dis- *Corresponding author (email: [email protected]) © Science China Press and Springer-Verlag Berlin Heidelberg 2010 phys.scichina.com www.springerlink.com 238 WANG Gang, et al. Sci China Phys Mech Astron February (2010) Vol. 53 No. 2 cussions on the numerical results reveal the applicability of advantages such as simple mesh structure, clear physics the above algorithms in computational explosion mechan- concept, easy to extend to three-dimensional cases and ics. convenient high-accuracy Taylor expansions. The physical variables Q, E and F at every mesh point (t, x, y)∈SE(j, k, n) are approximated by their second-order Taylor expansions 1 Numerical algorithms and models Q*, E* and F* as * Computational explosion mechanics refers to many areas of M (,δδδxyt ,)P′ study including computational mechanics, fluid mechanics, =+MMPP′′()xytδδδxyt + () M P ′ + () M P ′ solid mechanics, physics and chemistry. In this section, we 111222 mainly introduce the basic algorithms in SUPER CE/SE +++()()()()()()MMMxx PPP′′′δ xytyy δδtt including the improved CE/SE scheme, the local hybrid 222 2 (1) particle level set method, the common chemical reaction +++()MMMxy PPP′′′δδxy ()yt δδ yt ()xt δδ xt , models and the two-fluid model. The models of compressi- ble flow, incompressible flow and elastic-plastic flow can where M=Q, E, F, δ xxx= − P′ , δ yyy=−P′ , δ ttt= − P′ . be found in ref. [7,8]. A second-order accuracy CE/SE scheme can be constructed by eq. (1). 1.1 Improved CE/SE method 1.2 Local hybrid particle level set method The CE/SE method [9] originally proposed by Chang in 1995 is a new numerical framework for solving hyperbolic How to accurately trace the interface of materials is a chal- conservation laws. The CE/SE method substantially differs lenge confronted in the Eulerian method for multi-material from the existing well-established CFD methods in that it simulations. The Level Set method [12], originally proposed contains many excellent features such as a unified treatment by Osher and Sethian in 1988, is an effective Eulerian algo- of space and time, satisfying both local and global flux rithm for interface capturing. The level set method essen- conservation in space and time and simple treatments for tially refers to a level set function whose zero isoline ex- boundary conditions. As a result of its simplicity, generality actly determines the material interface is introduced to de- and accuracy, the CE/SE method has been successfully ap- scribe the development of the material interface. The level plied to scientific exploration and engineering. However set method can deal with very complicated topologic researchers have found defects in the original CE/SE changes of interfaces, because it is not necessary to recon- schemes, such as the complication of mesh structure, the struct interface fronts just like some “wave front tracing” absence of high-order accuracy schemes and the difficulty methods. to extend to three-dimensional situation. To solve these Adalsteinsson [13] proposed a fast local algorithm for the problems, Zhang et al. [10], Zhang et al. [11] and Liu et al. level set method which confines the solving region to a [5,6] have made some progress in this area. narrow band surrounding the interface. This algorithm has The improved CE/SE method used in SUPER CE/SE the time complexity of O(Nw), where N is the number of constructs SEs and CEs shown in Figure 1 on rectangular grid points and w is the width of the narrow band. Applying meshes and is extended to second-order accuracy. This con- the idea of fast local algorithm, Peng [14] proposed a fast struction of SE and CE on rectangular meshes has many level set method based on PDE (Partial Difference Equa- Figure 1 Mesh construction of the improved CE/SE method. (a) Solution element SE(P′); (b) conservation element CE(P′). WANG Gang, et al. Sci China Phys Mech Astron February (2010) Vol. 53 No. 2 239 tions), which has merits such as easy to be implemented and 1.4 Two-fluid model time complexity of O(N). In this paper, we assume the gas-droplet system has the fol- Enright [15] introduced non-mass Lagrangian particles to lowing characteristics: the flow is unsteady; the initial tem- the level set method to solve the problem of mass losing. perature of the gas and the droplets is the same; the radii of Based on this idea, he proposed a hybrid particle Level Set droplets are uniform; droplets are homogenously distributed method in 2002. initially; the phase of droplets is considered a continuous medium; the shape of droplets always remains spherical 1.3 Common chemical reaction models even in the process of separation, and evaporation; the tem- perature distribution in every droplet is uniform; the total The two-step reaction model simplifies a complicated volume of droplets can be ignored compared with the vol- chemical reaction into an induction reaction and an exo- ume of gas; the interactions between droplets are ignored; thermic reaction. The rates of a (progress parameter for in- when the fuel reaches the gaseous state, chemical reactions duction reaction) and β (progress parameter for exothermic occur and accomplish immediately; the chemical energy is reaction), ω and ω , are given as below [16]: α β absorbed only by gas, which is considered ideal. With these assumptions, the governing equations for gas ⎧ ⎛⎞Eα ⎪ωραα=−k exp⎜⎟ − , phase can be expressed as ⎪ ⎝⎠RT ⎪ 0 (α > 0), ∂∂∂QEF ⎧ S ⎪ ⎪ + +=, (4) ⎪ ⎡ E ∂∂∂txy ⎨ ⎪ +− 22 ⎛⎞β (2) ωωωββββ=+=−⎪ kp⎢ βexp⎜⎟ − ⎪ ⎢ RT ⎪ ⎨⎣⎝⎠ where ⎪ ⎪ EQ+ ⎤ Q = (ρ, ρu, ρv, E), ⎪ 2 ⎛⎞β 2 ⎪ −−(1βα ) exp⎜⎟ −⎥ (0), ≤ E = ( ρu, ρu +p, ρuv, (E+p)u), ⎪ RT 2 ⎩⎪ ⎩ ⎝⎠⎦⎥ F = (ρv, ρuv, ρv +p, (E+p)v), S = (Ip, −Fx+upIp, −Fy+vpIp, −(upFx+ vpFy) where ρ is the mass density, p the pressure, T the temperature, 2 2 +((up +up )/2+qr)Ip). R the gas constant, Q the heat release parameter, kα and kβ The governing equations for droplet phase can be expressed the constants of reaction rates, and Eα and Eβ the activation as energies. The two-step reaction model does not need large ∂∂∂QEF computing resource because it only takes the heat release pp++=p S , (5) parameter Q into account. However, the two-step reaction ∂∂∂txyp model can describe the essential characteristics of detona- tions. where The detailed chemical reaction model is extensively used Qp = (ρp, ρpup, ρpvp, N), E = (ρ u , ρ u 2, ρ u v Nu ), to describe the transformation of reactants into products at p p p p p p p p, p F = (ρ v , ρ u v ρ v 2, Nv ), the molecular level through a large number of elementary p p p p p p, p p p S = −(I , −F +u I , −F +v I , 0).
Recommended publications
  • What Technology Wants / Kevin Kelly
    WHAT TECHNOLOGY WANTS ALSO BY KEVIN KELLY Out of Control: The New Biology of Machines, Social Systems, and the Economic World New Rules for the New Economy: 10 Radical Strategies for a Connected World Asia Grace WHAT TECHNOLOGY WANTS KEVIN KELLY VIKING VIKING Published by the Penguin Group Penguin Group (USA) Inc., 375 Hudson Street, New York, New York 10014, U.S.A. Penguin Group (Canada), 90 Eglinton Avenue East, Suite 700, Toronto, Ontario, Canada M4P 2Y3 (a division of Pearson Penguin Canada Inc.) Penguin Books Ltd, 80 Strand, London WC2R 0RL, England Penguin Ireland, 25 St. Stephen's Green, Dublin 2, Ireland (a division of Penguin Books Ltd) Penguin Books Australia Ltd, 250 Camberwell Road, Camberwell, Victoria 3124, Australia (a division of Pearson Australia Group Pty Ltd) Penguin Books India Pvt Ltd, 11 Community Centre, Panchsheel Park, New Delhi - 110 017, India Penguin Group (NZ), 67 Apollo Drive, Rosedale, North Shore 0632, New Zealand (a division of Pearson New Zealand Ltd) Penguin Books (South Africa) (Pty) Ltd, 24 Sturdee Avenue, Rosebank, Johannesburg 2196, South Africa Penguin Books Ltd, Registered Offices: 80 Strand, London WC2R 0RL, England First published in 2010 by Viking Penguin, a member of Penguin Group (USA) Inc. 13579 10 8642 Copyright © Kevin Kelly, 2010 All rights reserved LIBRARY OF CONGRESS CATALOGING IN PUBLICATION DATA Kelly, Kevin, 1952- What technology wants / Kevin Kelly. p. cm. Includes bibliographical references and index. ISBN 978-0-670-02215-1 1. Technology'—Social aspects. 2. Technology and civilization. I. Title. T14.5.K45 2010 303.48'3—dc22 2010013915 Printed in the United States of America Without limiting the rights under copyright reserved above, no part of this publication may be reproduced, stored in or introduced into a retrieval system, or transmitted, in any form or by any means (electronic, mechanical, photocopying, recording or otherwise), without the prior written permission of both the copyright owner and the above publisher of this book.
    [Show full text]
  • Experiences Tuning a J90 for Engineering
    ExperiencesExperiences TuningTuning aa J90J90 forfor EngineeringEngineering ApplicationsApplications Sharan Kalwani [email protected] CUG 1997, San Jose MotivationMotivation ● Why do sites choose a J90? – Cost is usually the primary consideration – Computational needs satisfied by J90 – Starting point for future growth – Resources are also a major factor MotivationMotivation (contd...)(contd...) ● Field experience revealed: – Several sites out there (250+) – Site profiles were all different , eg: ● Human resources often limited ● Expertise scarce ● Learning curve was steep ● Balancing act of system versus “real” work BenefitsBenefits ● Share useful tips and techniques ● Gain insight into other possible scenarios ● Learn what works ● Learn what may NOT work CaveatCaveat EmptorEmptor ● Target Audience – Engineering Applications ● Computational Fluid Dynamics (CFD) – STAR-CD, KIVA-2, CHAD ● Metal Deformation Codes – (LS-DYNA) ● Finite Element Analysis (FEA) – NASTRAN, etc... TypicalTypical J90J90 EnvironmentEnvironment ● Usually small systems: – 4 to 8 CPUs – 64 Mwords – Single IOS – SCSI disks (sometimes also IPI disks) – small disk farm – no large online tapes (4 or 8mm only) J90J90 EnvironmentEnvironment (contd)(contd) ● Staff at sites: – single, but multi-tasking person :-) – has a real job in addition to Cray work – usually sophisticated end user, but – not a “pure” systems type TuningTuning RegionsRegions ● UNICOS Kernel ● Memory and ● Process Scheduling ● Disk Drive Strategy Tuning:Tuning: J90J90 KernelKernel ArenaArena ● UNICOS Kernel
    [Show full text]
  • Mechanical Engineering Magazine Will (Subject Line "Letters and Comments")
    Mechanical THE MAGAZINE OF ASME ENGINEERING No. 07 139 Technology that moves the world SMALL BUT MIGHTY Robots are boosting warehouse productivity. FASTEST WOMAN ON TWO WHEELS PAGE 18 FLEXIBLE POWER SUIT PAGE 38 BRAIN-CONTROLLED PROSTHETICS PAGE 44 ASME.ORG JULY 2017 GREAT NEWS! Thanks to our partnership with one of the leading engineering professional liability insurance carriers, Beazley Insurance &RPSDQ\ΖQF\RXKDYHDFFHVVWRDQH[SHUWUHVRXUFHIRU\RXUȴUPȇV professional liability needs. And because of this partnership, if you are an ASME member you are eligible for a 10% premium discount* exclusively through our program. ASME PROFESSIONAL LIABILITY PROGRAM DESIGNED SPECIFICALLY FOR ASME 7KH3URIHVVLRQDO/LDELOLW\3URJUDPLVR΍HUHGWKURXJKDSURYHQ industry leader in the engineering professional liability marketplace: ASME MEMBER EXCLUSIVE COVERAGE BENEFITS: Beazley is one of the top engineering insurance carriers that understands the unique risks engineers face. That’s why the ȏ 10% premium discount* through 3URIHVVLRQDO/LDELOLW\3URJUDPLQFOXGHVȵH[LEOHSROLF\OLPLWV our program DQGGHGXFWLEOHRSWLRQVΖWFRYHUVȴUPVRIDOOVL]HVDVZHOODV ȏ Consent to settle enhancement VHOIHPSOR\HGLQGLYLGXDOV ȏ Early resolution deductible credit Call today at 1.800.640.7637 or visit www.asmeinsurance.com/pl *Based on state filed minimum premium requirements. Administered by Mercer Health & Benefits Administration LLC In CA d/b/a Mercer Health & Benefits Insurance Services LLCt CA Insurance License #0G39709 t AR Insurance License #100102691 This plan is underwritten by Beazley Insurance Company, Inc. t CA Insurance License #0G55497 t WA Insurance License #1298 79128 (7/17) Copyright 2017 Mercer LLC. All rights reserved. LOG ON ASME.ORG MECHANICAL ENGINEERING 01 | JULY 2017 | P. For these articles and other content, visit asme.org. Hybrid Cars Compete for Speed The Formula Hybrid competition is entering its second decade, giving developers and students an opportunity to race hybrid and electric vehicles for speed and endurance.
    [Show full text]
  • Spray Combustion Cross-Cut Engine Research
    Spray Combustion Cross-Cut Engine Research Lyle M. Pickett, Scott A. Skeen Sandia National Laboratories FY 2017 DOE Vehicle Technologies Office Annual Merit Review Project ACS005, 11:30 AM, Wed. 7 June 2017 Sponsor: DOE Vehicle Technologies Office Program Managers: Gurpreet Singh and Leo Breton This presentation does not contain any proprietary, confidential, or otherwise restricted information. Overview Timeline Barriers ● Project provides fundamental ● Engine efficiency and emissions research that supports DOE/ ● Understanding direct-injection industry advanced engine sprays development projects. ● CFD model improvement for ● Project directions and engine design/optimization continuation are evaluated annually. Partners ● 15 Industry partners in MOU: Budget Advanced Engine Combustion ● Project funded by DOE/VTO: ● Engine Combustion Network FY17 - $950K – >20 experimental + >20 modeling – >100 participants attend ECN5 ● Project lead: Sandia – Lyle Pickett (PI), Scott Skeen 2 Engine efficiency gains require fuel (DI spray) delivery optimization ● Barriers for high-efficiency gasoline – Particulate emissions 8-hole, gasoline 80° total angle – Engine knock ~15mm – Slow burn rate or partial burn – Heat release control when using compression ignition – Lack of predictive CFD tools ● Barriers for high-efficiency diesel – Particulate emissions 573 K, 3.5 kg/m3 573 K, 3.5 kg/m3 – Heat release rate and phasing – Lack of predictive CFD, particularly for short and multiple injections θ=22.5° I/I0 0.8-1.0 High-speed microscopy at nozzle exit 3 Project
    [Show full text]
  • Improved Solvers for Advanced Engine Combustion Simulation
    Lawrence Livermore National Laboratory Improved Solvers for Advanced Engine Combustion Simulation M. J. McNenly (PI), S. M. Aceves, D. L. Flowers, N. J. Killingsworth, G. M. Oxberry, G. Petitpas and R. A. Whitesides Project ID # ACE076 2014 DOE Vehicle Technologies Program Annual Merit Review and Peer Evaluation Meeting June 17, 2014 - Washington, DC This presentation does not contain any proprietary, confidential or otherwise restricted information This work performed under the auspices of the U.S. Department of Energy by LLNL-PRES-653560 Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344 Overview Timeline Barriers . Ongoing project with yearly . Lack of fundamental knowledge of direction from DOE advanced engine combustion regimes . Lack of modeling capability for combustion and emission control . Lack of effective engine controls Budget Partners . FY12 funding: $340K . Ford, GM, Bosch, Volvo, Cummins, . FY13 funding: $340K Convergent Sciences & NVIDIA . FY14 funding: $475K . Argonne NL, Sandia NL, Oak Ridge NL . UC Berkeley, U. Wisc., U. Mich., LSU, Indiana U. & RWTH Aachen . FACE working group, AEC MOU, SAE, Combustion Inst., GPU Tech. Conf. & NSF HPC software planning Lawrence Livermore National Laboratory McNenly, et al. LLNL-PRES-653560 2 Relevance – Enhanced understanding of HECC requires expensive models that fully couple detailed kinetics with CFD Objective We want to use… Create faster and more accurate Detailed chemistry combustion solvers. Ex. Biodiesel component C20H42 (LLNL) Accelerates R&D on three major 7.2K species challenges identified in the VT multi- 53K reaction steps year program plan: A. Lack of fundamental knowledge in highly resolved 3D simulations of advanced engine combustion regimes C. Lack of modeling capability for combustion and emission control D.
    [Show full text]
  • Stochastic Spray and Chemically Reacting Flow in LS-DYNA®
    9th European LS-DYNA Conference 2013 _________________________________________________________________________________ Stochastic Spray and Chemically Reacting Flow In LS-DYNA® Kyoung-Su Im, Zen-Chan Zhang, and Grant O. Cook, Jr. Livermore Software Technology Corp. Livermore, CA 94551 1 Introduction The injection of fuel sprays into an automotive engine and liquid jets into a high-speed flow stream is an important process in modern automotive gasoline and diesel engines, and propulsion in gas turbine and supersonic vehicles. In such applications, the combustion performance depends strongly on spray atomization, penetration, and the mixing process between the free stream air and the liquid fuel. As a result, the study of liquid spray in such areas has become an important research topic. The spray can be represented by a stochastic equation giving the rate of change of the distribution function, f = , which at a given time changes the positions, x, velocities, v, Ý equilibrium radius, r, temperature, f x,v,r,T,y,y ,T,t distortion from sphericity, y, and the time rate of change dy/dt[1,2]. This equation is then coupled with ensemble averaged gas phase equations of mass, momentum, and energy conservation. In addition, for many engineering sprays, drop breakup and collisions must be considered when the drop Webber number is larger than a critical value. Thus the essential dynamics of spray and its interaction with a gas is an extremely complicated physical phenomenon. Much attention has recently been devoted to the chemically reactive simulation models in the fields of subsonic, supersonic, and hypersonic combustion. However, the modeling of chemically reacting flow is a challenging task because of the simultaneous contribution of a wide range of time scales to the system dynamics.
    [Show full text]
  • Large Eddy Simulation of Dispersed Multiphase Flow
    Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Dissertations, Master's Theses and Master's Reports - Open Reports 2012 Large Eddy Simulation of dispersed multiphase flow Yejun Gong Michigan Technological University Follow this and additional works at: https://digitalcommons.mtu.edu/etds Part of the Applied Mathematics Commons, and the Mechanical Engineering Commons Copyright 2012 Yejun Gong Recommended Citation Gong, Yejun, "Large Eddy Simulation of dispersed multiphase flow", Dissertation, Michigan Technological University, 2012. https://doi.org/10.37099/mtu.dc.etds/729 Follow this and additional works at: https://digitalcommons.mtu.edu/etds Part of the Applied Mathematics Commons, and the Mechanical Engineering Commons LARGE EDDY SIMULATION OF DISPERSED MULTIPHASE FLOW By Gong, Yejun A DISSERTATION Submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY (Computational Science & Engineering) MICHIGAN TECHNOLOGICAL UNIVERSITY 2012 c 2012 Gong, Yejun This dissertation, "Large Eddy Simulation of Dispersed Multiphase Flow," is hereby ap- proved in partial fulfillment of the requirements for the Degree of DOCTOR OF PHILOS- OPHY IN COMPUTATIONAL SCIENCE & ENGINEERING. Department of Mathematical Science Signatures: Dissertation Advisor Dr. Tanner, Franz X Committee Member Dr. Feigl, Kathleen A Committee Member Dr. Yang, Song Lin Committee Member Dr. Kolkka, Robert W Department Chair Dr. Gockenbach, Mark S Date Dedication "Large Eddy Simulation of Dispersed Multiphase Flow" Ph.D. thesis by Yejun Gong, Dedicated to My parents and teachers and friends and in loving memory of my grandparents! . Contents List of Figures ..................................... xi List of Tables ...................................... xv Acknowledgments ...................................xvii ABBREVIATIONS ..................................xviii Abstract .........................................xxiii 1 Introduction ...................................
    [Show full text]
  • Notes on the Kiva-Ii Software and Chemically Reactive Fluid Mechanics
    NOTES ON THE KIVA-II SOFTWARE AND CHEMICALLY REACTIVE FLUID MECHANICS MICHAEL J. HOLST Numerical Mathematics Group Computing & Mathematics Research Division Lawrence Livermore National Laboratory Livermore, California NOTES ON THE KIVA-II SOFTWARE AND CHEMICALLY REACTIVE FLUID MECHANICS∗ Michael J. Holst Numerical Mathematics Group Computing & Mathematics Research Division Lawrence Livermore National Laboratory August 1, 1992 Abstract This report represents a set of working notes regarding the mechanics of chemically reactive fluids with sprays, and their numerical simulation with the KIVA-II software. KIVA-II is a large FORTRAN program developed at Los Alamos National Laboratory for internal combustion engine simulation. It is our hope is that these notes summarize some of the necessary background material in fluid mechanics and combustion, explain the numerical methods currently used in KIVA-II and similar combustion codes, and provide an outline of the overall structure of KIVA-II as a representative combustion program, in order to aid the researcher in the task of implementing KIVA-II or a similar combustion code on a massively parallel computer. The notes are organized into three parts as follows. In Part I, we give a brief intro- duction to continuum mechanics, to fluid mechanics, and to the mechanics of chemically reactive fluids with sprays. In Part II, we take a close look at the governing equations of KIVA-II, and discuss the methods employed in the numerical solution of these equations. We draw some conclusions and make some observations in Part III. ∗These notes were written while the author was visiting Lawrence Livermore National Laboratory during the summer of 1992, on loan from the Numerical Computing Group in the Department of Computer Science at the University of Illinois, Urbana-Champaign.
    [Show full text]
  • Download Article (PDF)
    International Conference on Intelligent Control and Computer Application (ICCA 2016) Grid generation research for intake-exhaust closing internal combustion engine using KIVA software Zhao Tianpena Liu Yongfengb*, Shi Yanc Beijing Research Centre of Sustainable Energy and Building Beijing Research Centre of Sustainable Energy and Building Beijing University of Civil Engineering and Architecture Beijing University of Civil Engineering and Architecture Beijing 100044, China Beijing 100044, China [email protected] [email protected], [email protected] Abstract—In order to simulate the combustion process of the of simple blocks. Therefore, the program can calculate very internal combustion engine, a new grid generation was carried. complex geometries. It can not only save the storage space, This paper developed the KIVA-3V software to establish a but also improves the calculation speed. In 1997, KIVA-3V gasoline engine combustion chamber model. Firstly, the which is KIVA-3's improved version was complete. KIVA-3V schematic of grid generation and block of combustion zone is increased oil film flow model, turbulent combustion model given out. The data file (IPREP) and other files are built up, and and soot model. the grid data of gasoline engine combustion chamber is obtained. Secondly, the output files are obtained by running the KIVA-3V. The development of simulation technology in China is Finally, the new three dimensional grids of the intake-exhaust relatively late. Two dimensional simulation is first applied. closing internal combustion engine are obtained. It gives a new The most representative is Professor Wang Rongsheng and way to build up simulation grid for internal engine. Professor Xie Maozhao of Dalian University of Technology.
    [Show full text]
  • The Projection Method Truegrid® Meshes Are Formed More Quickly and Easily Using the Projection Method
    Trademarks: are registered trademarks of Dassault Systems. CADDS™ is a trademark and Pro/ENGINEER® is a registered trademark of PTC. ANSYS®, AUTODYN®, and FLUENT® are registered trademarks of ANSYS, Inc. TASCflow® is a registered trademark of Advanced Scientific Computing, Ltd. NASTRAN® is a registered trademark of NASA. MSC NASTRAN® and PATRAN® are registered trademarks MSC.Software Corporation. Solid Edge® is a registered trademark Siemens Corporation. NX is a trademark or registered trademark of Siemens Product Lifecycle Management Software, Inc. LS-DYNA® is a registered trademark Livermore Software Technology Corporation. Windows® is a registered trademark of Microsoft Corporation. SolidWorks® is a registered trademark of Dassault Systems SolidWorks Corporation. SurfCam® CAD/CAM is a registered trademark of Surfware, Inc. Sun® and Solaris® are registered trademarks of Sun Microsystems, Inc. FIDAP™ is a trademark of Fluent, Inc. STAR-CD® is a registered trademark of CD-adapco. Unix® is a registered trademark of X/Open Company Limited. IBM® and AIX™ are trademarks or registered trademarks of IBM Corporation. DEC® is a registered trademark of Hewlett- Packard. CFX is a trademark of Sony Corporation. Intel® and Itanium® are registered trademarks of Intel Corporation. Red Hat® is a registered trademark of Red Hat, Inc. Apple® is a registered trademark of Apple, Inc. AMD™ and Opteron™ are trademarks of Advanced Micro Devices Inc. SUSE™ is a trade mark of NOVELL Inc. LINUX® is a registered trademark owned by Linus Torvalds. All other brand, product, service and feature names or trademarks are the property of their respective owners. The TrueGrid® Advantage TrueGrid® software is advanced meshing technology for FEA and CFD simulation codes.
    [Show full text]
  • NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector
    NASA Numerical and Experimental Evaluation of UTRC Low Emissions Injector Yolanda R. Hicks,1 Sarah A. Tedder,2 Robert C. Anderson,3 and Anthony C. Iannetti4 NASA Glenn Research Center, Cleveland, Ohio, 44135, USA Lance L. Smith5 and Zhongtao Dai6 United Technologies Research Center, East Hartford, CT, 06108, USA Computational and experimental analyses of a PICS—Pilot-In-Can-Swirler technology injector, developed by United Technologies Research Center (UTRC) are presented. NASA has defined technology targets for near term (called “N+1”, circa 2015), midterm (“N+2”, circa 2020) and far term (“N+3”, circa 2030) that specify realistic emissions and fuel efficiency goals for commercial aircraft. This injector has potential for application in an engine to meet the Pratt & Whitney N+3 supersonic cycle goals, or the subsonic N+2 engine cycle goals. Experimental methods were employed to investigate supersonic cruise points as well as select points of the subsonic cycle engine; cruise, approach, and idle with a slightly elevated inlet pressure. Experiments at NASA employed gas analysis and a suite of laser-based measurement techniques to characterize the combustor flow downstream from the PICS dump plane. Optical diagnostics employed for this work included Planar Laser-Induced Fluorescence of fuel for injector spray pattern and Spontaneous Raman Spectroscopy for relative species concentration of fuel and CO2. The work reported here used unheated (liquid) Jet-A fuel for all fuel circuits and cycle conditions. The initial tests performed by UTRC used vaporized Jet-A to simulate the expected supersonic cruise condition, which anticipated using fuel as a heat sink. Using the National Combustion Code a PICS-based combustor was modeled with liquid fuel at the supersonic cruise condition.
    [Show full text]
  • 2016 KIVA-Hpfe Development: a Robust and Accurate Engine Modeling Software David Carrington Los Alamos National Laboratory June 7, 2016 3:15 P.M
    2016 DOE Merit Review 2016 KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software David Carrington Los Alamos National Laboratory June 7, 2016 3:15 p.m. This presentation does not contain any proprietary, confidential, or otherwise restricted information LA-UR-16-22144 Project ID # ACE014 UNCLASSIFIED Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA 2016 DOE Overview Merit Review Timeline Barriers . Improve understanding of the fundamentals of fuel injection, fuel-air mixing, thermodynamic . 10/01/09 combustion losses, and in-cylinder combustion/ emission formation processes over a range of . combustion temperature for regimes of interest by 9/31/18 adequate capability to accurately simulate these processes . 85% complete (the last 10% in CFD . Engine efficiency improvement and engine-out development takes the most effort) emissions reduction . Minimization of time and labor to develop engine technology – User friendly (industry friendly) software, robust, accurate, more predictive, & quick meshing Budget Partners • Total project funding to date: • Dr. Jiajia Waters, LANL – 4000K • University of New Mexico- Dr. Juan Heinrich – 705K in FY 16 • University of Purdue, Calumet - Dr. Xiuling Wang – Contractor (Universities) share • Informal collaboration Reactive-Design/ANSYS ~15% • Informal collaboration GridPro Inc. Slide 2 UNCLASSIFIED Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA Objectives 2016 DOE FY 10 to FY 15 KIVA-Development Merit Review
    [Show full text]