WO2024253047A1 - ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法 - Google Patents

ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法 Download PDF

Info

Publication number
WO2024253047A1
WO2024253047A1 PCT/JP2024/020120 JP2024020120W WO2024253047A1 WO 2024253047 A1 WO2024253047 A1 WO 2024253047A1 JP 2024020120 W JP2024020120 W JP 2024020120W WO 2024253047 A1 WO2024253047 A1 WO 2024253047A1
Authority
WO
WIPO (PCT)
Prior art keywords
knocking
prediction
fuel
engine control
prediction method
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.)
Ceased
Application number
PCT/JP2024/020120
Other languages
English (en)
French (fr)
Inventor
薫 丸田
雄飛 森井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tohoku University NUC
Original Assignee
Tohoku University NUC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tohoku University NUC filed Critical Tohoku University NUC
Priority to CN202480037229.4A priority Critical patent/CN121241194A/zh
Priority to JP2025526095A priority patent/JPWO2024253047A1/ja
Publication of WO2024253047A1 publication Critical patent/WO2024253047A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02D—CONTROLLING COMBUSTION ENGINES
    • F02D45/00—Electrical control not provided for in groups F02D41/00 - F02D43/00

Definitions

  • the present disclosure relates to a knocking prediction method, a knocking prediction device, an engine control device, a vehicle, a method for designing a combustion device, a method for designing a power generating device, a method for designing a vehicle, and a method for designing a fuel.
  • Knocking is a phenomenon accompanied by shock waves generated by the spontaneous ignition of unburned gas in the combustion chamber, and can damage the engine.
  • the pressure and temperature (compression ratio) of the combustion chamber are lowered to prevent knocking from occurring.
  • the Livengood-Wu integral is a well-known method for predicting knocking. This method estimates the ignition delay time when the pressure and temperature in the combustion chamber are changed, and integrates it over time. When the time integral value reaches "1", autoignition, or knocking, occurs.
  • the Livengood-Wu integral is widely used in the design of actual engines, as shown in Patent Documents 1 and 2 below.
  • the "Livengood-Wu integral" is a simple method that can predict the occurrence of knocking (autoignition) to a certain extent, but the accuracy of the prediction is not necessarily sufficient. In order to make the most of the capabilities of the internal combustion engine as a power generation device, it is highly desirable to establish a knocking prediction method with higher accuracy than the "Livengood-Wu integral.”
  • the present disclosure has been made in consideration of the above-mentioned circumstances, and aims to provide a knocking prediction method and device, an engine control device, and a vehicle that can further improve accuracy.
  • the present disclosure aims to provide a design method for a combustion device, a power generating device, a vehicle, and a fuel that can improve the substantial performance of the fuel.
  • the knocking prediction method is a knocking prediction method for predicting knocking that may occur in a premixed combustion field, and predicts knocking based on the knowledge that in the premixed combustion field, the temporal changes in normalized fuel mass fraction and normalized temperature in zero-dimensional homogeneous ignition are equivalent to the spatial changes in a one-dimensional laminar premixed flame.
  • the knocking prediction method is a knocking prediction method for predicting knocking that may occur in a premixed combustion field, and in a characteristic diagram showing the relationship between normalized temperature and normalized fuel mass fraction of the fuel, a fuel with a Lewis number greater than 1 is deemed to have a possibility of knocking, and a fuel with a Lewis number less than 1 is deemed to have no possibility of knocking.
  • the knocking prediction method is a knocking prediction method for predicting knocking that may occur in a premixed combustion field, and includes a calculation step for sequentially calculating the combustion speed when the inlet temperature is gradually increased by one-dimensional steady-state analysis of the premixed combustion field, a determination step for determining whether the combustion speed cannot be obtained because a flame cannot exist in the premixed combustion field, and an occurrence condition acquisition step for acquiring the inlet temperature corresponding to the last acquired combustion speed as the occurrence condition of the knocking.
  • the knocking prediction method according to the fourth aspect of the present disclosure is the knocking prediction method according to the third aspect, in which the calculation step calculates the size of the area including the unburned mixture, the flame reaction zone, and the burned gas in the premixed combustion field by reducing it to the limit that the flame can contain.
  • the knocking prediction method according to the fifth aspect of the present disclosure is the knocking prediction method according to the third or fourth aspect, which predicts fuels with Lewis numbers greater than 1.
  • the knocking prediction device is a knocking prediction device that predicts knocking that may occur in a premixed combustion field, and includes a calculation means that sequentially calculates the combustion speed when the inlet temperature is gradually increased by one-dimensional steady-state analysis of the premixed combustion field, a determination means that determines whether the combustion speed cannot be obtained because a flame cannot exist in the premixed combustion field, and an occurrence condition acquisition means that acquires the inlet temperature corresponding to the last acquired combustion speed as the occurrence condition of the knocking.
  • the knocking prediction device is the knocking prediction device according to the first aspect, in which the calculation means calculates the size of the area including the unburned mixture, the flame reaction zone, and the burned gas in the premixed combustion field by reducing it to the limit that can be contained by the flame.
  • the knocking prediction device is the knocking prediction device according to the first or second aspect, in which the calculation means calculates fuels with Lewis numbers greater than 1.
  • the engine control device includes a knocking prediction device according to the first or second aspect, and generates an engine control signal based on the knocking prediction result of the knocking prediction device.
  • the engine control device is the engine control device according to the first aspect, further comprising a communication device capable of freely communicating with the outside, and transmits at least the knocking prediction result and the engine operation amount indicated by the engine control signal to the outside using the communication device.
  • the engine control device is the engine control device according to the second aspect, which transmits engine control variables and other engine-related information other than the engine control variables and the engine operation variables to the outside.
  • the vehicle according to the first aspect of the present disclosure includes an engine control device according to the first or second aspect, and an engine controlled by the engine control device.
  • the method for designing a combustion device according to the first aspect of the present disclosure is to design a combustion device based on the knocking occurrence conditions estimated based on the knocking prediction method according to the third or fourth aspect.
  • the method for designing a power generating device according to the first aspect of the present disclosure is to design a power generating device based on the knocking occurrence conditions estimated based on the knocking prediction method according to the third or fourth aspect.
  • the vehicle design method according to the first aspect of the present disclosure is to design a vehicle based on the knocking occurrence conditions estimated based on the knocking prediction method according to the third or fourth aspect.
  • the fuel design method according to the first aspect of the present disclosure is to design the fuel based on the knocking occurrence conditions estimated based on the knocking prediction method according to the third or fourth aspect.
  • a knocking prediction method it is possible to provide a knocking prediction method, a knocking prediction device, a combustion device design method, a power generation device design method, a vehicle design method, and a fuel design method that can further improve accuracy. Furthermore, according to the present disclosure, it is possible to provide a design method for a combustion device, a power generating device, a vehicle, and a fuel, which are capable of improving the substantial performance of the fuel.
  • FIG. 4 is a flowchart illustrating a knock prediction method according to an embodiment of the present disclosure.
  • 1 is a block diagram showing a functional configuration of a knocking prediction device according to an embodiment of the present disclosure. 1 is a list of various physical quantities used in the theoretical explanation of an embodiment of the present disclosure.
  • FIG. 4 is a first characteristic diagram showing a knocking prediction characteristic in the embodiment of the present disclosure.
  • FIG. 4 is a second characteristic diagram showing a knocking prediction characteristic in the embodiment of the present disclosure.
  • FIG. 11 is a third characteristic diagram showing a knocking prediction characteristic in the embodiment of the present disclosure.
  • FIG. 4 is a characteristic diagram showing a comparative example of a knocking prediction method according to an embodiment of the present disclosure.
  • 1 is a block diagram showing an engine control device and a vehicle according to an embodiment of the present disclosure.
  • FIG. 2 is a schematic diagram illustrating external communication of a vehicle according to an embodiment of the present disclosure.
  • FIG. 1 is a flowchart showing a knocking prediction method according to the present embodiment.
  • This flowchart shows a procedure for acquiring a flame propagation limit temperature, which is a condition for the occurrence of knocking, that is, an initial temperature at which knocking occurs, for a premixed combustion field (combustion chamber) in which a premixed mixture of fuel and air (oxidizer) is ignited by a compulsory means such as spark ignition and combusted.
  • This flowchart also shows an overview of the processing steps of a knocking prediction program (simulation program) executed by a specific computer.
  • the knocking prediction method according to this embodiment is realized by information processing (simulation) using a computer.
  • the computer that executes the knocking prediction method according to this embodiment i.e., the computer on which the knocking prediction program is installed, is the knocking prediction device according to this embodiment.
  • the knocking prediction device according to this embodiment is a device that outputs the conditions for knocking to occur in a premixed combustion field by executing the knocking prediction program.
  • FIG. 2 is a block diagram showing the functional configuration of the knocking prediction device A according to this embodiment.
  • the knocking prediction device A includes at least a storage unit 1 that stores a knocking prediction program in advance, a calculation unit 2 that executes the knocking prediction program, an operation unit 3 that inputs an operator's operation instructions to the calculation unit 2, and an output unit 4 that outputs the calculation result of the calculation unit 2, i.e., the above-mentioned flame propagation limit temperature (initial temperature), to the outside.
  • the above-mentioned flame propagation limit temperature initial temperature
  • memory unit 1, calculation unit 2, and operation unit 3 constitute the calculation means in this disclosure. Details will be described later, but memory unit 1, calculation unit 2, and operation unit 3 work together to sequentially calculate the combustion speed when the inlet temperature is gradually increased by one-dimensional steady-state analysis of the premixed combustion field.
  • the memory unit 1, the calculation unit 2, and the operation unit 3 in the knocking prediction device A constitute the determination means in this disclosure. That is, the memory unit 1, the calculation unit 2, and the operation unit 3 work together to determine whether or not it has become impossible to obtain the combustion speed because a flame can no longer exist in the premixed combustion field.
  • the memory unit 1, the calculation unit 2, and the operation unit 3 in the knocking prediction device A constitute the occurrence condition acquisition means in this disclosure. That is, the memory unit 1, the calculation unit 2, and the operation unit 3 cooperate to acquire the inlet temperature (flame propagation limit temperature) corresponding to the last acquired combustion speed as the occurrence condition of knocking.
  • the calculation means of the knocking prediction device A according to this embodiment corresponds to the calculation step in the knocking prediction method according to this embodiment.
  • the determination means of the knocking prediction device A according to this embodiment corresponds to the determination step in the knocking prediction method according to this embodiment.
  • the raw condition acquisition means of the knocking prediction device A according to this embodiment corresponds to the occurrence condition acquisition step in the knocking prediction method according to this embodiment.
  • the premixed combustion field in this embodiment is a combustion field that corresponds to the combustion chamber of an internal combustion engine (prime mover), such as a gasoline engine installed in a vehicle as a power generating device.
  • an internal combustion engine such as a gasoline engine installed in a vehicle as a power generating device.
  • measures such as lowering the compression ratio are taken to prevent knocking, that is, to prevent spontaneous ignition of unburned gas.
  • the above-mentioned internal combustion engine includes direct injection engines and engines with a pre-chamber.
  • the knocking prediction method according to this embodiment relates to a design method for a combustion device equipped with a premixed combustion field (combustion chamber), a power generating device (prime mover), a vehicle, etc. Furthermore, the knocking prediction method according to this embodiment also relates to a design method for a fuel that burns in a premixed combustion field.
  • the knocking prediction method is based on the premise that when a time-space transformation of a flame is applied to a premixed combustion field, the temporal changes in normalized fuel mass fraction and normalized temperature in zero-dimensional homogeneous ignition are equivalent to the spatial changes in a one-dimensional laminar premixed flame.
  • this knocking prediction method is based on the theoretical basis of the equivalence between the temporal changes in normalized fuel mass fraction and normalized temperature in zero-dimensional homogeneous ignition and the spatial changes in one-dimensional laminar premixed flames.
  • the change in flame propagation over time is about the same magnitude as the compressed premixed mixture or shock waves in a premixed combustion field. It is also clear that the Lewis number of the fuel in a premixed combustion field is a physical quantity that determines the existence of a premixed flame structure.
  • FIG. 3 is a table listing the various physical quantities used in the following theoretical explanation.
  • Equation (1) and (2) are the governing equations for zero-dimensional homogeneous ignition using normalized fuel mass fraction and temperature. These governing equations (1) and (2) are the conservation laws of energy and chemical species under constant pressure and constant enthalpy, and are given by the mass fraction of the kth species, the chemical production rate of the kth species, the molecular weight of the kth species, the mass density, the temperature, the total number of chemical species, the enthalpy of the kth species, the average specific heat, and time.
  • the mass density is given by the gas equation of state used for the interconversion of pressure, density, and temperature, as shown in the following equation (3). Note that R in this equation of state (3) is the universal gas constant.
  • equation (5) the multivariate function g shown in equation (5) below is obtained.
  • the total differential of this multivariate function g is expressed as equation (6) below, with the variables being a, mass density, mass fraction, and time.
  • the total differential of the multivariate function g is then expressed as equation (7) below based on equation (6).
  • equation (8) is obtained based on a comparison between equation (6) and equation (7). Furthermore, since equation (5) above can be defined as shown in equation (9), the relationship shown in equation (10) holds between multivariate function f and multivariate function g.
  • a candidate for the multivariate function f in equations (8) and (10) is the normalized temperature.
  • a candidate for the multivariate function g can be the normalized fuel mass fraction in the case of the single-step chemical reaction model in equation (8), or the normalized progress variable in the case of the multi-step chemical reaction model in equation (10).
  • the above equation (8) is expressed as the following equation (11), which shows the relationship between the time change of the normalized temperature and the time change of the normalized fuel mass fraction.
  • the normalized temperature is defined by the following equation (12), which has the initial temperature and final temperature as parameters.
  • the normalized fuel mass fraction is defined by the following equation (13), which has the initial fuel mass fraction and final fuel mass fraction as parameters.
  • the normalized temperature, normalized fuel mass fraction, and normalized progress variable must be convex functions, and the total derivative of the normalized temperature must be greater than 0, the total derivative of the normalized fuel mass fraction must be less than 0, and the total derivative of the normalized progress variable must be greater than 0.
  • the normalized fuel mass fraction is taken as the multivariate function g, and the signs of the normalized progress variables are simply inverted. Once either the normalized temperature or the normalized fuel mass fraction is determined, all remaining variables are determined. Based on this premise, equations (1) and (2) can be rewritten as the following equations (15) and (16) using the normalized temperature and normalized fuel mass fraction.
  • equation (11) the relationship between normalized temperature and normalized fuel mass fraction is independent of time. Therefore, equation (11) can be rewritten using residence time as follows. Note that residence time is defined by equation (18) below, which consists of the velocity of the fluid parcel and the position of the fluid parcel, and is the total time that the fluid parcel spends within the control volume. The total derivative of such residence time is expressed by equation (19) below.
  • equations (15) and (16) are transformed using equation (19), the following equations (20) and (21) are obtained. Furthermore, equation (17) can be rewritten as the following equation (22) using the parcel's position. Furthermore, the relationship shown in equation (23) holds based on equation (22).
  • equations (24) to (26) consisting of the laminar burning velocity, the mixing diffusion coefficient of the kth species, and the thermal conductivity of the premixed gas.
  • equations (25) and (26) can be rewritten as the following equations (27) and (28) using the normalized temperature and normalized fuel mass fraction.
  • equation (29) By substituting equations (22) and (23) into equation (27), we obtain the following equation (29). As can be seen by comparing equation (29) with equation (28), the only difference between the two equations is the first term on the right-hand side.
  • the Lewis number of the fuel is defined by the following equation (30). Comparing equations (28) and (29) with a focus on this Lewis number, it can be seen that both equations are equivalent when the Lewis number is 1. In other words, for a fuel with a Lewis number of 1, the temporal changes in the normalized fuel mass fraction and normalized temperature in zero-dimensional homogeneous ignition are equivalent to the spatial changes in a one-dimensional laminar premixed flame.
  • Figure 4 is a characteristic diagram showing the relationship between zero-dimensional homogeneous ignition and one-dimensional laminar premixed flame in a premixed combustion field when hydrogen, methane, propane and SNIa are used as fuel.
  • SNIa represents the carbon state in the nuclear combustion process of Type Ia supernova and has a very high Lewis number.
  • Hydrogen, methane and propane have Lewis numbers of 0.36, 0-96 and 1.95.
  • zero-dimensional homogeneous ignition as indicated by the normalized fuel mass fraction and normalized temperature, is equivalent to a one-dimensional laminar premixed flame, as indicated by the normalized fuel mass fraction and normalized temperature.
  • the profile of a one-dimensional laminar premixed flame is downwardly convex compared to the profile of a fuel with a Lewis number of 1 (zero-dimensional homogeneous ignition).
  • the profile is convex above the profile for fuels with Lewis numbers of 1 (zero-dimensional homogeneous ignition) and may intersect with the profile for fuels with Lewis numbers of 1 (zero-dimensional homogeneous ignition) in the high temperature regions (reaction zones) of the flame where heat release is high.
  • the consumption of the normalized fuel mass fraction at a given normalized temperature can be considered as the progress of the combustion reaction.
  • the progress of the combustion reaction is in the following order: fuel with a Lewis number less than 1, fuel with a Lewis number of 1 (zero-dimensional homogeneous ignition), and fuel with a Lewis number greater than 1.
  • gasoline one of the fuels, has a Lewis number greater than 1, so it is expected that self-ignition (knocking) will occur whenever the unburned gas region rises to a temperature and pressure that makes flame propagation impossible.
  • the characteristic time of the flame behind the shock wave is similar to that at the time of ignition, which can affect the transition from flame propagation to detonation.
  • gasoline since gasoline has a Lewis number greater than 1, there may be an area behind the shock wave where flame propagation is not possible.
  • the normalized temperature and normalized fuel mass fraction must be convex functions in order to maintain the Legendre transformation, so the calculation domain must be as short as possible. Therefore, in this simulation, the size of the domain containing the unburned mixture, the flame reaction zone, and the burned gas in the premixed combustion field is set to the smallest limit that can contain the flame.
  • the numerical conditions (initial conditions) for the pressure and equivalence ratio are fixed values, for example, 0.1 MPa and 1.0.
  • Such a simulation corresponds to the knocking prediction method according to this embodiment. Furthermore, this simulation is the processing operation of the knocking prediction device according to this embodiment, that is, information processing executed by the calculation unit 2 based on the knocking prediction program pre-stored in the storage unit 1.
  • the operator first operates the operation unit 3 to input the initial conditions required for the simulation, such as temperature, pressure, and the type of fuel, such as hydrogen fuel or n-heptane fuel, and the equivalence ratio, into the calculation unit 2 (step S1).
  • the initial conditions required for the simulation such as temperature, pressure, and the type of fuel, such as hydrogen fuel or n-heptane fuel, and the equivalence ratio
  • Calculation unit 2 uses the temperature, pressure, fuel type, and equivalence ratio input in step S1 and gradually increases the inlet temperature (calculation start temperature) to perform a one-dimensional steady-state analysis (calculation process) to calculate the flame propagation limit temperature for each fuel (step S2).
  • This one-dimensional steady-state analysis is a convergence calculation that shrinks the flame front and flame rear until the one-dimensional premixed flame structure becomes the narrowest calculation region that results in a unique solution.
  • the calculation unit 2 sequentially calculates the burning velocity when the inlet temperature is gradually increased using a one-dimensional steady-state analysis of the premixed combustion field, but the burning velocity calculation may converge or diverge depending on the set value of the inlet temperature.
  • the calculation unit 2 performs a convergence judgment (judgment process) for each calculation of the burning rate for a certain inlet temperature (calculation start temperature) (step S3). That is, in step S3, the calculation unit 2 judges whether or not the burning rate cannot be obtained because a flame can no longer exist in the premixed combustion field. Then, if the judgment result in step S3 is "convergence," the calculation unit 2 sets the inlet temperature even higher and executes the burning rate calculation.
  • step S4 occurrence condition acquisition process.
  • the calculation unit 2 acquires the flame propagation limit temperature (inlet temperature) corresponding to the last acquired combustion speed as the occurrence condition for knocking.
  • the calculation unit 2 then outputs the simulation result including the flame propagation limit temperature acquired in this manner as the occurrence condition for knocking to the output unit 4.
  • Figure 5 shows an example of this simulation result (one-dimensional steady-state analysis result), that is, a calculation result including the knocking occurrence conditions that is output by the calculation unit 2 to the output unit 4.
  • Figure 6 shows an example of the relationship between the inlet temperature and the combustion speed of the premixed combustion field based on such a simulation result.
  • the black circles are the simulation results in this embodiment (one-dimensional steady-state analysis result), and the squares are the results of direct numerical simulation (DNS) shown as a comparative example.
  • DMS direct numerical simulation
  • Figure 6 shows that for hydrogen fuel, the burning rate increases with increasing inlet or initial temperature, and a propagating flame exists even at an inlet temperature of 3000 K.
  • Figure 6 also shows that for n-heptane fuel, the burning rate increases with increasing inlet temperature, but a propagating flame cannot exist above 1270 K.
  • the theoretical basis for this embodiment is that it theoretically shows that the relationship between normalized fuel mass fraction and normalized temperature is equivalent for zero-dimensional homogeneous ignition and one-dimensional laminar premixed flames, whose Lewis number after space-time transformation is 1, and that the Lewis number is a physical quantity that determines the existence of a premixed flame structure near the explosive transition of the flame in a premixed combustion field.
  • Figure 7 shows a comparative example of this embodiment.
  • this comparative example is the result of calculating the "Livengood-Wu integral" using direct numerical simulation (DNS).
  • DNS direct numerical simulation
  • the "Livengood-Wu integral” causes auto-ignition (knocking) when the time integral value becomes 1.
  • this embodiment is based on the knowledge that for a fuel with a Lewis number of 1, the temporal changes in normalized fuel mass fraction and normalized temperature in zero-dimensional homogeneous ignition are equivalent to the spatial changes in a one-dimensional laminar premixed flame, and therefore has higher accuracy than the "Livengood-Wu integral.”
  • this embodiment can provide a knocking prediction method and knocking prediction device that can further improve accuracy.
  • knocking prediction method and knocking prediction device have high knocking prediction accuracy, which provides great benefits to the design of combustion devices, power generation devices (internal combustion engines), vehicles, fuels, etc.
  • the knocking prediction method and knocking prediction device make it possible to estimate the margin for knocking that should be considered in the design of combustion devices, power generating devices, vehicles, etc. more narrowly than ever before.
  • the knocking prediction method and knocking prediction device can provide a method for designing a combustion device, a method for designing a power generation device, and a method for designing a vehicle that can improve the performance of a combustion device, a power generation device, a vehicle, etc.
  • the knocking prediction method and knocking prediction device can predict knocking with high accuracy for each type of fuel, which can contribute to the design of fuels that avoid knocking. Therefore, the knocking prediction method and knocking prediction device according to this embodiment can provide a fuel setting method that can improve the performance of the fuel.
  • this embodiment predicts knocking of the fuel to be evaluated based on a one-dimensional steady-state analysis (simulation) of the fuel to be evaluated. Therefore, this disclosure is not limited to pure-substance fuels, but can be applied to multi-substance fuels that are mixtures of multiple components, as long as one-dimensional steady-state analysis (simulation) is possible.
  • fuels for which chemical species and elementary reactions have been identified include surrogate fuels consisting of up to five components defined in the SIP (Innovative Combustion Technology) surrogate chemical reaction mechanism.
  • SIP Innovative Combustion Technology
  • chemical species and elementary reactions have not yet been identified for multi-substance fuels consisting of six or more components.
  • the engine control device is an ECU (Electronic Control Unit) 11 shown in FIG. 8.
  • the vehicle according to this embodiment is an automobile M equipped with a vehicle operation device 10, an ECU 11, an engine 12, and a communication device 13, as also shown in FIG. 8.
  • the vehicle operation device 10 is an operation panel of the automobile M, which generates an operation instruction signal based on the input operation of the driver and outputs it to the ECU 11.
  • the operation instruction signal is a signal that indicates the driver's operation instruction for the automobile M.
  • the ECU 11 is a so-called embedded computer, and includes an arithmetic circuit, a memory circuit, and an input/output circuit.
  • the ECU 11 controls the engine 12 by executing an engine control program that is pre-stored in the memory circuit.
  • the ECU 11 generates various engine control signals by executing an engine control program in an arithmetic circuit.
  • These engine control signals are signals that indicate various operation amounts (engine operation amounts) for the engine 12, such as the opening degree of the air intake valve and the opening degree and timing of the fuel injection valve.
  • the ECU 11 controls the engine 12 by outputting such engine control signals to the engine 12 using an input/output circuit.
  • the engine 12 is an internal combustion engine that generates power by burning fuel in a combustion chamber, and is equipped with a sensor group 12a consisting of multiple sensors.
  • This sensor group 12a detects various state quantities related to the engine 12 as engine control quantities, such as the opening degree of the air intake valve, the operation of the fuel injection valve, and the carbon dioxide contained in the exhaust gas, and outputs these engine control quantities to the ECU 11 as sensor signals.
  • the ECU 11 When generating the engine control signal, the ECU 11 refers to the operation instruction signal input from the vehicle operation device 10 and various sensor signals input from the sensor group 12a. The ECU 11 generates the engine control signal (engine operation amount) by processing the operation instruction signal and engine control amount according to an engine control program.
  • the ECU 11 in this embodiment also includes the above-mentioned knocking prediction device A, as shown in FIG. 8.
  • the knocking prediction device A is a computer that predicts the occurrence of knocking by executing, in the calculation unit 2, a knocking prediction program that is pre-stored in the storage unit 1, as described above.
  • knocking prediction device A in ECU 11 derives the fuel composition from big data based on information already collected from surrounding vehicles using date, time and location information, and obtains the knocking occurrence prediction result by executing a knocking prediction program on that fuel composition.
  • the timing for obtaining the knocking occurrence prediction result is immediately after fuel is refueled into automobile M and before automobile M starts to drive.
  • the fuel composition is constantly updated based on the driving conditions while automobile M is driving, and when there is a large deviation from the initially predicted fuel composition, the knocking occurrence prediction result is also updated.
  • the updated information is transmitted by data communication such as wirelessly, and is stored as big data that is useful for surrounding vehicles and for development and improvement.
  • the ECU 11 is an embedded computer that includes the knocking prediction device A as an additional function associated with the engine control function by pre-storing a knocking prediction program in addition to the engine control program.
  • Such an ECU 11 refers to the knocking prediction results obtained by executing the knocking prediction program when generating an engine control signal. For example, when the occurrence of knocking is predicted by executing the knocking prediction program, the ECU 11 generates an engine control signal to avoid knocking.
  • the ECU 11 provides the knocking prediction results based on the knocking prediction program, the engine operation amount indicated by the engine control signal, and the engine control amount obtained by the sensor group 12a, together with additional information, to an external data collection device D.
  • the additional information is engine-related information other than the engine control amount and the engine operation amount, such as the date, time, outside temperature, season, fuel composition, and the traveling position of the automobile M.
  • the ECU 11 sequentially stores in the memory circuit as learning data (time-series data) the knocking prediction results obtained before the vehicle M starts to move, and the engine operation amounts, engine control amounts, and associated information obtained in chronological order while the vehicle M is moving.
  • the ECU 11 reads out such learning data from the memory circuit at a predetermined timing and outputs it to the communication device 13.
  • the communication device 13 is a wireless communication device that is provided between the ECU 11 and the data collection device D and mediates communication between the ECU 11 and the data collection device D. In other words, the communication device 13 converts the learning data input from the ECU 11 into a wireless signal and transmits it to the data collection device D.
  • such learning data is provided to the data collection device D from multiple automobiles M1 to Mn (n: a natural number of 2 or more) equipped with an ECU 11 incorporating a knocking prediction device A.
  • the data collection device D sequentially accumulates learning data related to multiple automobiles M1 to Mn as big data.
  • This kind of big data can be used in a variety of ways, as each learning data item contains knocking prediction results corresponding to engine operation amounts, engine control amounts, and additional information.
  • knocking is avoided by switching to the motor under conditions that would cause knocking, or by increasing motor assist, and the big data in this embodiment can be useful data for considering designs to avoid such knocking.
  • the present invention makes it possible to omit this data collection and input process during automobile manufacturing. Furthermore, if an automobile is equipped with the functions of the present invention, the knocking suppression effect can be achieved without going through the process of inputting data in advance, even for new fuel compositions such as e-fuel that will be released after the sale. Furthermore, when a sub-chamber is provided in the combustion chamber of the engine 12, it is possible to flexibly control ignition properties, thereby improving thermal efficiency.
  • the present disclosure can be used in knocking prediction methods, knocking prediction devices, engine control devices, vehicles, combustion device design methods, power generation device design methods, vehicle design methods, and fuel design methods.
  • a Knocking prediction device D Data collection device M Automobile (vehicle) REFERENCE SIGNS LIST 1 Memory unit 2 Calculation unit 3 Operation unit 4 Output unit 10 Vehicle operation device 11 ECU (engine control unit) 12 engine 12a sensor group 13 communication device

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本開示は、精度をより向上させることが可能なノッキング予測方法の提供を目的とする。　本開示は、予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、予混合燃焼場において、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化とが等価であるという知見に基づいてノッキングを予測する、という手段を採用する。

Description

ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法
　本開示は、ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法に関する。
　周知のようにガソリン・エンジン等、予混合燃焼方式の内燃機関では、ノッキングの抑制が設計上の重要課題である。ノッキングは、燃焼室内の未燃ガスが自着火することによって発生する衝撃波を伴う現象であり、エンジンを破損する可能性がある。内燃機関の設計では、燃焼室の圧力及び温度(圧縮比)をノッキングが発生しないように下げることが行われる。
　一方、ノッキング予測方法として、「Livengood-Wu積分」が周知である。この「Livengood-Wu積分」は、燃焼室内の圧力及び温度を変化させた際の着火遅れ時間を見積もって時間積分していく手法であり、時間積分値が「１」になった時に自着火つまりノッキングが発生するというものである。「Livengood-Wu積分」は、下記特許文献１、２等に示すように実機エンジンの設計において広く利用されている
日本国特開２００４－３３２５８４号公報 日本国特開２００７－１２７００４号公報
　ところで、「Livengood-Wu積分」は簡便な手法であり、ノッキング（自着火）の発生をある程度予測できるものの、その予測精度は必ずしも十分なものとは言えない。動力発生装置として機能する内燃機関の能力を最大限利用するためには、「Livengood-Wu積分」よりも予測精度が高いノッキング予測方法の確立が切望されている。
　本開示は、上述した事情に鑑みてなされたものであり、精度をより向上させることが可能なノッキング予測方法及びノッキング予測装置、エンジン制御装置、車両の提供を目的とする。
　また、本開示は、燃焼装置、動力発生装置、車両の設計方法及び燃料の実質的な性能を高めることが可能な燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法の提供を目的とする。
　本開示の第１の態様に係るノッキング予測方法は、予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、前記予混合燃焼場において、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化とが等価であるという知見に基づいてノッキングを予測するものである。
　本開示の第２の態様に係るノッキング予測方法は、予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、正規化温度と燃料の正規化燃料質量分率との関係を示す特性図において、ルイス数が１より大きい燃料を前記ノッキングの可能性があるとし、前記ルイス数が１より小さい燃料を前記ノッキングの可能性がないとするものである。
　本開示の第３の態様に係るノッキング予測方法は、予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、前記予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算する計算工程と、前記予混合燃焼場に火炎が存在できなくなるために前記燃焼速度が取得不能になったか否かを判定する判定工程と、最後に取得された前記燃焼速度に対応する前記入口温度を前記ノッキングの発生条件として取得する発生条件取得工程とを有するものである。
　本開示の第４の態様に係るノッキング予測方法は、第３の態様に係るノッキング予測方法において、前記計算工程では、前記予混合燃焼場における未燃混合気、火炎の反応帯及び既燃焼気体を含む領域の大きさを、火炎が包含できる限界まで小さくして計算するものである。
　本開示の第５の態様に係るノッキング予測方法は、第３又は第４の態様に係るノッキング予測方法において、ルイス数が１より大きい燃料を予測対象とするものである。
　本開示の第１の態様に係るノッキング予測装置は、予混合燃焼場で発生し得るノッキングを予測するノッキング予測装置であって、前記予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算する計算手段と、前記予混合燃焼場に火炎が存在できなくなるために前記燃焼速度が取得不能になったか否かを判定する判定手段と、最後に取得された前記燃焼速度に対応する前記入口温度を前記ノッキングの発生条件として取得する発生条件取得手段とを備えるものである。
　本開示の第２の態様に係るノッキング予測装置は、第１の態様に係るノッキング予測装置において、前記計算手段は、前記予混合燃焼場における未燃混合気、火炎の反応帯及び既燃焼気体を含む領域の大きさを、火炎が包含できる限界まで小さくして計算するものである。
　本開示の第３の態様に係るノッキング予測装置は、第１又は第２の態様に係るノッキング予測装置において、前記計算手段は、ルイス数が１より大きい燃料を計算対象とするものである。
　本開示の第１の態様に係るエンジン制御装置は、第１の態様又は第２の態様に係るノッキング予測装置を備え、前記ノッキング予測装置のノッキング予測結果に基づいてエンジン制御信号を生成するものである。
　本開示の第２の態様に係るエンジン制御装置は、第１の態様に係るエンジン制御装置において、外部との通信が自在な通信装置をさらに備え、少なくとも前記ノッキング予測結果及び前記エンジン制御信号が示すエンジン操作量を前記通信装置を用いて外部に送信するものである。
　本開示の第３の態様に係るエンジン制御装置は、第２の態様に係るエンジン制御装置において、エンジン制御量に加え、該エンジン制御量及び前記エンジン操作量以外のエンジン関連情報を前記外部に送信するものである。
　本開示の第１の態様に係る車両は、第１の態様又は第２の態様に係るエンジン制御装置と、該エンジン制御装置によって制御されるエンジンと、を備えるものである。
　本開示の第１の態様に係る燃焼装置の設計方法は、第３又は第４の態様に係るノッキング予測方法に基づいて推定されたノッキングの発生条件に基づいて燃焼装置を設計するというものである。
　本開示の第１の態様に係る動力発生装置の設計方法は、第３又は第４の態様に係るノッキング予測方法に基づいて推定されたノッキングの発生条件に基づいて動力発生装置を設計するというものである。
　本開示の第１の態様に係る車両の設計方法は、第３又は第４の態様に係るノッキング予測方法に基づいて推定されたノッキングの発生条件に基づいて車両を設計するというものである。
　本開示の第１の態様に係る燃料の設計方法は、第３又は第４の態様に係るノッキング予測方法に基づいて推定されたノッキングの発生条件に基づいて燃料を設計するというものである。
　本開示によれば、精度をより向上させることが可能なノッキング予測方法、ノッキング予測装置、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法を提供することが可能である。
　また、本開示によれば、燃焼装置、動力発生装置、車両の設計方法及び燃料の実質的な性能を高めることが可能な燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法を提供することが可能である。
本開示の実施形態に係るノッキング予測方法を示すフローチャートである。 本開示の実施形態に係るノッキング予測装置の機能構成を示すブロック図である。 本開示の実施形態の理論説明で取り扱う各種物理量の一覧表である。 本開示の実施形態におけるノッキング予測特性を示す第１の特性図である。 本開示の実施形態におけるノッキング予測特性を示す第２の特性図である。 本開示の実施形態におけるノッキング予測特性を示す第３の特性図である。 本開示の実施形態におけるノッキング予測方法の比較例を示す特性図である。 本開示の実施形態に係るエンジン制御装置及び車両を示すブロック図である。 本開示の実施形態に係る車両の外部通信を示す模式図である。
　以下、図面を参照して、本開示の実施形態について説明する。
　図１は、本実施形態に係るノッキング予測方法を示すフローチャートである。このフローチャートは、燃料と空気（酸化剤）とからなる予混合気を火花点火のような強制的な手段により着火させることにより燃焼させる予混合燃焼場（燃焼室）について、ノッキングの発生条件である火炎伝播限界温度つまりノッキングが発生する初期温度の取得手順を示している。
　また、このフローチャートは、所定のコンピュータによって実行されるノッキング予測プログラム（シミュレーションプログラム）の処理手順の概要を示すものである。すなわち、本実施形態に係るノッキング予測方法は、コンピュータを用いた情報処理（シミュレーション）によって実現されるものである。
　本実施形態に係るノッキング予測方法を実行するコンピュータつまり上記ノッキング予測プログラムが搭載されたコンピュータは、本実施形態に係るノッキング予測装置である。すなわち、本実施形態に係るノッキング予測装置は、ノッキング予測プログラムを実行することにより、予混合燃焼場におけるノッキングの発生条件を出力する装置である。
　図２は、本実施形態に係るノッキング予測装置Ａの機能構成を示すブロック図である。ノッキング予測装置Ａは、図２に示すように、ノッキング予測プログラムを予め記憶する記憶部１、ノッキング予測プログラムを実行する演算部２、演算部２に作業者の操作指示を入力する操作部３、演算部２の演算結果つまり上記火炎伝播限界温度（初期温度）を外部に出力する出力部４を少なくとも備えている。
　ノッキング予測装置Ａの構成要素のうち、記憶部１、演算部２及び操作部３は、本開示における計算手段を構成している。詳細については後述するが、記憶部１、演算部２及び操作部３は、協働することにより、予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算する。
　また、ノッキング予測装置Ａにおける記憶部１、演算部２及び操作部３は、本開示における判定手段を構成している。すなわち、記憶部１、演算部２及び操作部３は、協働することにより、予混合燃焼場に火炎が存在できなくなるために燃焼速度が取得不能になったか否かを判定する。
　さらに、ノッキング予測装置Ａにおける記憶部１、演算部２及び操作部３は、本開示における発生条件取得手段を構成している。すなわち、記憶部１、演算部２及び操作部３は、協働することにより、最後に取得された燃焼速度に対応する入口温度（火炎伝播限界温度）をノッキングの発生条件として取得する。
　ここで、本実施形態に係るノッキング予測装置Ａの計算手段は、本実施形態に係るノッキング予測方法における計算工程に対応する。また、本実施形態に係るノッキング予測装置Ａの判定手段は、本実施形態に係るノッキング予測方法における判定工程に対応する。さらに、本実施形態に係るノッキング予測装置Ａの生条件取得手段は、本実施形態に係るノッキング予測方法における発生条件取得工程に対応する。
　本実施形態における予混合燃焼場は、例えば車両に動力発生装置として搭載されるガソリン・エンジン等の内燃機関（原動機）の燃焼室に相当する燃焼場である。ガソリン・エンジン等の設計においては、ノッキングが発生しないように、つまり未燃ガスの自着火が発生しないように圧縮比を下げる等の対策が施される。なお、上記内燃機関には、直噴エンジンや副室を有するエンジンが含まれる。
　すなわち、本実施形態に係るノッキング予測方法は、予混合燃焼場（燃焼室）を備える燃焼装置、動力発生装置（原動機）及び車両等の設計方法に関係するものである。また、本実施形態に係るノッキング予測方法は、予混合燃焼場で燃焼する燃料の設計方法にも関係するものである。
　このような本実施形態に係るノッキング予測方法は、予混合燃焼場に火炎の時空間変換が適用される場合において、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化が１次元層流予混合火炎における空間的変化と等価であることを前提としている。
　すなわち、このノッキング予測方法は、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化との等価性に関する理論的根拠に基づくものである。
　この理論的根拠によれば、予混合燃焼場の予熱域における正規化燃料質量分率の低下の度合いは、燃料のルイス数が１未満、ルイス数が１(＝点火)、ルイス数が１より大の順で速くなる。さらに、爆発的遷移付近の予熱域における自着火は、１より大きいルイス数を有する１次元層流予混合火炎に対してのみ起こり得ることが導き出される。
　すなわち、火炎伝播の時間変化は、予混合燃焼場において圧縮された予混合気又は衝撃波と同程度の大きさである。また、予混合燃焼場における燃料のルイス数は、予混合火炎構造の存在を決定づける物理量であることがわかる。
　以下では、本実施形態に係るノッキング予測方法の理論的根拠について説明する。なお、図３は、以下の理論説明で取り扱う各種物理量の一覧表である。
　下式（１）、（２）は、正規化燃料質量分率及び温度を用いたゼロ次元均質点火の支配方程式である。この支配方程式（１）、（２）は、一定圧力及び一定エンタルピー下におけるエネルギー及び化学種の保存則であり、ｋ番目の種の質量分率、k番目の種の化学生成速度、k番目の種の分子量、質量密度、温度、化学種の総数、k番目の種のエンタルピー、平均比熱及び時間によって与えられる。
Figure JPOXMLDOC01-appb-M000001
　また、質量密度は、下式（３）に示すように、圧力、密度及び温度の相互変換に使用される気体の状態方程式として与えられる。なお、この状態方程式（３）におけるＲは普遍的なガス定数である。
Figure JPOXMLDOC01-appb-M000002
　ここで、多変量関数ｆからのLegendre変換を考える。温度、質量密度、質量分率及び時間を変数としてガス混合物の特性を示す多変量関数ｆの全微分は下式（４）で与えられる。
Figure JPOXMLDOC01-appb-M000003
　この式（４）において多変量関数ｆの温度に対する偏微分を「ａ」と置くと、下式（５）に示す多変量関数ｇが得られる。また、この多変量関数ｇの全微分は、上記ａ、質量密度、質量分率、時間を変数とする下式（６）のように表される。そして、多変量関数ｇの全微分は、式（６）に基づいて下式（７）のように表される。
Figure JPOXMLDOC01-appb-M000004
　さらに、このような式（６）と式（７）との比較に基づいて式（８）が得られる。また、上記式（５）については式（９）に示すように定義することができるので、多変量関数ｆと多変量関数ｇとの間には式（１０）に示す関係が成立する。
Figure JPOXMLDOC01-appb-M000005
　ここで、式（８）及び式（１０）における多変量関数ｆの候補は、正規化温度である。また、多変量関数gの候補は、式（８）のシングルステップ化学反応モデルの場合における正規化燃料質量分率、又は式（１０）のマルチステップ化学反応モデルの場合における正規化進捗変数が考えられる。
　ある時間における多変量関数f又は多変量関数gが決定されると、予混合燃焼場における予混合気の特性を表す変数（時間、質量密度及び燃料質量分率）が全て決定可能であると仮定すると、すなわち多変量関数f及び多変量関数gが時間方向にのみ依存すると仮定すると、上式（８）、（１０）は偏微分式から常微分式に書き換えることができる。
　この結果、正規化温度が多変量関数fとして適用され、また正規化燃料質量分率が多変量関数gとして適用される場合、上式（８）は正規化温度の時間変化と正規化燃料質量分率の時間変化との関係を示す下式（１１）のように表される。なお、正規化温度は、初期温度と最終温度とをパラメータとする下式（１２）によって定義される。正規化燃料質量分率は、初期燃料質量分率と最終燃料質量分率とをパラメータとする下式（１３）によって定義される。
Figure JPOXMLDOC01-appb-M000006
　また、正規化温度が多変量関数fとして適用され、多段階化学反応モデルの正規化進行変数が多変量関数gとして適用される場合には、上式（１０）は正規化温度の時間変化と正規化進行変数の時間変化との関係を示す下式（１４）のように表される。
Figure JPOXMLDOC01-appb-M000007
　ここで、Legendre変換を可能にするためには、正規化温度、正規化燃料質量分率及び正規化進行変数は凸関数でなければならず、また正規化温度の全微分は０より大きな値、正規化燃料質量分率の全微分は０より小さな値、また正規化進行変数の全微分は０より大きな値でなければならない。この制約条件は、常にシングルステップ化学反応モデルに当てはまり、通常ではマルチステップ化学反応モデルにも当てはまる。
　以下の説明では、説明を簡潔にするために正規化燃料質量分率を多変量関数g とし、正規化進行変数については符号を単に反転させる。正規化温度及び正規化燃料質量分率は、いずれかが決定されると残りの全ての変数が決定される。この前提に基づくと、式（１）、（２）は、正規化温度及び正規化燃料質量分率を用いた下式（１５）、（１６）のように書換えられる。
Figure JPOXMLDOC01-appb-M000008
　上述した式（１１）に示されるように正規化温度と正規化燃料質量分率との関係は時間とは無関係である。したがって、式（１１）は、滞留時間を用いることにより以下のように書き換えることができる。なお、滞留時間は、流体パーセルの速度と流体パーセルの位置からなる下式（１８）によって定義されるものであり、流体小包が制御容積内で費やした総時間である。このような滞留時間の全微分は、下式（１９）によって表される。
Figure JPOXMLDOC01-appb-M000009
　この式（１９）を用いて式（１５）、（１６）を変形すると、下式（２０）、（２１）が得られる。また、式（１７）は、パーセルの位置を用いることにより下式（２２）のように書き換えることができる。さらに、この式（２２）に基づいて式（２３）に示す関係が成立する。
Figure JPOXMLDOC01-appb-M000010
　一方、１次元層流予混合火炎の保存方程式は、層流燃焼速度、ｋ番目の種の混合拡散係数及び予混合気の熱伝導率からなる下式（２４）～（２６）によって与えられる。また、式（２５）、（２６）については、正規化温度及び正規化燃料質量分率を用いることにより下式（２７）、（２８）のように書き換えることができる。
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
　式（２２）、（２３）を式（２７）に代入することにより下式（２９）が得られる。この式（２９）を式（２８）と比較するとわかるように、両式の違いは右辺の第1項のみである。
Figure JPOXMLDOC01-appb-M000013
　ここで、燃料のルイス数は下式（３０）によって定義される。このルイス数に着目して式（２８）、（２９）を対比すると、両式はルイス数が１のとき等しくなることがわかる。すなわち、ルイス数が１の燃料について、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化とは等価である。
Figure JPOXMLDOC01-appb-M000014
　図４は、水素、メタン、プロパン及びＳＮＩａを燃料とした場合の予混合燃焼場におけるゼロ次元均質点火と1次元層流予混合火炎との関係を示す特性図である。ＳＮＩａは、Ia型超新星の核燃焼過程における炭素状態を示すものであり、非常に高いルイス数を有する。水素、メタン及びプロパンは、０．３６、０－９６及び１．９５のルイス数を有する。
　上述したように、ルイス数が１の燃料では、正規化燃料質量分率及び正規化温度によって示されるゼロ次元均質点火は、正規化燃料質量分率及び正規化温度によって示される１次元層流予混合火炎に等しい。また、ルイス数が１より小さい燃料は、１次元層流予混合火炎のプロファイルは、ルイス数が１の燃料(ゼロ次元均質点火)のプロファイルよりも下に凸状である。
　一方、ルイス数が１より大きい燃料では、プロファイルは、ルイス数が１の燃料(ゼロ次元均質点火)のプロファイルより上に凸状であり、熱放出が高い火炎の高温領域(反応域)において、ルイス数が１の燃料(ゼロ次元均質点火)のプロファイルと交差する可能性がある。
　これは、ルイス数が１より小さい燃料では、燃料の物質拡散係数が熱拡散係数よりも大きいために、燃料は加熱されるよりも優先的に拡散されることを示している。これに対して、ルイス数が１より大きい燃料では、燃料の熱拡散係数が物質拡散係数よりも大きいために、燃料は拡散されるよりも優先的に加熱される。
　すなわち、ルイス数が１より大きい燃料において拡散されるよりも優先的に加熱されるという事項は、周知の技術的知見に合致する。温度と燃料質量分率との間の関係に対するルイス数の影響は、周知のBuckmasterの書籍において理論的に論じられており、図４の特性はその理論と一致している。
　以下では、予混合燃焼場についてルイス数の影響が明らかな比較的低温の領域(予熱域)について検討する。なお、この検討では周知のソレット効果及びデュファー効果は無視できると仮定する。
　図４の特性図において、所定の正規化温度における正規化燃料質量分率の消費は、燃焼反応の進行度合いとみなすことができる。すなわち、燃焼反応の進行度合いは、ルイス数が１より小さい燃料、ルイス数が１の燃料(ゼロ次元均質点火)、ルイス数が１より大きい燃料の順となる。
　ルイス数が１より小さい燃料では、予熱域における自己着火（ノッキング）は発生しない。すなわち、予混合燃焼場における初期温度及び初期圧力が非常に高い場合であっても、1次元層流予混合火炎が存在し得る。これに対して、ルイス数が１より大きい燃料では、予熱域において自己着火（ノッキング）が発生する可能性がある。すなわち、予混合燃焼場における初期温度及び初期圧力によっては、１次元層流予混合火炎が存在し得ない。
　例えば、燃料の1つであるガソリンは、ルイス数が１より大きいので、未燃ガス領域が火炎伝播を不可能にする温度及び圧力まで上昇するたびに自己着火（ノッキング）が発生することが予想される。すなわち、衝撃波の背後の火炎の特性時間が点火時と同程度の大きさであるため、火炎伝播から爆轟への遷移に影響を及ぼし得る。また、ガソリンでは、ルイス数が１を超えているので、火炎伝播ができない衝撃波の背後の領域が存在する可能性がある。
　ゼロ次元均質点火及び1次元層流予混合火炎について、定圧及び定エンタルピーを条件として２種類の予混合気、つまりルイス数が１よりも小さい水素と空気との予混合気（水素燃料）及びルイス数が１よりも大きいｎ－ヘプタンと空気との予混合気（ｎ－ヘプタン燃料）についてシミュレーション（１次元定常解析）を行った。
　このシミュレーションでは、周知の多段階化学反応モデルを用いることによりゼロ次元均質点火及び1次元層流予混合火炎の１次元定常解析を行った。水素燃料の多段階化学反応モデルは、周知のUT-JAXAモデルであり、ｎ－ヘプタン燃料の多段階化学反応モデルは、周知のReduced（簡略化された）SIPモデルである。
　このシミュレーションでは、Legendre変換を維持するために、正規化温度及び正規化燃料質量分率は凸関数でなければならないので、計算領域をできるだけ短くしなければない。したがって、このシミュレーションでは、予混合燃焼場における未燃混合気、火炎の反応帯及び既燃焼気体を含む領域の大きさを、火炎が包含できる限界まで小さく設定する。なお、圧力及び当量比の数値条件（初期条件）は、例えば０．１ＭＰａ及び１．０の固定値である。
　また、このシミュレーションでは、１次元層流予混合火炎の１次元定常解析を入口温度を３００～３０００Ｋまで１００Ｋ刻みで変化させることによって行った。計算領域を減少させても計算の失敗が確認された場合には、入口温度を１Ｋ刻みで変化させることによって数値計算を実行した。
　さらに、このシミュレーションでは、１次元層流予混合火炎の１次元定常解析を実施した後、１次元層流予混合火炎を計算し得る初期温度を使用して、ゼロ次元均質点火の数値計算を実施した。
　このようなシミュレーション（１次元定常解析）は、本実施形態に係るノッキング予測方法に相当するものである。また、このシミュレーションは、本実施形態に係るノッキング予測装置の処理動作つまり記憶部１に予め記憶されたノッキング予測プログラムに基づいて演算部２が実行する情報処理である。
　すなわち、本実施形態に係るノッキング予測方法では、最初に作業者が操作部３を操作することによってシミュレーションに必要な初期条件である温度、圧力、また水素燃料やｎ－ヘプタン燃料等の燃料種及び当量比を演算部２に入力する（ステップＳ１）。
　演算部２は、ステップＳ１で入力された温度、圧力、燃料種及び当量比を用いるとともに入口温度（計算開始温度）を徐々に増加させて1次元定常解析（計算工程）を実行することにより、各燃料に関する火炎伝播限界温度を計算する（ステップＳ２）。
　この1次元定常解析は、１次元予混合火炎構造を一意の解になる最も狭い計算領域になるまで火炎前方及び火炎後方を縮める収束計算である。すなわち、演算部２は、予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算するが、入口温度の設定値に応じて燃焼速度の計算が収束する場合と発散する場合とがある。
　演算部２は、ある入口温度（計算開始温度）に関する対する燃焼速度の計算毎に収束判定（判定工程）を行う（ステップＳ３）。すなわち、ステップＳ３において、演算部２は、予混合燃焼場に火炎が存在できなくなるために燃焼速度が取得不能になったか否かを判定する。そして、演算部２は、ステップＳ３の判定結果が「収束」の場合、入口温度をさらに高く設定して燃焼速度の計算を実行する。
　一方、演算部２は、ステップＳ３の判定結果が「発散」の場合には、前回の燃焼速度の計算に対応する入口温度を火炎伝播限界温度に設定する（ステップＳ４：発生条件取得工程）。すなわち、演算部２は、最後に取得された燃焼速度に対応する火炎伝播限界温度（入口温度）をノッキングの発生条件として取得する。そして、演算部２は、このようにして取得した火炎伝播限界温度をノッキングの発生条件として含むシミュレーション結果を出力部４に出力する。
　図５は、このシミュレーション結果（１次元定常解析結果）つまり演算部２が出力部４に出力するノッキングの発生条件を含む計算結果の一例である。また。図６は、このようなシミュレーション結果に基づく予混合燃焼場の入口温度と燃焼速度との関係の一例である。この図６において、黒丸は本実施形態におけるシミュレーション結果（１次元定常解析結果）であり、四角は比較例として示す直接数値計算（ＤＮＳ：Direct Numerical Simulation）の結果である。
　図６は、水素燃料の場合、燃焼速度は入口または初期温度の上昇と共に増加し、３０００Ｋの入口温度でも伝播火炎が存在することを示している。また、この図６は、ｎ－ヘプタン燃料の場合、燃焼速度は入口温度の上昇と共に増加するが、伝播火炎は１２７０Ｋを超えて存在することができないことを示している。
　すなわち、本実施形態における理論的根拠は、正規化燃料質量分率と正規化温度との関係が空間－時間変換後のルイス数が１を有するゼロ次元均質点火及び1次元層流予混合火炎について等価であることが理論的に示すものであり、またルイス数が予混合燃焼場における火炎の爆発的遷移付近の予混合火炎構造の存在を決定づける物理量であることを理論的に示すものである。
　ルイス数が１よりも小さい燃料の１次元層流予混合火炎について、予熱域において点火は起こらず、よって火炎構造は常に存在する。これに対して、ルイス数が１よりも大きい燃料の１次元層流予混合火炎については、予熱域で発火する可能性があり、温度がある閾値より高い場合には、火炎構造が存在しない。
　換言すると、点火及び火炎の特徴的な時間が衝撃波の背後またはSIエンジンの爆発的遷移と同程度の大きさである場合、点火及び火炎についてはルイス数を介して論じられるべきである。
　図７は、本実施形態の比較例である。すなわち、この比較例は、直接数値計算（ＤＮＳ）による「Livengood-Wu積分」の計算結果である。なお、直接数値計算（ＤＮＳ）が予混合燃焼場の燃焼状態を精度よくシミュレートするものであることは周知事項である。また、周知のように、「Livengood-Wu積分」は、時間積分値が１になった時に自着火（ノッキング）が発生するというものである。
　しかしながら、この計算結果によれば、時間積分値が１になった時点（ＣＦＩ）ではノッキングが発生していない。この計算結果によれば、ノッキングの発生は、時間積分値が１を超えた以降の時点（ＫＯ）である。したがって、「Livengood-Wu積分」に基づくノッキング予測方法は、精度の面で問題を含んでいる。
　これに対して、本実施形態によれば、ルイス数が１の燃料についてゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化とは等価であるという知見に基づくので、精度が「Livengood-Wu積分」より高い。すなわち、本実施形態によれば、精度をより向上させることが可能なノッキング予測方法及びノッキング予測装置を提供することが可能である。
　また、本実施形態に係るノッキング予測方法及びノッキング予測装置は、ノッキングの予測精度が高いので燃焼装置、動力発生装置（内燃機関）、車両及び燃料等の設計に対して多大な恩恵をもたらすものである。
　例えば燃焼装置の設計分野では、ノッキングによる装置の破損を回避するために熱エネルギーの回収率をある程度犠牲にせざるを得ない。また、動力発生装置（内燃機関）及び車両の設計分野では、ノッキングによる装置の破損を回避するために出力をある程度犠牲にせざるを得ない。
　このような燃焼装置、動力発生装置（内燃機関）及び車両等の設計分野では、ノッキングの予測精度が向上することによって熱エネルギーの回収率及び出力の犠牲分を縮小することが可能となる。すなわち、本実施形態に係るノッキング予測方法及びノッキング予測装置によれば、燃焼装置、動力発生装置及び車両等の設計で考慮されるべきノッキングに対する余裕度をこれまでよりも狭く見積もることが可能となる。
　したがって、本実施形態に係るノッキング予測方法及びノッキング予測装置で推定されたノッキングの発生条件に基づいて燃焼装置、動力発生装置（内燃機関）及び車両等を設計することにより、燃焼装置、動力発生装置（内燃機関）及び車両等の性能を向上させる設計を実現することが可能である。
　このような本実施形態に係るノッキング予測方法及びノッキング予測装置によれば、燃焼装置、動力発生装置及び車両等の性能を高めることが可能な燃焼装置の設計方法、動力発生装置の設計方法及び車両の設計方法を提供することが可能である。
　なお、このような燃焼装置の設計方法、動力発生装置の設計方法及び車両の設計方法では、ルイス数が１より大きい燃料を予測対象（計算対象）とすることが考えられる。すなわち、図４の特性図に示すように、ルイス数が１より小さい燃料はノッキングが発生する虞がないので、予測対象（計算対象）から除外してもよい。
　さらには、本実施形態に係るノッキング予測方法及びノッキング予測装置によれば、燃料の種別毎にノッキングを精度高く予測することが可能なので、ノッキングを回避するような燃料の設計に寄与することができる。したがって、本実施形態に係るノッキング予測方法及びノッキング予測装置によれば、燃料の性能を高めることが可能な燃料の設定方法を提供することが可能である。
　ここで、本実施形態では、ノッキング予測の評価対象燃料の一例として、「Hydrogen（水素）」や「ｎ-Heptane（へプタン）」等、単成分の燃料（純物質燃料）について本開示の有効性を説明した。しかしながら、本開示は、メタンやプロパン等、他の純物質燃料についても適用可能である。
　また、本実施形態は、評価対象燃料に関する１次元定常解析（シミュレーション）に基づいて評価対象燃料のノッキングを予測するものである。したがって、本開示は、１次元定常解析（シミュレーション）が可能なものであれば、純物質燃料に限定されず、複数の成分が混合する多物質燃料に適用することが可能である。
　周知のように、燃料の１次元定常解析（シミュレーション）を実行するためには、燃料の燃焼現象に関与する化学種と当該化学種に基づく素反応とを特定する必要がある。現時点で化学種と素反応とが特定されている燃料には、例えばSIP（革新的燃焼技術）サロゲート化学反応機構において定義された５種までの成分からなるサロゲート燃料がある。SIPサロゲート化学反応機構では、６種以上の成分からなる多物質燃料については化学種と素反応とが特定されていない状況である。
　しかしながら、６種以上の成分からなる多物質燃料についても、今後、化学種及び素反応の特定が十分に期待される。すなわち、６種以上の成分からなる多物質燃料についても、化学種と素反応とが特定されて１次元定常解析（シミュレーション）が可能となることが十分に想定され、本開示に基くノッキング予測が十分に可能になると推測される。
　続いて、本実施形態に係るエンジン制御装置（ＥＣＵ：Electronic Control Unit）及び車両（自動車）について、図８及び図９を参照して説明する。
　本実施形態に係るエンジン制御装置は、図８に示すＥＣＵ（Electronic Control Unit）１１である。また、本実施形態に係る車両は、同じく図８に示すように、車両操作装置１０とＥＣＵ１１とエンジン１２と通信装置１３とを備える自動車Ｍである。
　車両操作装置１０は、自動車Ｍの操作パネルである、運転者の入力操作に基づいて操作指示信号を生成してＥＣＵ１１に出力する。すなわち、操作指示信号は、自動車Ｍに対する運転者の操作指示を示す信号である。
　ＥＣＵ１１は、所謂組込型コンピュータであり、演算回路、記憶回路及び入出力回路を備える。このＥＣＵ１１は、記憶回路に予め記憶されたエンジン制御プログラムを実行することによってエンジン１２を制御する。
　すなわち、このＥＣＵ１１は、エンジン制御プログラムを演算回路で実行することによって各種のエンジン制御信号を生成する。このエンジン制御信号は、例えば給気バルブの開度及び燃料噴射バルブの開度やタイミング等、エンジン１２に対する各種の操作量（エンジン操作量）を示す信号である。ＥＣＵ１１は、このようなエンジン制御信号を入出力回路を用いてエンジン１２に出力することによってエンジン１２を制御する。
　エンジン１２は、燃料を燃焼室で燃焼させることによって動力を発生する内燃機関であり、複数のセンサからなるセンサ群１２ａを備えている。このセンサ群１２ａは、例えば給気バルブの開度、燃料噴射バルブの動作及び排ガスに含まれる二酸化炭素等、エンジン１２に関する各種の状態量をエンジン制御量として検出し、このエンジン制御量をセンサ信号としてＥＣＵ１１に出力する。
　ＥＣＵ１１は、上記エンジン制御信号の生成に際して、車両操作装置１０から入力される操作指示信号と上記センサ群１２ａから入力される各種のセンサ信号とを参照する。ＥＣＵ１１は、このような操作指示信号及びエンジン制御量をエンジン制御プログラムに従って情報処理することによって上記エンジン制御信号（エンジン操作量）を生成する。
　また、本実施形態におけるＥＣＵ１１は、図８に示すように、上述したノッキング予測装置Ａを備えている。ノッキング予測装置Ａは、上述したように記憶部１に予め記憶されたノッキング予測プログラムを演算部２で実行することによってノッキングの発生を予測するコンピュータである。
　ＥＣＵ１１におけるノッキング予測装置Ａは、自動車Ｍに燃料が給油されると、日時や位置情報から周囲の車輌から収集済みの情報に基づくビッグデータから燃料組成を導出し、その燃料組成についてノッキング予測プログラムを実行することによってノッキング発生予測結果を取得する。つまり、ノッキング発生予測結果の取得タイミングは、自動車Ｍに燃料が給油された直後であり、自動車Ｍが走行を開始する前のタイミングである。さらに、燃料組成の再検討は自動車Ｍが走行中にその運転条件から常に更新し、最初に予測した燃料組成との乖離が大きくなった際には、ノッキング発生予測結果も更新される。また、更新された情報は無線等のデータ通信により送信され、周囲の車輌や開発・改良に資するビッグデータとして蓄積される。
　本実施形態におけるＥＣＵ１１は、エンジン制御プログラムに加えてノッキング予測プログラムを予め記憶することによって、ノッキング予測装置Ａをエンジン制御機能に付帯する追加機能として内包する組込型コンピュータである。
　このようなＥＣＵ１１は、ノッキング予測プログラムの実行によって得られたノッキング予測結果をエンジン制御信号の生成において参照する。例えば、ＥＣＵ１１は、ノッキング予測プログラムの実行によってノッキングの発生が予測された場合、ノッキングを回避するようにエンジン制御信号を生成する。
　さらに、本実施形態におけるＥＣＵ１１は、ノッキング予測プログラムに基づくノッキング予測結果、エンジン制御信号が示すエンジン操作量及びセンサ群１２ａによって得られるエンジン制御量を付帯情報とともに外部のデータ収集装置Ｄに提供する。上記付帯情報は、エンジン制御量及びエンジン操作量以外のエンジン関連情報であり、例えば日付、時刻、外気温、季節、燃料組成及び自動車Ｍの走行位置である。
　すなわち、ＥＣＵ１１は、自動車Ｍの走行前に取得したノッキング予測結果と自動車Ｍの走行時において時系列的に順次得られるエンジン操作量、エンジン制御量及び付帯情報を学習データ（時系列データ）として記憶回路に順次記憶させる。ＥＣＵ１１は、このような学習データを所定のタイミングで記憶回路から読み出して通信装置１３に出力する。
　通信装置１３は、ＥＣＵ１１とデータ収集装置Ｄとの間に設けられ、ＥＣＵ１１とデータ収集装置Ｄとの通信を仲介する無線通信装置である。すなわち、通信装置１３は、ＥＣＵ１１から入力された学習データを無線信号に変換してデータ収集装置Ｄに送信する。
　このような学習データは、図９に示すように、ノッキング予測装置Ａを内包するＥＣＵ１１が備えられた複数の自動車Ｍ１～Ｍｎ（ｎ：２以上の自然数）からデータ収集装置Ｄに提供される。すなわち、データ収集装置Ｄには、複数の自動車Ｍ１～Ｍｎに関する学習データがビッグデータとして順次蓄積される。
　このようなビッグデータは、個々の学習データにエンジン操作量、エンジン制御量及び付帯情報に対応するノッキング予測結果が含まれているので、様々な活用方法が考えられる。
　例えば、ハイブリッドエンジンではノッキングが起こる条件でモータに切り替えたり、モーターアシストを強化してノッキングを回避するが、本実施形態におけるビッグデータは、このようなノッキング回避の設計検討に有益なデータとなり得る。
　また、エンジン性能のビッグデータ解析と本実施形態におけるビッグデータとを比較することによって、予測よりも性能が出ない場合に圧力センサの不具合を含む故障を検知し、エンジン１２の経時変化や故障・劣化の情報を提供することが可能となり、以ってメンテナンス性を向上させることが可能である。
　また、本実施形態におけるビッグデータを共有することによって、クラウドベースの車両診断や最適なエンジン運転パラメータに活かすことが可能である。また、ハイブリッド燃料システムとしてガソリンに、軽油等の炭化水素や異種燃料を混ぜることで着火性を上げ、以って燃焼速度を上げることで安全性の向上に繋げることが可能である。
　また、燃料が販売されている国や地域毎に燃料の種類が異なるため最適な燃焼条件を記述したデータをそろえることでノッキングを抑制しているが、本発明では自動車製造時にこうしたデータ採取や入力プロセスを省くことも可能になる。また、本発明に係る機能を搭載した自動車であれば、販売後に登場した今後e-fuelなどの新組成燃料でも、予めデータ入力するプロセスを経ることなく、ノッキング抑制の効果が得られる。さらには、エンジン１２の燃焼室に副室を設けた場合において着火性のコントロールを柔軟にし、以って熱効率を向上させることが可能である。
　本開示は、ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法に利用することができる。
　Ａ　ノッキング予測装置
　Ｄ　データ収集装置
　Ｍ　自動車（車両）
　１　記憶部
　２　演算部
　３　操作部
　４　出力部
　１０　車両操作装置
　１１　ＥＣＵ（エンジン制御装置）
　１２　エンジン
　１２ａ　センサ群
　１３　通信装置
　

Claims (16)

  1. 　予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、
    　前記予混合燃焼場において、ゼロ次元均質点火における正規化燃料質量分率及び正規化温度の時間的変化と１次元層流予混合火炎における空間的変化とが等価であるという知見に基づいてノッキングを予測するノッキング予測方法。
  2. 　予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、
    　正規化温度と燃料の正規化燃料質量分率との関係を示す特性図において、ルイス数が１より大きい燃料を前記ノッキングの可能性があるとし、前記ルイス数が１より小さい燃料を前記ノッキングの可能性がないとするノッキング予測方法。
  3. 　予混合燃焼場で発生し得るノッキングを予測するノッキング予測方法であって、
    　前記予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算する計算工程と、
    　前記予混合燃焼場に火炎が存在できなくなるために前記燃焼速度が取得不能になったか否かを判定する判定工程と、
    　最後に取得された前記燃焼速度に対応する前記入口温度を前記ノッキングの発生条件として取得する発生条件取得工程と
    　を有するノッキング予測方法。
  4. 　前記計算工程では、前記予混合燃焼場における未燃混合気、火炎の反応帯及び既燃焼気体を含む領域の大きさを、火炎が包含できる限界まで小さくして計算する請求項３に記載のノッキング予測方法。
  5. 　ルイス数が１より大きい燃料を予測対象とする請求項３又は４に記載のノッキング予測方法。
  6. 　予混合燃焼場で発生し得るノッキングを予測するノッキング予測装置であって、
    　前記予混合燃焼場の１次元定常解析によって入口温度を徐々に増加させた場合の燃焼速度を順次計算する計算手段と、
    　前記予混合燃焼場に火炎が存在できなくなるために前記燃焼速度が取得不能になったか否かを判定する判定手段と、
    　最後に取得された前記燃焼速度に対応する前記入口温度を前記ノッキングの発生条件として取得する発生条件取得手段と
    　を備えるノッキング予測装置。
  7. 　前記計算手段は、前記予混合燃焼場における未燃混合気、火炎の反応帯及び既燃焼気体を含む領域の大きさを、火炎が包含できる限界まで小さくして計算する請求項６に記載のノッキング予測装置。
  8. 　前記計算手段は、ルイス数が１より大きい燃料を計算対象とする請求項６又は７に記載のノッキング予測装置。
  9. 　請求項３又は４に記載のノッキング予測方法に基づいて推定された前記ノッキングの発生条件に基づいて燃焼装置を設計する燃焼装置の設計方法。
  10. 　請求項３又は４に記載のノッキング予測方法に基づいて推定された前記ノッキングの発生条件に基づいて動力発生装置を設計する動力発生装置の設計方法。
  11. 　請求項３又は４に記載のノッキング予測方法に基づいて推定された前記ノッキングの発生条件に基づいて車両を設計する車両の設計方法。
  12. 　請求項３又は４に記載のノッキング予測方法に基づいて推定された前記ノッキングの発生条件に基づいて燃料を設計する燃料の設計方法。
  13. 　請求項６又は７に記載のノッキング予測装置を備え、
    　前記ノッキング予測装置のノッキング予測結果に基づいてエンジン制御信号を生成するエンジン制御装置。
  14. 　外部との通信が自在な通信装置をさらに備え、
    　少なくとも前記ノッキング予測結果及び前記エンジン制御信号が示すエンジン操作量を前記通信装置を用いて外部に送信する請求項１３に記載のエンジン制御装置。
  15. 　エンジン制御量に加え、該エンジン制御量及び前記エンジン操作量以外のエンジン関連情報を前記外部に送信する請求項１４に記載のエンジン制御装置。
  16. 　請求項１３又は１４に記載のエンジン制御装置と、
    　該エンジン制御装置によって制御されるエンジンと、
    　を備える車両。
PCT/JP2024/020120 2023-06-06 2024-05-31 ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法 Ceased WO2024253047A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202480037229.4A CN121241194A (zh) 2023-06-06 2024-05-31 爆震预测方法、爆震预测装置、发动机控制装置、车辆、燃烧装置的设计方法、动力产生装置的设计方法、车辆的设计方法以及燃料的设计方法
JP2025526095A JPWO2024253047A1 (ja) 2023-06-06 2024-05-31

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
PCT/JP2023/020935 WO2024252510A1 (ja) 2023-06-06 2023-06-06 ノッキング予測方法、ノッキング予測装置、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法
JPPCT/JP2023/020935 2023-06-06

Publications (1)

Publication Number Publication Date
WO2024253047A1 true WO2024253047A1 (ja) 2024-12-12

Family

ID=93795329

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/JP2023/020935 Ceased WO2024252510A1 (ja) 2023-06-06 2023-06-06 ノッキング予測方法、ノッキング予測装置、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法
PCT/JP2024/020120 Ceased WO2024253047A1 (ja) 2023-06-06 2024-05-31 ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
PCT/JP2023/020935 Ceased WO2024252510A1 (ja) 2023-06-06 2023-06-06 ノッキング予測方法、ノッキング予測装置、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法

Country Status (3)

Country Link
JP (1) JPWO2024253047A1 (ja)
CN (1) CN121241194A (ja)
WO (2) WO2024252510A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119647344B (zh) * 2024-12-17 2025-09-02 中国科学技术大学 一种基于进化算法的简化化学反应动力学模型构建方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006200396A (ja) * 2005-01-19 2006-08-03 Hitachi Ltd 内燃機関の点火時期制御装置および点火時期制御ロジック作成装置
JP2007332934A (ja) * 2006-06-19 2007-12-27 Toyota Central Res & Dev Lab Inc 火花点火式内燃機関の最適点火時期設定方法及び火花点火式内燃機関の最適点火時期設定装置
JP2008031916A (ja) * 2006-07-28 2008-02-14 Toyota Motor Corp 内燃機関の自着火予測装置
CN112345259A (zh) * 2020-10-20 2021-02-09 同济大学 一种基于爆震自识别的汽油机虚拟标定方法
CN113673089A (zh) * 2021-07-23 2021-11-19 东风汽车集团股份有限公司 一种发动机性能确定方法、装置以及电子设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006200396A (ja) * 2005-01-19 2006-08-03 Hitachi Ltd 内燃機関の点火時期制御装置および点火時期制御ロジック作成装置
JP2007332934A (ja) * 2006-06-19 2007-12-27 Toyota Central Res & Dev Lab Inc 火花点火式内燃機関の最適点火時期設定方法及び火花点火式内燃機関の最適点火時期設定装置
JP2008031916A (ja) * 2006-07-28 2008-02-14 Toyota Motor Corp 内燃機関の自着火予測装置
CN112345259A (zh) * 2020-10-20 2021-02-09 同济大学 一种基于爆震自识别的汽油机虚拟标定方法
CN113673089A (zh) * 2021-07-23 2021-11-19 东风汽车集团股份有限公司 一种发动机性能确定方法、装置以及电子设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LAPALME DENIS; LEMAIRE ROMAIN; SEERS PATRICE: "Assessment of the method for calculating the Lewis number of H2/CO/CH4mixtures and comparison with experimental results", INTERNATIONAL JOURNAL OF HYDROGEN ENERGY, ELSEVIER, AMSTERDAM, NL, vol. 42, no. 12, 11 March 2017 (2017-03-11), AMSTERDAM, NL, pages 8314 - 8328, XP085008031, ISSN: 0360-3199, DOI: 10.1016/j.ijhydene.2017.01.099 *
YOUHI MORII; KAORU MARUTA: "What connects ignition and deflagration? -- On explosive transition of deflagration", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 18 January 2023 (2023-01-18), 201 Olin Library Cornell University Ithaca, NY 14853, XP091415012 *

Also Published As

Publication number Publication date
JPWO2024253047A1 (ja) 2024-12-12
CN121241194A (zh) 2025-12-30
WO2024252510A1 (ja) 2024-12-12

Similar Documents

Publication Publication Date Title
Lapuerta et al. Autoignition of blends of n-butanol and ethanol with diesel or biodiesel fuels in a constant-volume combustion chamber
Yıldız et al. Zero-dimensional single zone engine modeling of an SI engine fuelled with methane and methane-hydrogen blend using single and double Wiebe Function: A comparative study
Pourkhesalian et al. Alternative fuel and gasoline in an SI engine: A comparative study of performance and emissions characteristics
Irdmousa et al. Data-driven modeling and predictive control of combustion phasing for RCCI engines
Gravalos et al. Detection of fuel type on a spark ignition engine from engine vibration behaviour
WO2024253047A1 (ja) ノッキング予測方法、ノッキング予測装置、エンジン制御装置、車両、燃焼装置の設計方法、動力発生装置の設計方法、車両の設計方法及び燃料の設計方法
Hu et al. A study of liquid fuel injection and combustion in a constant volume vessel at diesel engine conditions
Ortiz-Soto et al. A comprehensive engine to drive-cycle modelling framework for the fuel economy assessment of advanced engine and combustion technologies
Desantes et al. Validity of the Livengood & Wu correlation and theoretical development of an alternative procedure to predict ignition delays under variable thermodynamic conditions
Hann et al. A Quasi-Dimensional SI Burn Rate Model for Predicting the Effects of Changing Fuel, Air-Fuel-Ratio, EGR and Water Injection
Vidal et al. Performance map of a LPG-diesel dual-fuel engine based on experimental and non-linear least squares determined wiebe function
JP2019157652A (ja) 内燃機関の制御装置
Chin et al. Development and validation of a reduced DME mechanism applicable to various combustion modes in internal combustion engines
Tougri et al. Knocking prediction in internal combustion engines via thermodynamic modeling: preliminary results and comparison with experimental data
Zhang et al. Reduced mechanism generation for methanol-based toluene reference fuel with combined reduction methods
Welscher et al. A Detailed Reaction Kinetics-Based Calculation Tool for Internal Combustion Engine-Related Ignition Processes
Rublewski et al. Modeling NO formation in spark ignition engines with a layered adiabatic core and combustion inefficiency routine
Nishii et al. Study on automatic adaptation for Control-oriented model of advanced diesel engine
Sadiq Al-Baghdadi Development of a pre-ignition submodel for hydrogen engines
Karim The onset of knock in gas fuelled spark ignition engines prediction and experiment
Fonseca et al. Tuning the parameters of ECFM-3Z combustion model for CFD 3D simulation of a two valve engine fueled with ethanol
Prince Development of an N-Propanol Mechanism and Its Application on HCCI Engine Simulations
Voicu et al. A numerical simulation of the influence of injection characteristics on performance and emissions of a tractor diesel engine
Floweday Two contrasting approaches to auto-ignition modelling for HCCI engines
Zarghami et al. Performance enhancement of spark ignition engines by using fractional order controller

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 24819277

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2025526095

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2025526095

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE