CN114004179A - Heat release rate rapid prediction method of marine diesel engine based on phenomenological process - Google Patents
Heat release rate rapid prediction method of marine diesel engine based on phenomenological process Download PDFInfo
- Publication number
- CN114004179A CN114004179A CN202111299997.1A CN202111299997A CN114004179A CN 114004179 A CN114004179 A CN 114004179A CN 202111299997 A CN202111299997 A CN 202111299997A CN 114004179 A CN114004179 A CN 114004179A
- Authority
- CN
- China
- Prior art keywords
- fuel
- formula
- cylinder
- diesel engine
- rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
The invention provides a method for quickly predicting the heat release rate of a marine diesel engine based on a phenomenological process, which is based on a phenomenological modeling mechanism, considers penetration of spraying, air entrainment, evaporation and combustion from the moment of oil injection, reduces the real combustion process of the diesel engine by a formula as simple as possible, and realizes quick prediction of the heat release rate. According to the adjustment of the oil injection condition and the geometric parameters of the diesel engine, the model prediction method is not limited to the diesel engine with a specific model and a specific operation condition. Besides the heat release rate, the method can also predict other combustion characteristic parameters such as a stagnation period, cylinder pressure, air entrainment rate and the like.
Description
Technical Field
The invention relates to the field of diesel engine combustion performance analysis, in particular to a method for quickly predicting heat release rate of a marine diesel engine based on a phenomenological process.
Background
The marine diesel engine is one of main power devices of various large ocean-going ships in the world, and the performance of the marine diesel engine has a great influence on the benefit of the shipping industry. The heat release rate of the diesel engine is one of core parameters representing combustion and working performance, a method capable of rapidly predicting the heat release rate is found, and the method plays an important role in improving the transient performance prediction capability and the operation reliability of the marine diesel engine.
At present, the method for predicting the heat release rate of the diesel engine mainly comprises zero-dimensional model prediction and CFD model prediction. The zero-dimensional model is simple to build, the prediction speed and the response are high, but the combustion part is replaced by a Weber or double Weber semi-empirical formula, so that the physicochemical process in the diesel engine cannot be reflected, and the predicted heat release rate is lack of authenticity. The CFD model has a detailed structure which is most similar to the working process of a real diesel engine, but because the factors considered in calculation are excessive, the calculation load is large, the operation efficiency is low, and the rapid prediction of the heat release rate is difficult to realize.
Disclosure of Invention
The invention aims to realize a rapid heat release rate prediction method of a marine diesel engine based on a phenomenological process, which is used for rapidly calculating the heat release rate and predicting the performance of a working condition based on a real spray development and combustion mechanism and providing a better numerical analysis method for transient performance prediction of the marine diesel engine.
A method for rapidly predicting the heat release rate of a marine diesel engine based on a phenomenological process comprises the following steps:
step 1, determining boundary conditions and initial conditions according to geometric parameters and working condition parameters of the diesel engine, dividing the whole working process of the diesel engine into a plurality of small calculation steps from the beginning of compression to the end of combustion, and calculating the temperature change in the cylinder of each crank angle step from the beginning of compression to the beginning of oil injection by adopting an energy conservation equation
In the formulaIs the in-cylinder heat transfer amount in one step,volume work for the working medium, m mass of the working medium, cv1Is the specific heat capacity of the working medium.
According to the above formula, temperature T according to initial conditions1Solving for the temperature T of the end of the first step2Will T2And performing iterative calculation on the initial temperature serving as the next step to obtain the temperature of each step length from the compression to the oil injection starting moment. Knowing the temperature of each step, the pressure p can be obtained from the ideal gas state equation pV ═ RT for each step.
Step 2, after oil injection is started, calculating the spray penetration distance S and the air entrainment rate m in each time step by adopting the following formulaa:
In the formula,. DELTA.pinjIs the difference between the injection pressure and the in-cylinder pressure, rhogIs the density of the working medium in the cylinder, t is the time from the beginning of the injection to the calculation step length, dnozIs the diameter of the orifice, T is the temperature in the cylinder, mfIs the mass flow of the injected fuel.
Step 3, calculating the change of the turbulence energy k in the cylinder brought by the development of the spray in the step 2:
wherein M represents the total mass of the spray zone, ujThe nozzle outlet speed belongs to turbulence energy dissipation rate, and the following formula is adopted for calculation:
in the formula, LjFor the in-cylinder turbulence energy scale, the following equation is used to calculate:
in the formula, ρfIs the fuel density.
And 4, calculating the evaporation rate of the spray in each step, wherein only the evaporated spray can participate in final combustion:
in the formula, raIs the rate of heat absorption per unit mass of fuel, p is the in-cylinder pressure, rvM is the evaporation rate of fuel oilliqAnd E is the mass of the liquid fuel in the current step length, and the heat required by the evaporation of the unit mass of the fuel.
Step 5, calculating the stagnation period tau by adopting the following formulai:
In the formula, SpThe average speed of the piston under the current working condition, R is a gas constant, EaThe activation energy of the fuel can be calculated by the following formula:
Ea=618840/(CN+25)
in the formula, CN is the octane number of the fuel oil.
And 6, calculating the combustion rate of the evaporated fuel in the step 4, wherein a laminar-turbulent characteristic time model is adopted in the calculation process as follows:
wherein x represents the mass fraction of fuel, τcIs the characteristic time of combustion. Characteristic time taucFrom laminar characteristic time τlWith characteristic time τ of turbulencetAnd a delay factor f consisting of:
τc=τl+fτt
characteristic time tau of laminar flowlCalculated using the formula:
wherein A is a calibration coefficient, [ Fuel]Is the molar fraction of fuel oil, [ O ]2]Is the oxygen molar fraction.
Characteristic time of turbulence τtThen it is calculated from the turbulence energy k calculated in step 3 and the turbulence dissipation factor ∈:
τt=0.1k/∈
the delay coefficient f is calculated by the formula:
where r is the ratio of the mass of all combustion products to the mass of all reactants, expressed as:
and 7, calculating the heat release rate dQ/dt of the whole cylinder according to the combustion rate calculated in the step 6:
where LHV is the lower heating value of the fuel, NnozThe number of the spray holes on the diesel injector.
Compared with the prior art, the invention has the beneficial effects that:
the method is based on a phenomenological modeling mechanism, considers the penetration of spray, air entrainment, evaporation and combustion from the moment of oil injection, reduces the real combustion process of the diesel engine by a formula as simple as possible, and realizes the rapid prediction of the heat release rate. According to the adjustment of the oil injection condition and the geometric parameters of the diesel engine, the model prediction method is not limited to the diesel engine with a specific model and a specific operation condition. Besides the heat release rate, the method can also predict other combustion characteristic parameters such as a stagnation period, cylinder pressure, air entrainment rate and the like.
Drawings
FIG. 1 is a flow chart of a prediction method of the present invention;
FIG. 2 is a graph comparing the calculated results and experimental results for the cylinder pressure and heat release rate of a diesel engine of a certain model;
FIG. 3 shows parameters for a model of diesel engine of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
Aiming at the problems that the existing heat release rate prediction method of the diesel engine is difficult to reflect a real process, the calculation speed is too slow, the transient prediction performance is poor and the like, the method is based on a phenomenological modeling theory, considers the processes of penetration of spray, air entrainment, evaporation and combustion from the moment of oil injection, and provides the heat release rate rapid prediction method based on the phenomenological process. The method has the advantages of high calculation speed, high precision and wide application range, and is convenient for predicting and evaluating the transient working performance of the marine diesel engine.
The purpose of the invention is realized by the following technical scheme:
step 1: according to geometric parameters and working conditions of the diesel engineThe parameters determine boundary conditions and initial conditions, the whole diesel engine working process is divided into a plurality of small calculation steps from the beginning of compression to the end of combustion, and the following energy conservation equation is adopted to calculate the temperature change in the cylinder of each crank angle step from the beginning of compression to the beginning of oil injection
In the formulaIs the in-cylinder heat transfer amount in one step,volume work for the working medium, m mass of the working medium, cv1Is the specific heat capacity of the working medium. According to the above formula, temperature T according to initial conditions1Solving for the temperature T of the end of the first step2Will T2And performing iterative calculation on the initial temperature serving as the next step to obtain the temperature of each step length from the compression to the oil injection starting moment. Knowing the temperature of each step, the pressure p can be obtained from the ideal gas state equation pV ═ RT for each step.
Step 2: after the oil injection is started, the spray penetration distance S and the air entrainment rate m in each time step are calculated by adopting the following formulaa:
In the above formula,. DELTA.pinjIs the difference between the injection pressure and the in-cylinder pressure, rhogThe density of the working medium in the cylinder, and t is the time from the beginning of oil injection to the calculation step lengthM, dnozIs the diameter of the orifice, T is the temperature in the cylinder, mfIs the mass flow of the injected fuel.
And step 3: calculating the change of the turbulence energy k in the cylinder brought by the development of the spray in the step 2:
wherein M represents the total mass of the spray zone, ujThe nozzle outlet speed belongs to turbulence energy dissipation rate, and the following formula is adopted for calculation:
in the formula, LjFor the in-cylinder turbulence energy scale, the following equation is used to calculate:
in the formula, ρfIs the fuel density.
And 4, step 4: the evaporation rate of the spray in each step is calculated, and only the evaporated spray can participate in the final combustion:
in the formula, raIs the rate of heat absorption per unit mass of fuel, p is the in-cylinder pressure, rvM is the evaporation rate of fuel oilliqAnd E is the mass of the liquid fuel in the current step length, and the heat required by the evaporation of the unit mass of the fuel.
And 5: the combustion lag period tau is calculated by the following formulai:
In the formula, SpThe average speed of the piston under the current working condition, R is a gas constant, EaThe activation energy of the fuel can be calculated by the following formula:
Ea=618840/(CN+25)
in the formula, CN is the octane number of the fuel oil. When the in-cylinder crank angle exceeds the calculated value of the stagnation period formula, the fuel evaporated in step 4 starts to burn.
Step 6: and (4) calculating the combustion rate of the evaporated fuel in the step (4), wherein a laminar-turbulent characteristic time model is adopted in the calculation process as follows:
wherein x represents the mass fraction of fuel, τcIs the characteristic time of combustion. Characteristic time taucFrom laminar characteristic time τlWith characteristic time τ of turbulencetAnd a delay factor f consisting of:
τc=τl+fτt
characteristic time tau of laminar flowlCalculated using the formula:
wherein A is a calibration coefficient, [ Fuel]Is the molar fraction of fuel oil, [ O ]2]Is the oxygen molar fraction.
Characteristic time of turbulence τtThen it is calculated from the turbulence energy k calculated in step 3 and the turbulence dissipation factor ∈:
τt=0.1k/∈
the delay coefficient f is calculated by the formula:
where r is the ratio of the mass of all combustion products to the mass of all reactants, expressed as:
step 7, calculating the heat release rate dQ/dt of the entire cylinder from the combustion rate calculated in step 6:
where LHV is the lower heating value of the fuel, NnozThe number of the spray holes on the diesel injector.
The foregoing is merely a description of preferred embodiments of the present invention and does not limit the spirit and scope of the invention. Any modification, variation, improvement, etc. made by those skilled in the art based on the technical solution of the present invention should fall into the protection scope of the present invention without departing from the innovative concept of the present invention and the protection scope of the claims. The technical content of the invention is fully described in the claims.
Claims (1)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111299997.1A CN114004179B (en) | 2021-11-04 | 2021-11-04 | Rapid prediction method for heat release rate of marine diesel engine based on phenomenological process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111299997.1A CN114004179B (en) | 2021-11-04 | 2021-11-04 | Rapid prediction method for heat release rate of marine diesel engine based on phenomenological process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114004179A true CN114004179A (en) | 2022-02-01 |
| CN114004179B CN114004179B (en) | 2024-06-04 |
Family
ID=79927334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111299997.1A Active CN114004179B (en) | 2021-11-04 | 2021-11-04 | Rapid prediction method for heat release rate of marine diesel engine based on phenomenological process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114004179B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114962088A (en) * | 2022-06-09 | 2022-08-30 | 北京理工大学 | A turbulent-phase transition synergistic supercritical combustion enhancement method |
| CN118428185A (en) * | 2024-04-28 | 2024-08-02 | 哈尔滨工程大学 | A calculation method for the spray combustion process of a high-pressure direct injection dual-fuel marine internal combustion engine |
| CN119761023A (en) * | 2024-12-17 | 2025-04-04 | 哈尔滨工程大学 | A fast calculation method and system for cylinder pressure of compression ignition methanol engine based on phenomenological spray characteristics |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002195084A (en) * | 2000-12-25 | 2002-07-10 | Nippon Soken Inc | Fuel injection control system of diesel engine |
| CN1871416A (en) * | 2003-10-01 | 2006-11-29 | 韦斯特波特研究公司 | Method and apparatus for controlling combustion quality of a gaseous-fuelled internal combustion engine |
| CN101571090A (en) * | 2009-06-04 | 2009-11-04 | 北京航空航天大学 | Fuel injection rule measuring device of diesel engine and measuring method thereof |
| US20120004826A1 (en) * | 2010-06-30 | 2012-01-05 | Mazda Motor Corporation | Diesel engine and method of controlling the diesel engine |
| CN112304623A (en) * | 2020-10-28 | 2021-02-02 | 哈尔滨工程大学 | Effective thermal efficiency prediction method of marine diesel engine based on fuel components |
-
2021
- 2021-11-04 CN CN202111299997.1A patent/CN114004179B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002195084A (en) * | 2000-12-25 | 2002-07-10 | Nippon Soken Inc | Fuel injection control system of diesel engine |
| CN1871416A (en) * | 2003-10-01 | 2006-11-29 | 韦斯特波特研究公司 | Method and apparatus for controlling combustion quality of a gaseous-fuelled internal combustion engine |
| CN101571090A (en) * | 2009-06-04 | 2009-11-04 | 北京航空航天大学 | Fuel injection rule measuring device of diesel engine and measuring method thereof |
| US20120004826A1 (en) * | 2010-06-30 | 2012-01-05 | Mazda Motor Corporation | Diesel engine and method of controlling the diesel engine |
| CN112304623A (en) * | 2020-10-28 | 2021-02-02 | 哈尔滨工程大学 | Effective thermal efficiency prediction method of marine diesel engine based on fuel components |
Non-Patent Citations (2)
| Title |
|---|
| 刘博;韩长福;丁虎;刘龙: "重油现象学喷雾模型建模研究", 船舶工程, no. 0, 31 December 2019 (2019-12-31) * |
| 孟维;邓康耀: "直喷柴油机燃烧的现象学快速预测模型", 柴油机, vol. 41, no. 005, 31 December 2019 (2019-12-31) * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114962088A (en) * | 2022-06-09 | 2022-08-30 | 北京理工大学 | A turbulent-phase transition synergistic supercritical combustion enhancement method |
| CN114962088B (en) * | 2022-06-09 | 2023-10-20 | 北京理工大学 | A turbulence-phase change synergistic supercritical combustion enhancement method |
| CN118428185A (en) * | 2024-04-28 | 2024-08-02 | 哈尔滨工程大学 | A calculation method for the spray combustion process of a high-pressure direct injection dual-fuel marine internal combustion engine |
| CN118428185B (en) * | 2024-04-28 | 2025-02-11 | 哈尔滨工程大学 | A calculation method for the spray combustion process of a high-pressure direct injection dual-fuel marine internal combustion engine |
| CN119761023A (en) * | 2024-12-17 | 2025-04-04 | 哈尔滨工程大学 | A fast calculation method and system for cylinder pressure of compression ignition methanol engine based on phenomenological spray characteristics |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114004179B (en) | 2024-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114004179B (en) | Rapid prediction method for heat release rate of marine diesel engine based on phenomenological process | |
| Xin et al. | Computer modeling of diesel spray atomization and combustion | |
| Sharma et al. | Effect of fuel injection pressure and engine speed on performance, emissions, combustion, and particulate investigations of gasohols fuelled gasoline direct injection engine | |
| Liu et al. | An examination of performance deterioration indicators of diesel engines on the plateau | |
| Zhou et al. | Scaling fuel sprays for different size diesel engines | |
| Zhou et al. | The similarity ratio effects in design of scaled model experiments for marine diesel engines | |
| Feng et al. | Effect of nozzle geometry on combustion of a diesel-methanol dual-fuel direct injection engine | |
| Liu et al. | Investigation of the effect of simulated atmospheric conditions at different altitudes on the combustion process in a heavy-duty diesel engine based on zero-dimensional modeling | |
| Yan et al. | Numerical simulation study of spray wall impingement and adhered fuel film under different altitude conditions during cold start of diesel engines | |
| Jia et al. | A micro-variable circular orifice fuel injector for HCCI-conventional engine combustion-Part I numerical simulation of cavitation | |
| Richards et al. | Modeling the effects of EGR and injection pressure on emissions in a high-speed direct-injection diesel engine | |
| Guangxin et al. | Effects of fuel temperature on injection process and combustion of dimethyl ether engine | |
| Lee et al. | Effects of injection strategies on the flow and fuel behavior characteristics in port dual injection engine | |
| Mohan et al. | Numerical simulation on spray characteristics of ether fuels | |
| Geng et al. | Numerical study of effect of post injection coupled with EGR on combustion and emission performances of CRDI engine fueled with biodiesel-ethanol blends | |
| CN114810447A (en) | A test device and test method for gas-liquid two-phase ratio during transient injection of ammonia fuel with phase change | |
| Das et al. | Factors affecting heat transfer in a diesel engine: low heat rejection engine revisited | |
| Hountalas et al. | Phenomenological modelling of oxygen-enriched combustion and pollutant formation in heavy-duty diesel engines using exhaust gas recirculation | |
| Hountalas et al. | Validation of multi-zone combustion model ability to predict two stroke diesel engine performance and NOx emissions using on board measurements | |
| Liu et al. | Effects of Scavenging-Port Semi-Direct Injection Strategy on Mixture Formation for a Two-Stroke Spark-Ignition Kerosene Piston Engine | |
| Goryntsev et al. | Towards large eddy simulation of spray combustion in direct injection spark ignition engine | |
| Ghahremani et al. | Experimental and theoretical investigation on spray characteristics of bio-ethanol blends using a direct injection system | |
| Bella et al. | Combustion and spray simulation of a DI turbocharged diesel engine | |
| Kowalski | The theoretical study on influence of fuel injection pressure on combustion parameters of the marine 4-stroke engine | |
| Zhang et al. | Simulation research on matching of spray and combustion chamber geometry in diesel engine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |