1. Introduction
Approximately 20–40% of the total provided energy after combustion is lost in the exhaust gasses [1,2]. However, the residual energy content of the exhaust gasses can be split into two parts, thermal energy and pressure energy. The exhaust gasses in the cylinder after combustion are hot and the gasses contain therefore potential usable energy, which can be extracted in a further process [3]. Thus, turbocharger are often used to expand in part the heat and the pressure potential [4], in order to increase the specific engine efficiency. The exhaust stroke of a four stroke engine can be divided into two main periods, the blow-down and the scavenging pulse. The pressure difference between the in-cylinder pressure and the pressure in the exhaust tract forces the exhaust gas expulsion at the
early stage, which is referred as the blow-down phase. When the valves start to open, only a small slid between the valve head and the valve seat represents the discharge passage. Nevertheless, during the blow-down phase, the mass flow rate is high and decreases rapidly with the pressure dropping in the cylinder pod. Thereafter, when the pressure difference between the cylinder chamber and the exhaust manifold is balanced, the scavenging phase takes place, in which the piston cleans out the rest of the residual gas by its motion towards the top dead center. During this phase, the mass flow rate is lower than during the blow-down phase. However, the response of the mass flow rate during the scavenging phase depends highly on the pressure conditions downstream in the collecting manifold, the possible interactions with other cylinders over the manifold, and the back-pressure provoked by the turbocharger. The turbine of the turbocharger extractsflow energy and generatestherebyaresistancetotheflowdownstreamoftheexhaustport. Theflowresistanceinducesaback-pressure,which complicates the
gasexpulsionfromthecylindersincethisresultsinadditionalwork demanded from the piston. In order to reduce the back-pressure effect of the turbocharger, several approaches can be employed. A wastegatetobypasstheturbochargerduringthescavengingphase canbeusedtoreducetheback-pressure.Withvariablevalveactuation and divided exhaust periods two separated exhaust ports are used for the blow-down phase and the scavenging phase, where thevalvesareactuatedatdifferenttimes[5].Variablevalvetiming, i.e.regulatingadaptivelythevalveopeningtime,the openingduration,andthevalvelift,iscrucialfortheengineperformance[6]and canbeusedtoreducethefuelconsumptionoftheinternalcombustion engine [7]. However, the alteration of the valve timing influences the emission production. Simulating the entire process of an internal combustion engine with all degrees of freedom is computationally unaffordable. Therefore, fast simplified models areutilizedfortheinitialdesignprocessoftheinternalcombustion engine[8].However,thephysicsandthebehavioroftheenginemust be represented by the model. An oversimplified model may not represent the real engine any more. The inflow conditions into the turbocharger turbine can crucially effect its performance [9]. Usually, the valve timing and the valve opening-speed is studied and optimized by one-dimensional simulations, where the quality of the flow field and therefore the inflow conditions into the turbochargerturbine arenot accountedfor. Inone-dimensionalmodelingofthe internalcombustionengine process, tables, databases, and operation maps are used to describe the performance behavior of the individual constituents, such as the piping system, constrictions, or bends [10,11]. The tables or databases are based on empirical formulas with experimentally evaluated coefficients. For the exhaust port, the discharge coefficient is evaluated and employed in the one-dimensional engine analysis to evaluate the total pressure drop [12]. A so-called flow bench experiment is performed, where the exhaust port discharge coefficient is evaluated at fixed valve lifts and constant specified total pressure drops. These experiments are usually performed on the real cylinder head at room temperatures. Hence, the flow streaming through the port is cold in contrast to a real engine case [13]. The discharge coefficient of sharp edged flow constrictions was evaluated experimentally [14], with the conclusion that the steady discharge coefficients are generally lower than those under unsteady flow conditions. The effect of pulsating flow on the exhaust port flow coefficients has been investigated numerically for different geometries regarding the blow down pulse and room temperatures [15]. It was found that the maximum difference in the flow coefficient was at most a 6% increase or a 7% decrease. The averaged flow coefficient over several valve lifts varied between 0:5% and 2:5% for the different geometries, which resembles a low impact on the total engine performance. Therefore, Bohac and Landfahrer [15] conclude that flow bench measurements at constant valve lift with constant mass flow rate are adequate. Continuous ordinary differential equations can be used in onedimensional simulation codes for more accurate modeling of engine components [16] or power plants [17]. The transfer functions for a given geometry can be obtained by impulse excitation to determine the characteristics of an engine transmission component. An experimental approach relating the mass flow rate with thedynamicpressurefluctuationscanbeusedtoabstractanordinary differential equation model for an engine part [18]. A flow based model can be obtained reducing the Navier–Stokes equations via Galerkin projection onto representative flow modes [19]. Thus, the original governing system of partial differential equations is replaced by a set of ordinary differential equations, which are computationally inexpensive to compute. The Proper Orthogonal Decomposition (POD) method decomposes the flow field into a series of modes, where the modes are constructed such
that the least number of modes is needed to reproduce the kinetic energy of the flow field [20]. The POD modes are also used to characterize the flow [21] and identify coherent structures with large energy contend. Also Fourier mode decomposition or dynamic mode decomposition have been used for constructing reduced order models [22], where certain frequencies are important [23]. The experimental assessment of the flow field in a complex, confined geometry, such as the exhaust port, is a challenging task [24]. However, experimental flow visualizations of the flow field development inside complex geometries, such as a closed cylinder with moving piston and combustion, have been performed [25]. Further, experimental flow visualization in an exhaust port can be performed with a considerable effort and assumptions [26], such as fixed valve lifts. Nevertheless, numerical computations of the flow in a complex geometry with non-moving boundaries can be easily performed. The computational effort of a numerical simulation depends on the complexity, amount of modeling and flow scales considered. A time-averaged simulation approach, such as steady state Reynolds Averaged Navier–Stokes (RANS) simulations, can be achieved with a rather low computational afford. However, the required time resolution of the flow phenomena for POD is not computed. To resolve all flow scales, i.e. Direct Numerical Simulation (DNS), is computationally expensive and the present purpose does not justify this the excessive usage of resources. With a reasonably resolved LES simulation, a substantial proportion of the inertial subrange in the turbulence spectra is resolved and only the smallest flow scales are modeled. The smallest flow scales dissipate the flow energy to heat and this behavior is of universal character, independent of the geometry. Generally, large energetic turbulent structures imply a substantial proportion of turbulent dissipation [27]. Hence, the LES approach reassembles a suitable trade between accuracy, reliability, and effort, which has beensuccessfully appliedto optimizetheportgeometryofinternal combustion engines [28]. Many numerical investigations analyze the flowinthe intakeportand theconsequencesonthe in-cylinder flow. Additionally, the flow in the ports of an internal combustion engine with pulsating boundary conditions has been simulated, e.g.[28].However,themoststudiesemphasizeonthecomputation models, performance parameter, or optimization of the geometry ratherthan givinginsight into the generated flow structures occurring in the exhaust port. This analysis is part of a larger investigation study treating the gas exchange process in an internal combustion engine. Within this paper the importance of dynamic boundary conditions on the flow generated losses in a realistic exhaust port of an internal combustion engine is investigated numerically. The three-dimensional compressible Navier–Stokes equations are simulated using an LES approach. A comparison of two cases, i.e. using constant inflow and outflow boundary conditions, and using engine cycle dependent inflow boundary conditions, at a constant valve lift are performed. These assumptions allow to isolate the effects caused due to flow pulsation. Therefore, the focus of contrasting the cases is held on the quantities, which are characteristic for the flow induced losses, as e.g. turbulent dissipation and friction losses. The change of the coherent flow structures are analyzed using the POD method.
1. บทนำประมาณ 20 – 40% ของผลรวมที่ให้พลังงานจากการเผาไหม้จะหายไปใน gasses ไอเสีย [1, 2] อย่างไรก็ตาม เนื้อหาพลังงานส่วนที่เหลือของ gasses ไอสามารถแบ่งเป็นสองส่วน พลังงานความร้อน และพลังงานความดัน Gasses ไอเสียในกระบอกหลังจากเผาไหม้ร้อนและ gasses ประกอบด้วยดังนั้นอาจใช้พลังงาน ซึ่งสามารถแยกกระบวนการต่อไป [3] ดังนั้น เทอร์โบมักใช้เพื่อขยายในส่วนความร้อนและความดันอาจเกิดขึ้น [4], เพื่อเพิ่ม efficiency เครื่องยนต์ specific จังหวะไอเสียของเครื่องยนต์สี่จังหวะสามารถแบ่งออกเป็นสองรอบระยะเวลาที่หลัก ลงเป่า และชีพจร scavenging ผลต่างความดันระหว่างความดันในลูกสูบและความดันในระบบทางเดินไอเสียขับก๊าซไอเสียที่บังคับการระยะแรก ๆ ซึ่งอ้างอิงเป็นระยะระเบิดลง เมื่อวาล์วเริ่มเปิด ขนาดเล็กฝ่อระหว่างหัววาล์ว และวาล์วนั่งแสดงกาลปล่อย อย่างไรก็ตาม ระยะเป่าลง อัตรามวล flow สูง และลดลงอย่างรวดเร็ว ด้วยความดันที่ปล่อยในฝักรูปทรงกระบอก หลังจากนั้น เมื่อความสมดุลความแตกต่างความดันระหว่างหอถังและอเนกไอเสีย ระยะ scavenging เกิด ซึ่งลูกสูบทำความสะอาดส่วนเหลือของก๊าซส่วนที่เหลือ โดยการเคลื่อนไหวไปยังศูนย์ตายบน ในระหว่างขั้นตอนนี้ อัตรา flow โดยรวมได้ต่ำกว่าระยะระเบิดลง อย่างไรก็ตาม การตอบสนองของอัตรามวล flow ระยะ scavenging นั้นสูงความดันสภาพน้ำในการเก็บรวบรวมช่องทางเข้าออก การโต้ตอบได้กับภาชนะบรรจุอื่น ๆ มากมาย และดันกลับ provoked โดยเทอร์โบ กังหันเทอร์โบ extractsflow พลังงานและ generatestherebyaresistancetotheflowdownstreamoftheexhaustport Theflowresistanceinducesaback-ความดัน ที่ complicates ในgasexpulsionfromthecylindersincethisresultsinadditionalwork ต้องจากลูกสูบ เพื่อลดผลดันกลับของเทอร์โบ สามารถทำงานหลายแนว การ wastegatetobypasstheturbochargerduringthescavengingphase canbeusedtoreducetheback ดัน ใช้สำหรับเป่าลงเฟสและเฟส scavenging, Withvariablevalveactuation และพอร์ตไอเสียแยกสองรอบไอเสียถูกแบ่งที่ thevalvesareactuatedatdifferenttimes [5] Variablevalvetiming, i.e.regulatingadaptivelythevalveopeningtime,the openingduration, andthevalvelift, iscrucialfortheengineperformance [6] และ canbeusedtoreducethefuelconsumptionoftheinternalcombustion เครื่องยนต์ [7] อย่างไรก็ตาม เปลี่ยน influences เวลาวาล์วผลิตมลพิษ การจำลองกระบวนการทั้งหมดของเครื่องยนต์สันดาปภายในกับทุกองศาความเป็นอิสระเป็น computationally ข้าวยากหมากแพง ดังนั้น อย่างรวดเร็ว simplified รุ่น areutilizedfortheinitialdesignprocessoftheinternalcombustion เครื่องยนต์ [8] อย่างไรก็ตาม thephysicsandthebehavioroftheenginemust แสดงตามแบบ แบบจำลอง oversimplified อาจไม่แสดงเครื่องยนต์จริงอีกด้วย เงื่อนไข inflow เป็นกังหันเทอร์โบยังสามารถมีผลประสิทธิภาพ [9] ปกติ เวลาวาล์วและวาล์วเปิดความเร็วในการศึกษา และเพิ่มประสิทธิภาพ โดยจำลอง one-dimensional เงื่อนไขที่คุณภาพ flow field และ inflow ที่เป็น turbochargerturbine arenot accountedfor Inone dimensionalmodelingofthe internalcombustionengine กระบวนการ ตาราง ฐานข้อมูล และการดำเนินการแผนที่จะใช้ในการอธิบายลักษณะการทำงานประสิทธิภาพของ constituents ละ เช่นระบบท่อ constrictions หรือจัดฟัน [10,11] ตารางหรือฐานข้อมูลอยู่บนสูตร empirical coefficients experimentally ค่า สำหรับพอร์ตไอเสีย coefficient จำหน่ายจะถูกประเมิน และทำงานในการวิเคราะห์เครื่องยนต์ one-dimensional ประเมินความดันรวมหล่น [12] ทดลองผู้พิพากษาเรียกว่า flow จะดำเนินการ ที่ coefficient จำหน่ายพอร์ตไอเสียจะถูกประเมินในลิฟท์ fixed วาล์ว และหยด specified คงความดันรวม ทดลองเหล่านี้มักจะดำเนินในรูปทรงกระบอก head จริงที่อุณหภูมิห้อง ดังนั้น flow สตรีมมิ่งผ่านพอร์ตเป็นเย็นตรงข้ามกรณีเครื่องยนต์จริง [13] มีประเมิน coefficient ปล่อยของชาร์ป edged flow constrictions experimentally [14], กับข้อสรุปที่ coefficients ปล่อยมั่นคงโดยทั่วไปต่ำกว่าสภาวะ unsteady flow ผลของการขยับ flow coefficients flow พอร์ตไอเสียได้ถูกสอบสวนเรียงตามตัวเลขในรูปทรงเรขาคณิตที่แตกต่างกันเกี่ยวกับระเบิดชีพจรและอุณหภูมิห้อง [15] จะพบว่า ความแตกต่างสูงสุดใน flow coefficient ถูกที่สุด 6% เพิ่มขึ้นหรือลดลง 7% Coefficient flow เฉลี่ยกว่าลิฟท์วาล์วหลายแตกต่างกันระหว่าง 0:5%และ% 2:5 สำหรับรูปทรงเรขาคณิตที่แตกต่างกัน ซึ่งคล้ายกับผลกระทบต่อประสิทธิภาพการทำงานของเครื่องยนต์รวมต่ำ ดังนั้น Bohac และ Landfahrer [15] สรุป flow วัดม้าที่ยกวาล์วคงที่มีอัตราคงโดยรวม flow เพียงพอ สามารถใช้ในรหัสการจำลอง onedimensional สมการเชิงอนุพันธ์สามัญอย่างต่อเนื่องสำหรับโมเดลถูกต้องมากขึ้น [16] ส่วนประกอบของเครื่องยนต์หรือไฟฟ้า [17] ฟังก์ชันโอนย้ายสำหรับเรขาคณิตที่กำหนดให้ได้ตามกระแสในการกระตุ้นเพื่อกำหนดลักษณะของคอมโพเนนต์เกียร์เครื่องยนต์ วิธีทดลองที่เกี่ยวข้องอัตรามวล flow กับ thedynamicpressurefluctuationscanbeusedtoabstractanordinary รูปแบบสมการเชิงอนุพันธ์สำหรับหนึ่ง [18] สามารถรับแบบ flow ที่ใช้ลดสมการ Navier-สโตกส์ผ่านฉาย Galerkin flow ตัวแทนโหมด [19] ดังนั้น ระบบควบคุมเดิมของสมการเชิงอนุพันธ์บางส่วนจะถูกแทนที่ ด้วยชุดปกติสมการเชิงอนุพันธ์ ซึ่งมีราคาไม่แพง computationally จะคำนวณ เหมาะสม Orthogonal แยกส่วนประกอบ (POD) วิธี decomposes flow field ในชุดโหมด ซึ่งมีวิธีการสร้างดังกล่าวthat the least number of modes is needed to reproduce the kinetic energy of the flow field [20]. The POD modes are also used to characterize the flow [21] and identify coherent structures with large energy contend. Also Fourier mode decomposition or dynamic mode decomposition have been used for constructing reduced order models [22], where certain frequencies are important [23]. The experimental assessment of the flow field in a complex, confined geometry, such as the exhaust port, is a challenging task [24]. However, experimental flow visualizations of the flow field development inside complex geometries, such as a closed cylinder with moving piston and combustion, have been performed [25]. Further, experimental flow visualization in an exhaust port can be performed with a considerable effort and assumptions [26], such as fixed valve lifts. Nevertheless, numerical computations of the flow in a complex geometry with non-moving boundaries can be easily performed. The computational effort of a numerical simulation depends on the complexity, amount of modeling and flow scales considered. A time-averaged simulation approach, such as steady state Reynolds Averaged Navier–Stokes (RANS) simulations, can be achieved with a rather low computational afford. However, the required time resolution of the flow phenomena for POD is not computed. To resolve all flow scales, i.e. Direct Numerical Simulation (DNS), is computationally expensive and the present purpose does not justify this the excessive usage of resources. With a reasonably resolved LES simulation, a substantial proportion of the inertial subrange in the turbulence spectra is resolved and only the smallest flow scales are modeled. The smallest flow scales dissipate the flow energy to heat and this behavior is of universal character, independent of the geometry. Generally, large energetic turbulent structures imply a substantial proportion of turbulent dissipation [27]. Hence, the LES approach reassembles a suitable trade between accuracy, reliability, and effort, which has beensuccessfully appliedto optimizetheportgeometryofinternal combustion engines [28]. Many numerical investigations analyze the flowinthe intakeportand theconsequencesonthe in-cylinder flow. Additionally, the flow in the ports of an internal combustion engine with pulsating boundary conditions has been simulated, e.g.[28].However,themoststudiesemphasizeonthecomputation models, performance parameter, or optimization of the geometry ratherthan givinginsight into the generated flow structures occurring in the exhaust port. This analysis is part of a larger investigation study treating the gas exchange process in an internal combustion engine. Within this paper the importance of dynamic boundary conditions on the flow generated losses in a realistic exhaust port of an internal combustion engine is investigated numerically. The three-dimensional compressible Navier–Stokes equations are simulated using an LES approach. A comparison of two cases, i.e. using constant inflow and outflow boundary conditions, and using engine cycle dependent inflow boundary conditions, at a constant valve lift are performed. These assumptions allow to isolate the effects caused due to flow pulsation. Therefore, the focus of contrasting the cases is held on the quantities, which are characteristic for the flow induced losses, as e.g. turbulent dissipation and friction losses. The change of the coherent flow structures are analyzed using the POD method.
การแปล กรุณารอสักครู่..
