CN104679962B - Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines - Google Patents
Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines Download PDFInfo
- Publication number
- CN104679962B CN104679962B CN201510117124.2A CN201510117124A CN104679962B CN 104679962 B CN104679962 B CN 104679962B CN 201510117124 A CN201510117124 A CN 201510117124A CN 104679962 B CN104679962 B CN 104679962B
- Authority
- CN
- China
- Prior art keywords
- pressure
- time step
- sound wave
- sound
- engine cylinder
- 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.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Engines (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
The embodiment of the present invention provides a kind of combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines.Sound wave propagates the first pressure renewal amount to grid cell in IC engine cylinder again in the inventive method, including use KIVA calculating N time step current combustion point sources and N time steps to initial time step-length；IC engine cylinder internal face reflection described in N time steps is calculated to the second pressure renewal amount of grid cell in IC engine cylinder using SYSNOISE；First pressure renewal amount described in superposition obtains the total rise of pressure of grid cell in N time step IC engine cylinders with the second pressure renewal amount；The pressure curve in the IC engine cylinder is determined according to the total rise of the pressure of each time step.The embodiment of the present invention more accurately calculates the combustion noise in IC engine cylinder.
Description
Technical field
The present embodiments relate to data simulation analysis field, more particularly to one kind is based in KIVA and SYSNOISE couplings
The combustion noise prediction analysis method of combustion engine.
Background technology
In engine combustion process, when piston is run near top dead center, in-cylinder pressure can be raised drastically, be formed in cylinder
The higher-order of oscillation of pressure.The engine in-cylinder pressure higher-order of oscillation is an important factor for producing engine combustion noise.Combustion noise
Power can be used as assessment method by the use of the spectrum curve of in-cylinder pressure.So the pressure curve of in-cylinder combustion is obtained to analysis
Combustion noise has great significance.At present, it is that the engine specimen manufactured after processing is installed on engine pedestal to carry out
The combustion noise of engine is researched and analysed in noise testing.If simulation meanses can be passed through in the design time of engine specimen
Carry out the authentic and valid pressure curve obtained in engine cylinder, there will be very important meaning to estimating combustion noise, and
The cost and manpower of manufacture processing can be saved.
Apply to the software of engine in-cylinder combustion numerical simulation now, such as KIVA, AVL_FIRE and FLUENT, be all
Establish on the basis of flow field model, its pressure curve being calculated all is a smooth curve.Because do not account for
Sound field in the coupled combustion chamber of burning influences.The curve of cyclical fluctuations phase that the smooth curve that simulation obtains obtains with test measurement
Than there is certain distortion.Combustion noise can not be accurately estimated if smooth cylinder directly using these buckles line.
The content of the invention
The embodiment of the present invention provides a kind of combustion noise forecast analysis side based on KIVA and SYSNOISE coupling internal combustion engines
Method, to overcome the combustion noise predictive analysis results distortion for using simulation softward to calculate internal combustion engine in the prior art is too high to ask
Topic.By the research to thermoacoustic coupling mechanism, the rule of the exploration in-cylinder pressure higher-order of oscillation can effectively solve " distortion " and ask
Topic.In KIVA softwares, using chemical reaction formed Combustion Source as sound source, its caused sound wave in cylinder propagation distance with
The increase of time and increase, sound wave influence in the range of in-cylinder pressure change therewith.When sound wave reaches combustion chamber wall surface, and from
Wall produces reflection, and back wave influences the distribution of in-cylinder pressure field again.For the communication process of back wave, SYSNOISE is used
Software is simulated.Interface routine being write using fortran language in KIVA and calling SYSNOSIE, simulation calculates engine cylinder
Interior pressure field.The reason for by disclosing the in-cylinder pressure higher-order of oscillation, explores combustion in IC engine mechanism of noise generation, to estimate combustion
Burn noise and theories integration is provided.
The embodiments of the invention provide a kind of combustion noise forecast analysis based on KIVA and SYSNOISE coupling internal combustion engines
Method, including：
Sound in N time step current combustion point sources and N time steps to initial time step-length is calculated using KIVA
Ripple propagates the first pressure renewal amount to grid cell in IC engine cylinder again；
IC engine cylinder internal face reflection described in N time steps is calculated to grid cell in IC engine cylinder using SYSNOISE
Second pressure renewal amount；
First pressure renewal amount described in superposition obtains N time step IC engine cylinder Intranets with the second pressure renewal amount
The total rise of pressure of lattice unit；
The pressure curve in the IC engine cylinder is determined according to the total rise of the pressure of each time step.
Further, it is described that IC engine cylinder internal face reflection described in N time steps is calculated to internal combustion using SYSNOISE
The second pressure renewal amount of grid cell in machine cylinder, including：
Flow field grid overall in IC engine cylinder is made into the solid solid of local flow with the wall residing for exciting force to couple, flowed
Admittedly the FEM model coupled；
Each current excitations power of each time step is calculated respectively using SYSNOISE to all grids in IC engine cylinder
Pressure change amount caused by unit；
The pressure change amount of the current excitations power of each time step is stored to local；
According to the Acoustic Wave Propagation distance of N time step current excitations power grid cell in the IC engine cylinder described
First pressure knots modification is extracted in pressure change amount；
Continue propagation distance according to sound wave corresponding to each time step in the N time steps to initial time step-length to exist
The second pressure knots modification that sound wave continues to propagate is extracted in the pressure change amount；
First pressure knots modification described in superposition obtains N time step IC engine cylinder inwalls with the second pressure rise
The second pressure renewal amount to grid cell in IC engine cylinder is reflected in face.
Further, it is described that N time step current combustion point sources and N time steps are calculated to initial using KIVA
Sound wave propagates the first pressure renewal amount to grid cell in IC engine cylinder again in time step, including：
Processor calculates each current combustion point source of each time step to owning in IC engine cylinder respectively using KIVA
Pressure change amount caused by grid cell；
The processor stores the pressure rise amount of the current combustion point source of each time step to local；
The processor passes first according to the sound wave of N time step current combustion point sources grid cell in IC engine cylinder
The distance broadcast extracts the 3rd pressure change amount that sound wave is propagated first in the pressure change amount；
Processor sound wave according to corresponding to each time step in N time steps to initial time step-length continues to propagate
Distance extracted in the pressure change amount sound wave continue propagate the 4th pressure change amount；
The 4th pressure change that 3rd pressure change amount described in the processor superposition continues to propagate with the sound wave measures
To the first pressure renewal amount in N time step IC engine cylinders；
The processor determines the pressure curve in IC engine cylinder according to the pressure renewal amount of each time step.
Further, the processor is according to N time step current combustion point sources grid cell in IC engine cylinder
The distance that sound wave is propagated first extracts the 3rd pressure change amount that sound wave is propagated first in the pressure change amount, including：
Processor determines the total number of current combustion point source in N time steps according to KIVA chemical reaction subprogram；
The spatial coordinates calculation of each node in processor grid cell according to corresponding to the current combustion point source
Locus coordinate in the N time steps of the current combustion point source；
The processor calculates the velocity of sound of the N burnings point source N time steps according to small amplitude three-dimensional acoustic wave equation,
And the distance propagated according to the velocity of sound calculating N burnings point source in N time steps internal space；
The processor is true according to the locus coordinate of the current combustion point source and the distance of the Acoustic Wave Propagation
Unit in the sound that the current combustion point source influences in the fixed N time steps；
The processor extracts in the pressure change amount stored that sound wave is propagated first according to unit in the sound
Three pressure change amounts.
Further, processor sound according to corresponding to each time step in N time steps to initial time step-length
The distance that ripple continues to propagate extracts the 4th pressure change amount that sound wave continues to propagate in the pressure change amount, including：
The processor according to the velocity of sound determines N time steps, and into initial time step-length, each time step sound wave is again
The distance of propagation；
The processor is according to each time step current combustion point source locus coordinate, the current combustion point source
Propagation distance determines unit in the sound that sound wave corresponding to each time step is propagated again again for the distance of propagation and the sound wave；
Unit extracts sound wave in the pressure change amount stored in the sound that the processor is propagated again according to the sound wave
Continue the 4th pressure change amount propagated.
Further, the 3rd pressure change amount described in the processor superposition continues the 4th pressure of propagation with the sound wave
After knots modification obtains the first pressure renewal amount in N time step IC engine cylinders, in addition to：
The processor stores file corresponding to the pressure renewal amount of N time steps to local, the file
Including：Unit in the acoustic propagation distance of the current combustion point source corresponding to the N time steps, the sound of the Acoustic Wave Propagation
Total number, the burning locus coordinate of point source, grid cell pressure rise, the internal combustion engine in the IC engine cylinder
The initial number of grid cell in cylinder, in the IC engine cylinder grid cell termination number.
Further, the processor stores file corresponding to the pressure renewal amount of N time steps to local
Afterwards, in addition to：
The processor judges any burning point source in file corresponding to pressure renewal amount described in the N time steps
And all sound waves corresponding to any exciting force propagate whether unit in sound is more than grid cell whole in IC engine cylinder again, if
It is then to delete file corresponding to the burning point source or exciting force.
Further, the processor according to each time step current combustion point source locus coordinate, it is described work as
Propagation distance determines list in the sound that the sound wave of each time step is propagated again again for the distance of preceding burning point source propagation and the sound wave
Member, including：
The processor judges that the space length of any grid cell and current combustion point source in the IC engine cylinder is
The no distance influenceed more than current combustion point source is simultaneously less than the distance propagated again of sound wave, if, it is determined that the grid cell is
Unit in the sound that sound wave is propagated again, if not, it is determined that sound wave described in the grid cell N time steps is propagated again
The outer unit of sound.
The present embodiment by N time steps burning point source, N time steps to each time step of initial time step-length
Sound wave is being propagated to first pressure renewal amount caused by grid cell in IC engine cylinder and IC engine cylinder internal face exciting force, corresponding
Propagated again to second pressure renewal amount caused by grid cell in IC engine cylinder, the superposition first pressure in the sound wave of the exciting force
Renewal amount and second pressure renewal amount, obtain the total rise of pressure of grid cell in N time step IC engine cylinders, and according to
The total rise of pressure corresponding to each time step determines the pressure curve in IC engine cylinder, realizes and internal combustion engine in-cylinder combustion is made an uproar
The prediction of sound, improve the accuracy for calculating internal combustion engine in-cylinder pressure.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are this hairs
Some bright embodiments, for those of ordinary skill in the art, without having to pay creative labor, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the combustion noise prediction analysis method flow chart that the present invention couples internal combustion engine based on KIVA and SYSNOISE；
Fig. 2 is inventive combustion engine inside wall of cylinder face back wave schematic diagram；
Fig. 3 is that KIVA of the present invention runs prompt statement schematic diagram；
Fig. 4 is sound wave of the present invention to interior cell schematics at the sound；
Fig. 5 is processor storage file schematic diagram of the present invention；
Fig. 6 is that sound wave of the present invention propagates schematic diagram again；
Fig. 7 is the schematic diagram that present invention burning point source forms Acoustic Wave Propagation；
Fig. 8 sends the change schematic diagram that sound wave is propagated again for present invention burning point source；
Fig. 9 is that back wave hesitation propagates schematic diagram at inventive combustion engine wall；
Figure 10 is the in-cylinder pressure field distribution schematic diagram that prior art uses KIVA to obtain；
Figure 11 is that combustion noise prediction analysis method of the present invention based on KIVA and SYSNOISE coupling internal combustion engines obtains
In-cylinder pressure field distribution schematic diagram；
Figure 12 is that acoustic simulation of the present invention is calculated with testing the in-cylinder pressure curve comparison schematic diagram measured；
Figure 13 is in-cylinder pressure curve Spectrum Analysis Comparison schematic diagram of the present invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
Part of the embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Fig. 1 is the combustion noise prediction analysis method flow chart that the present invention couples internal combustion engine based on KIVA and SYSNOISE,
As shown in figure 1, this method embodiment, including：
Step 101, N time step current combustion point sources and N time steps calculated to initial time using KIVA
Sound wave propagates the first pressure renewal amount to grid cell in IC engine cylinder again in step-length；
Step 102, using SYSNOISE calculate N time steps described in IC engine cylinder internal face reflection in IC engine cylinder
The second pressure renewal amount of grid cell；
First pressure renewal amount described in step 103, superposition obtains N time step internal combustions with the second pressure renewal amount
The total rise of the pressure of grid cell in machine cylinder；
Step 104, determine according to the total rise of the pressure of each time step pressure in the IC engine cylinder
Curve.
Specifically, Theory of Combustion is combined by the present invention with acoustic theory, on the basis of KIVA, by combustion in IC engine
Sound source of the grid cell as generation combustion noise corresponding to the burning point source of the chemical reaction occurred in the space lattice of room, those
Sound source produces sound wave in combustion.According to the coverage of Acoustic Wave Propagation by sound source peripheral region be divided into sound unit and
The outer unit of sound, the pressure of grid cell in IC engine cylinder is calculated according to the different zones residing for it.And the Acoustic Wave Propagation
Reflection sound wave can be produced during to IC engine cylinder internal face, so as to which the second pressure produced to grid cell in IC engine cylinder updates
Amount, the second pressure renewal amount, last superposition first pressure renewal amount and the second pressure renewal amount are calculated using SYSNOISE,
The total rise of pressure of grid cell in N time step IC engine cylinders is obtained, so as to always be risen according to the pressure of each time step
A large amount determines the pressure curve in IC engine cylinder, and carries out combustion in IC engine noise by the spectrum analysis to in-cylinder pressure curve
Prediction.
The sound wave that in-cylinder combustion noise source is sent, increased over time using Combustion Source as starting point distally empty step by step
Between propagate.With the change of engine crankshaft corner, i.e., as the accumulation of time step dt in KIVA, sound wave are distally propagated
.KIVA calling SYSNOISE service charts, the sound wave sent by Combustion Source can be with slow after the fuel-air mixture in cylinder starts burning
Residual effect should be gradually transferred to the grid cell at combustion chamber wall surface, can be sent out close to the pressure of grid cell at combustion chamber wall surface this moment
Changing.The change of pressure, that is, the exciting force to combustion chamber wall surface is formd, the power taps wall, and sound wave starts to reflect afterwards.
As shown in Figure 2.After sound wave reaches combustion chamber wall surface, reflected, back wave continues progressively to propagate in cylinder, when back wave reaches
To after combustion chamber wall surface, then reflect again.So the communication process of sound wave is several combustion noise sources in engine cylinder
The sound wave sent the sound wave incident that space is carried out in cylinder and the process of reflection.By that analogy, with the progress of burning, burning
Before the pressure wave that source is sent reaches combustion chamber peripheral wall surface, the pressure wave of Combustion Source is all propagated in cylinder.Work as Combustion Source
After the pressure wave sent reaches wall, back wave is all propagated in cylinder.So the unit of some opening position is being worked as in cylinder
In preceding time step, both by multiple incident wave actions of Combustion Source, also influenceed by multiple back waves at combustion chamber wall surface.
As shown in figure 3, KIVA calls SYSNOISE service charts after the fuel-air mixture in cylinder starts burning, by Combustion Source
The sound wave sent can be gradually transferred to the grid cell at combustion chamber wall surface with hesitation, this moment close to net at combustion chamber wall surface
The pressure of lattice unit can change.The change of pressure, that is, the exciting force to combustion chamber wall surface being formd, the power taps wall,
Sound wave starts to reflect afterwards.
In order to realize when KIVA is run, the instant percussion wall of grid cell at the combustion chamber wall surface in each time step
The effect in face, therefore SYSNOISE is vibrated-noise software is combined with KIVA programs.This need to write corresponding interface routine with
Realize the communication between them.Because KIVA programs are programmed using formula translation, therefore writing interface routine
When be necessary also to use formula translation.Windows API are operating system application programming interfaces, as progress application program
The System Programming interface write.The major function of interface routine is to judge that SYSNOISE is according to the window ejected on computer desktop
No normal startup.When starting successfully, prompt statement can be embodied in the dos pages of KIVA operations.
The present embodiment by N time steps burning point source, N time steps to each time step of initial time step-length
Sound wave is being propagated to first pressure renewal amount caused by grid cell in IC engine cylinder and IC engine cylinder internal face exciting force, corresponding
Propagated again to second pressure renewal amount caused by grid cell in IC engine cylinder, the superposition first pressure in the sound wave of the exciting force
Renewal amount and second pressure renewal amount, obtain the total rise of pressure of grid cell in N time step IC engine cylinders, and according to
The total rise of pressure corresponding to each time step determines the pressure curve in IC engine cylinder, realizes and internal combustion engine in-cylinder combustion is made an uproar
The prediction of sound, improve the accuracy for calculating internal combustion engine in-cylinder pressure.
Further, it is described that N time step current combustion point sources and N time steps are calculated to initial using KIVA
Sound wave propagates the first pressure renewal amount to grid cell in IC engine cylinder again in time step, including：
Processor calculates each current combustion point source of each time step to owning in IC engine cylinder respectively using KIVA
Pressure change amount caused by grid cell；
The processor stores the pressure rise amount of the current combustion point source of each time step to local；
The processor passes first according to the sound wave of N time step current combustion point sources grid cell in IC engine cylinder
The distance broadcast extracts the 3rd pressure change amount that sound wave is propagated first in the pressure change amount；
Processor sound wave according to corresponding to each time step in N time steps to initial time step-length continues to propagate
Distance extracted in the pressure change amount sound wave continue propagate the 4th pressure change amount；
The 4th pressure change that 3rd pressure change amount described in the processor superposition continues to propagate with the sound wave measures
To the first pressure renewal amount in N time step IC engine cylinders；
The processor determines the pressure curve in IC engine cylinder according to the pressure renewal amount of each time step.
Specifically, Acoustic Wave Propagation has hesitation, in the characteristic the present embodiment, processor uses KIVA first
Each current combustion point source corresponding to random time step-length N to each time step in initial time step-length is calculated respectively to internal combustion
Pressure change amount caused by all grid cells in machine cylinder, and the pressure change amount is stored to local.If the N time steps
It is not initial time step-length, the sound wave that corresponding current combustion point source is sent in N-1 time steps is inevitable in the N times
Had occurred and that during step-length and continue to propagate.Similarly, corresponding to N-2 time steps there is also sound wave propagation phenomenon again, by that analogy,
Processor calculate N-1 time steps to initial time step-length in corresponding to each time step in the IC engine cylinder by sound wave
The pressure rise amount of unit in the sound of formation is propagated again.Then processor is pressed according to corresponding to N time step current combustion point sources
When power rise and N time steps to initial time step-length sound wave propagate corresponding pressure rise amount Iterative calculation method N again
Between step-length to the pressure renewal amount of grid cell in the IC engine cylinder.Finally, processor is according to N time steps to initial time
The pressure renewal amount of each time step draws in-cylinder pressure curve in step-length, and FFT spectrum analysis is carried out to the curve so as to predict
Analyze combustion in IC engine noise.
Further, the processor is according to N time step current combustion point sources grid cell in IC engine cylinder
The distance that sound wave is propagated first extracts the 3rd pressure change amount that sound wave is propagated first in the pressure change amount, including：
Processor determines the total number of current combustion point source in N time steps according to KIVA chemical reaction subprogram；
Specifically, judged by the thermal discharge of previous cycle in IC engine cylinder with the difference of previous circulation thermal discharge
Whether chemical reaction occurs.If the difference of chemical heat release amount is more than 0, show that oil droplet starts to burn in cylinder, and determine to work as
The total number of preceding burning point source.
The spatial coordinates calculation of each node in processor grid cell according to corresponding to the current combustion point source
Locus coordinate in the N time steps of the current combustion point source；
The processor calculates the velocity of sound of the N burnings point source N time steps according to small amplitude three-dimensional acoustic wave equation,
And the distance propagated according to the velocity of sound calculating N burnings point source in N time steps internal space；
The processor is true according to the locus coordinate of the current combustion point source and the distance of the Acoustic Wave Propagation
Unit in the sound that the current combustion point source influences in the fixed N time steps；
The processor extracts in the pressure change amount stored that sound wave is propagated first according to unit in the sound
Three pressure change amounts.
Specifically, the heat release to be chemically reacted in IC engine cylinder, it is exactly the oil droplet combustion heat release of fuel-air mixture in fact,
The oil droplet number to chemically react is exactly the number of current combustion point source, so as to the oil droplet that chemically reacts in cylinder it is three-dimensional
Correspond to several units in space lattice, can be calculated often by the space coordinates of each node in those grid cells
The locus coordinate of one burning point source.According to formula
The velocity of sound of the burning point source is calculated, wherein, the γ is specific heat ratio, the p0For the pressure of gaseous mixture in cylinder, institute
State ρ0For the density of gaseous mixture, KIVA can be used to calculate the cumulative volume of Flow field, the gross mass and cylinder Intranet of Flow field
The summation of lattice unit acting.The average pressure in cylinder is obtained by total work divided by cumulative volume, is obtained using gross mass divided by cumulative volume
Averag density in cylinder.And the distance of Acoustic Wave Propagation is tried to achieve according to the velocity of sound.
As shown in figure 4, the center of circle of spheroid corresponding to unit is this in the sound that current combustion point source influences in N time steps
Burn point source locus coordinate, and radius is the distance of Acoustic Wave Propagation.
Further, processor sound according to corresponding to each time step in N time steps to initial time step-length
The distance that ripple continues to propagate extracts the 4th pressure change amount that sound wave continues to propagate in the pressure change amount, including：
The processor according to the velocity of sound determines N time steps, and into initial time step-length, each time step sound wave is again
The distance of propagation；
The processor is according to each time step current combustion point source locus coordinate, the current combustion point source
Propagation distance determines unit in the sound that sound wave corresponding to each time step is propagated again again for the distance of propagation and the sound wave；
Specifically, the processor judges the sky of any grid cell and current combustion point source in the IC engine cylinder
Between distance whether be more than the distance that current combustion point source influences and the distance propagated again less than sound wave, if, it is determined that the net
Unit in the sound that lattice unit is propagated again for sound wave, if not, it is determined that the grid cell is sound wave described in the N time steps
The outer unit of sound propagated again.
Unit extracts sound wave in the pressure change amount stored in the sound that the processor is propagated again according to the sound wave
Continue the 4th pressure change amount propagated.
Specifically, the present embodiment is that sound wave is being propagated to the net in IC engine cylinder in processor calculating N time steps
The pressure rise amount of lattice unit：As shown in fig. 6, because sound wave has hesitation, it is therefore, current in N-1 time steps
Burning point source sound wave in N time steps is propagated again equally produces pressure rise amount to grid cell in IC engine cylinder.Similarly,
Current combustion point source in the N-2 time steps sound wave in N time steps is propagated equally to internal combustion engine Intranet lattice unit again
Produce pressure rise amount.The pressure rise amount superposition that the pressure rise amount is propagated first with current combustion point source sound wave is to internal combustion engine
Grid cell produces pressure renewal amount in cylinder.As shown in fig. 7, embodying burning point source forms the influence that sound wave continues to propagate.
In current time step, the sound wave coverage that burning point source I is sent is a-quadrant, and the sound wave that burning point source II is sent influences model
Enclose for B area, in previous time step, the sound wave coverage that burning point source II is sent is C regions.Sheng Nei cells Ds area
Domain representation is influenceed and burnt corresponding to point source II by the first time acoustic propagation of burning point source I simultaneously in current time step
The influence propagated again of sound wave；In sound unit E region representations in current time step only by burning point source II sound wave again
The influence of propagation；Unit F region representation is also not affected by the influence of burning point source Acoustic Wave Propagation in current time step in sound.It is right
The pressure change amount step that the pressure rise amount method and the sound wave of current combustion point source propagated again in calculating sound wave are propagated first
Difference is that element method is in the sound that sound wave is propagated again：Processor judges any grid cell and sound wave in IC engine cylinder
Whether the space length propagated again is more than the distance that N-1 time step sound waves are propagated again and is less than N time steps sound wave again
The distance of propagation, if, it is determined that unit in the sound that the grid cell is propagated again for N time step sound waves, if not, it is determined that
The grid cell is the outer unit of sound that the N time step sound waves are propagated again.So as to processor superposition current combustion point source
The pressure rise that first pressure rise and sound wave are propagated again measures corresponding current combustion point in the N time steps
Propagate the pressure renewal amount to grid cell in the transmitter cylinder again with sound wave in source.
Further, the 3rd pressure change amount described in the processor superposition continues the 4th pressure of propagation with the sound wave
After knots modification obtains the first pressure renewal amount in N time step IC engine cylinders, in addition to：
The processor stores file corresponding to the pressure renewal amount of N time steps to local, the file
Including：Unit in the acoustic propagation distance of the current combustion point source corresponding to the N time steps, the sound of the Acoustic Wave Propagation
Total number, the burning locus coordinate of point source, grid cell pressure rise, the internal combustion engine in the IC engine cylinder
The initial number of grid cell in cylinder, in the IC engine cylinder grid cell termination number.
Specifically, after the processor obtains internal combustion engine in-cylinder pressure renewal amount corresponding to N time steps, place
Reason device can not store the time to local into corresponding N files, as shown in figure 5, when N corresponds to different value, can exist
Corresponding to the file of different time steps.The storage content of this document is：N time steps correspond to the burning point source and propagated
Distance, Acoustic Wave Propagation sound in unit total number, the locus coordinate of point source of burning, grid cell pressure in IC engine cylinder
The initial number of grid cell in rise, the IC engine cylinder, in the IC engine cylinder grid cell termination number.This article
Part propagates the pressure renewal amount to grid cell in IC engine cylinder for future time step size computation sound wave again.
Further, it is described that IC engine cylinder internal face reflection described in N time steps is calculated to internal combustion using SYSNOISE
The second pressure renewal amount of grid cell in machine cylinder, including：
Flow field grid overall in IC engine cylinder is made into the solid solid of local flow with the wall residing for exciting force to couple, flowed
Admittedly the FEM model coupled；
Each current excitations power of each time step is calculated respectively using SYSNOISE to all grids in IC engine cylinder
Pressure change amount caused by unit；
The pressure change amount of the current excitations power of each time step is stored to local；
According to the Acoustic Wave Propagation distance of N time step current excitations power grid cell in the IC engine cylinder described
First pressure knots modification is extracted in pressure change amount；
Continue propagation distance according to sound wave corresponding to each time step in the N time steps to initial time step-length to exist
The second pressure knots modification that sound wave continues to propagate is extracted in the pressure change amount；
First pressure knots modification described in superposition obtains N time step IC engine cylinder inwalls with the second pressure rise
The second pressure renewal amount to grid cell in IC engine cylinder is reflected in face.
Specifically, for the hesitation of real simulated reflections sound wave, the fluid structure interaction mode used strikes to be local
Striking.Calculated first by KIVA can search out it is a series of close to the unit for having obvious pressure change at combustion chamber wall surface.
The pressure change of these hexahedral elements establishes a solid wall on the surface that unit contacts with combustion chamber wall surface as exciting force
Face, Flow field grid couple with this solid wall surface, as shown in Figure 9.So, it can simulate and calculate one close to burning
To the influence of sound field in cylinder when back wave caused by unit at locular wall face is propagated with hesitation.Meanwhile set in cylinder and mix
The net boundary impedance of gas is the impedance value of gaseous mixture, and such back wave just can be projected away completely, without producing combustion chamber
The reverberation effect of wall reflection.
IC engine cylinder internal face exciting force specifically calculated for the second pressure renewal amount of grid cell in IC engine cylinder
Journey and current combustion point source are similar, and here is omitted.
Further, the processor stores file corresponding to the pressure renewal amount of N time steps to local
Afterwards, in addition to：
The processor judges any burning point source in file corresponding to pressure renewal amount described in the N time steps
And all sound waves corresponding to any exciting force propagate whether unit in sound is more than grid cell whole in IC engine cylinder again, if
It is then to delete file corresponding to the burning point source or exciting force.
Specifically, for sound wave is propagated again, when reaching sometime step-length, scope meeting that the sound wave is propagated again
The grid cell in IC engine cylinder, therefore, the sound wave is propagated again not to be changed to internal combustion engine in-cylinder pressure, such as Fig. 8 institutes
Show, respectively the scope of N, N+1, N+2 and N+3 time step burning point source Acoustic Wave Propagation.As can be seen here, pushing away over time
Move, it is increasing that sound wave propagates the coverage sent again.Now, if processor continues to propagate the sound wave again corresponding file
Carry out calculating the waste that can cause resource, therefore, processor judges that those sound waves are again in N time steps to initial time step-length
Whether unit includes grid cell whole in IC engine cylinder in whole sound corresponding to propagation, if so, then deleting the time step
The interior sound wave propagates corresponding file again, so as to improve operating efficiency.
As shown in Figure 10, obtained using former KIVA programs, color is than more uniform in cylinder, each position in this explanation cylinder
Relatively, in-cylinder pressure is distributed than more uniform the grid cell pressure value at place.Figure 11 is calculated using acoustic simulation program
The in-cylinder pressure field pattern arrived.As seen from the figure, in the region of oil droplet burning, it may appear that the face different from other positions in cylinder
Color, this shows partial high pressure occur in oil droplet burning near zone.Figure 10 and Figure 11 contrasts understand that Combustion Source energy is in sluggishness
Propagation and Combustion Source energy being evenly distributed in cylinder under effect, the influence to in-cylinder pressure field is clearly.
As shown in figure 12, what is provided is the song that the in-cylinder pressure of acoustic simulation calculating and the in-cylinder pressure of test measurement contrast
Line.The pressure that acoustic simulation calculates is the hesitation for considering sound wave incident ripple and back wave in communication process.This two
All there is obvious pressure oscillation in curve, and oscillation amplitude occurs mainly near top dead center.In combustion process of the internal-combustion engine,
Constant volume burning can be approximately considered as only near top dead center.With descent of piston, although burning is still being carried out, due to burning
Room volume increase, pressure oscillation are no longer fairly obvious.The result of acoustic simulation and the result of test measurement have preferable uniformity.
As shown in figure 13, the comparison diagram that the sound pressure level curve after spectrum analysis is carried out to in-cylinder pressure curve is given.This three curves
Respectively from KIVA simulations, acoustic simulation and test measurement.With the cylinder internal pressure of sample frequency 100KHz collection experiment diesel engines
Force signal, and FFT spectrum analysis is done to cylinder pressure signal, it is 20 to take reference pressure, and in-cylinder pressure level spectrum curve is as schemed
In 12 shown in blue curve.The in-cylinder pressure spectrum curve of test measurement is in 3.8KHz；4.4KHz；6.1KHz；7.0KHz；
8.1KHz；Occurs peak value at 9.3KHz.Red curve is the result of spectrum analysis of acoustic simulation in Figure 13, itself and test measurement
As a result it is coincide substantially in crest frequency point corresponding to spectrum curve.The spectrum analysis curve HFS that KIVA is simulated in Figure 13 does not have
There is obvious peak value, and HFS is the deciding factor of combustion noise.Therefore, predict that burning is made an uproar using KIVA simulations
Sound can bring larger error.In whole frequency range, contrast in-cylinder pressure level amplitude is visible, in low-frequency range, acoustic mode analog values
Slightly above test measurements；In Mid Frequency, acoustic mode analog values and test measurements mutual are intersected；In high band, acoustic mode analog values
Slightly above test measurements.And in full frequency band, the analogue value of the test measurements obviously higher than KIVA.The low-frequency range of spectrum curve
(frequency<In the range of 100Hz) amplitude size is depending on the maximum combustion pressure in cylinder, and, acoustic simulation obtains from Figure 12 also
The peak value of in-cylinder pressure curve be slightly above the result of test measurement, therefore be slightly above test measurement in low-frequency range, acoustic mode analog values
Value.Mid Frequency (the 100Hz of spectrum curve<Frequency<In the range of 2500Hz) amplitude size depends on the rate of pressure rise, 2500~
The amplitude of 4000Hz frequency bands is mainly together decided on by the rate of pressure rise and pressure high-frequency oscillation.And the frequency band higher than 4000Hz
It is due to caused by the higher-order of oscillation of in-cylinder pressure.It can be proved from the medium-high frequency section in Figure 13, the acoustic simulation in Figure 12
Frequency of oscillation is higher than the frequency of oscillation of test measurement, and test measurement is higher than what KIVA was simulated.Spectrum curve is in test measurement
It is because KIVA is not combined with CHEMKIN softwares, in chemical reaction part with the difference of acoustic simulation between the two
Computation model shows slightly coarse, and this can cause to produce error when calculating combustion noise source.In addition, SYSNOISE sound field nets can only be set
The integral damping of lattice model, and the fluid damping of each grid cell of setting that can not be in detail, this will necessarily be to back wave complete
Energy loss when being propagated in the grid cell of field produces certain influence.All in all, acoustic simulation result and test measurements
The degree of approach is preferable.
The in-cylinder pressure sound pressure level frequency spectrum obtained according to numerical simulation calculation can be estimated effectively to combustion noise.
Because the transmission letter of engine combustion noise can be passed through between combustion noise sound pressure level frequency spectrum and in-cylinder pressure sound pressure level frequency spectrum
Number opening relationships.Combustion noise sound pressure level frequency spectrum is equal to the product of in-cylinder pressure sound pressure level frequency spectrum and transmission function.And engine fires
The transmission function for burning noise does not change with factors such as rotating speed, load and igniting (oil spout) advance angles, only by internal combustion engine itself
Structures shape.So effectively burning can be estimated by accurately calculating in-cylinder pressure sound pressure level frequency spectrum and transmission function
The sound pressure level of noise.
Finally it should be noted that：Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations；To the greatest extent
The present invention is described in detail with reference to foregoing embodiments for pipe, it will be understood by those within the art that：Its according to
The technical scheme described in foregoing embodiments can so be modified, either which part or all technical characteristic are entered
Row equivalent substitution；And these modifications or replacement, the essence of appropriate technical solution is departed from various embodiments of the present invention technology
The scope of scheme.
Claims (7)
- A kind of 1. combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines, it is characterised in that including：Using sound wave in KIVA calculating N time step current combustion point sources and N time steps to initial time step-length again Propagate the first pressure renewal amount to grid cell in IC engine cylinder；IC engine cylinder internal face reflection described in N time steps is calculated to the of grid cell in IC engine cylinder using SYSNOISE Two pressure renewal amounts, including：Wall residing for flow field grid overall in IC engine cylinder and exciting force is made into the solid solid of local flow Coupling, obtains the FEM model of fluid structurecoupling；Calculate each current excitations of each time step respectively using SYSNOISE Power is to pressure change amount caused by all grid cells in IC engine cylinder；By the pressure of the current excitations power of each time step Knots modification is stored to local；According to the Acoustic Wave Propagation of N time step current excitations power grid cell in the IC engine cylinder Distance extracts first pressure knots modification in the pressure change amount；According in the N time steps to initial time step-length Sound wave corresponding to each time step continues propagation distance and the second pressure that sound wave continues to propagate is extracted in the pressure change amount Knots modification；First pressure knots modification described in superposition obtains N time step IC engine cylinder inwalls with the second pressure knots modification The second pressure renewal amount to grid cell in IC engine cylinder is reflected in face；First pressure renewal amount described in superposition obtains grid list in N time step IC engine cylinders with the second pressure renewal amount The total rise of pressure of member；The pressure curve in the IC engine cylinder is determined according to the total rise of the pressure of each time step.
- 2. according to the method for claim 1, it is characterised in that described that N time step current combustions are calculated using KIVA Sound wave is propagated to the first pressure of grid cell in IC engine cylinder more again in point source and N time steps to initial time step-length New amount, including：Processor calculates each current combustion point source of each time step to all grids in IC engine cylinder respectively using KIVA Pressure change amount caused by unit；The processor stores the pressure rise amount of the current combustion point source of each time step to local；The processor is propagated first according to N time step current combustion point sources sound wave of grid cell in IC engine cylinder Distance extracts the 3rd pressure change amount that sound wave is propagated first in the pressure change amount；Processor sound wave according to corresponding to each time step in N time steps to initial time step-length continue propagate away from From in the pressure change amount extract sound wave continue propagate the 4th pressure change amount；The 4th pressure change that 3rd pressure change amount described in the processor superposition continues to propagate with the sound wave is measured to N First pressure renewal amount in time step IC engine cylinder；The processor determines the pressure curve in IC engine cylinder according to the pressure renewal amount of each time step.
- 3. according to the method for claim 2, it is characterised in that the processor is according to N time steps current combustion point The source distance that the sound wave of grid cell is propagated first in IC engine cylinder is extracted sound wave in the pressure change amount and propagated first The 3rd pressure change amount, including：Processor determines the total number of current combustion point source in N time steps according to KIVA chemical reaction subprogram；In processor grid cell according to corresponding to the current combustion point source described in the spatial coordinates calculation of each node Locus coordinate in the N time steps of current combustion point source；The processor calculates the velocity of sound of the N burnings point source N time steps, and root according to small amplitude three-dimensional acoustic wave equation The distance propagated according to the velocity of sound calculating N burnings point source in N time steps internal space；The processor determines institute according to the locus coordinate of the current combustion point source and the distance of the Acoustic Wave Propagation State unit in the sound that the current combustion point source influences in N time steps；The processor extracts the 3rd pressure that sound wave is propagated first according to unit in the sound in the pressure change amount stored Power knots modification.
- 4. according to the method for claim 2, it is characterised in that the processor is according to N time steps to initial time The distance that sound wave corresponding to each time step continues to propagate in step-length extracts sound wave in the pressure change amount and continues what is propagated 4th pressure change amount, including：The processor according to the velocity of sound determine that N time steps each time step sound wave into initial time step-length propagates again away from From；The processor is propagated according to each time step current combustion point source locus coordinate, the current combustion point source Distance and the sound wave propagation distance determines unit in the sound that sound wave corresponding to each time step is propagated again again；Unit extracts sound wave continuation in the pressure change amount stored in the sound that the processor is propagated again according to the sound wave The 4th pressure change amount propagated.
- 5. according to the method for claim 2, it is characterised in that the 3rd pressure change amount and institute described in the processor superposition The 4th pressure change that sound wave continues to propagate is stated to measure to after the first pressure renewal amount in N time step IC engine cylinders, Also include：The processor stores file corresponding to the pressure renewal amount of N time steps to local, and the file includes： Unit is always individual in the acoustic propagation distance of the current combustion point source corresponding to the N time steps, the sound of the Acoustic Wave Propagation Grid cell pressure rise in several, the described burning locus coordinate of point source, the IC engine cylinder, in the IC engine cylinder The termination number of grid cell in the initial number of grid cell and the IC engine cylinder.
- 6. according to the method for claim 2, it is characterised in that the processor by the pressure of N time steps more File corresponding to new amount is stored to local, in addition to：The processor judge any burning point source in file corresponding to pressure renewal amount described in the N time steps and All sound waves corresponding to any exciting force propagate whether unit in sound is more than grid cell whole in IC engine cylinder again, if so, Then delete file corresponding to the burning point source or exciting force.
- 7. according to the method for claim 4, it is characterised in that the processor is according to each time step current combustion point source Propagation distance determines each time step again for locus coordinate, the distance of current combustion point source propagation and the sound wave Unit in the sound that sound wave is propagated again, including：The processor judges whether the space length of any grid cell and current combustion point source in the IC engine cylinder is big The distance influenceed in current combustion point source and the distance propagated again less than sound wave, if, it is determined that the grid cell is sound wave Unit in the sound propagated again, if not, it is determined that the grid cell is outside the sound that sound wave described in the N time steps is propagated again Unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510117124.2A CN104679962B (en) | 2015-03-17 | 2015-03-17 | Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510117124.2A CN104679962B (en) | 2015-03-17 | 2015-03-17 | Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104679962A CN104679962A (en) | 2015-06-03 |
| CN104679962B true CN104679962B (en) | 2018-01-02 |
Family
ID=53315001
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510117124.2A Expired - Fee Related CN104679962B (en) | 2015-03-17 | 2015-03-17 | Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104679962B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109443778B (en) * | 2018-09-29 | 2020-11-10 | 江苏大学 | Engine sound quality prediction method |
| CN113343542B (en) * | 2021-07-14 | 2023-01-31 | 义乌吉利动力总成有限公司 | Optimal design method of intake manifold and intake manifold |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103147868A (en) * | 2013-03-07 | 2013-06-12 | 天津大学 | Cavity resonance-based control method for reducing combustion noises of combustion engines |
| CN104063553A (en) * | 2014-07-04 | 2014-09-24 | 上汽通用五菱汽车股份有限公司 | Optimization design method of combustion system of engine |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2898411B1 (en) * | 2006-03-08 | 2008-05-16 | Inst Francais Du Petrole | REAL-TIME ESTIMATION METHOD OF ENGINE COMBUSTION PARAMETERS FROM VIBRATORY SIGNALS |
-
2015
- 2015-03-17 CN CN201510117124.2A patent/CN104679962B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103147868A (en) * | 2013-03-07 | 2013-06-12 | 天津大学 | Cavity resonance-based control method for reducing combustion noises of combustion engines |
| CN104063553A (en) * | 2014-07-04 | 2014-09-24 | 上汽通用五菱汽车股份有限公司 | Optimization design method of combustion system of engine |
Non-Patent Citations (4)
| Title |
|---|
| The acoustic hysteresis effect in KIVA during the combustion process of internal-combustion engines;Wei Li等;《Institution of Mechanical Engineers》;20140127;第228卷(第5期);第1-10页 * |
| 内燃机燃烧过程中高频压力振荡的热声耦合机理;李维等;《内燃机学报》;20140925;第32卷(第5期);正文第455页第2段，第456页第2.2节-2.3节，图1-图2 * |
| 基于KIVA和SYSNOISE的柴油机燃烧振荡数值模拟计算;薛冬新等;《内燃机》;20120911(第4期);第8-12页 * |
| 基于声学的迟滞效应计算内燃机的热声耦合性质;李维等;《内燃机学报》;20130725;第31卷(第4期);第373-378页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104679962A (en) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bahri et al. | Real-time modeling of ringing in HCCI engines using artificial neural networks | |
| Pan et al. | LES analysis for auto-ignition induced abnormal combustion based on a downsized SI engine | |
| Liu et al. | Machine learning assisted analysis of an ammonia engine performance | |
| Lafossas et al. | Development and validation of a knock model in spark ignition engines using a CFD code | |
| Pan et al. | Research on in-cylinder pressure oscillation characteristic during knocking combustion in spark-ignition engine | |
| Fontanesi et al. | LES multi-cycle analysis of a high performance GDI engine | |
| Fontanesi et al. | CFD analysis of the acoustic behavior of a centrifugal compressor for high performance engine application | |
| Ewald et al. | A level set based flamelet model for the prediction of combustion in spark ignition engines | |
| Kosmadakis et al. | CFD-based method with an improved ignition model for estimating cyclic variability in a spark-ignition engine fueled with methane | |
| CN104679962B (en) | Combustion noise prediction analysis method based on KIVA and SYSNOISE coupling internal combustion engines | |
| Och et al. | Volumetric efficiency optimization of a single-cylinder DI diesel engine using differential evolution algorithm | |
| CN104123414A (en) | Parameterization method for optimization of air intake and exhaust channel flow field and heat transfer of pneumatic engine | |
| Eckert et al. | Modeling autoignition and engine knock under spark ignition conditions | |
| Ramajo et al. | In-cylinder flow computational fluid dynamics analysis of a four-valve spark ignition engine: comparison between steady and dynamic tests | |
| Pera et al. | Exploitation of multi-cycle engine LES to introduce physical perturbations in 1D engine models for reproducing CCV | |
| Zhang et al. | A numerical study of natural gas combustion in a lean burn engine | |
| CN103970935A (en) | Method for identifying noise source in waste heat boiler | |
| Rabeti et al. | Investigation of important semi-empirical heat transfer models for a natural gas-fueled HCCI engine | |
| Zou et al. | Effect of recess shape on combustion performance and knocking characteristics for a downsized gasoline rotary engine | |
| Iacovano et al. | Analysis and simulation of non-flamelet turbulent combustion in a research optical engine | |
| CN107491574A (en) | A kind of method and device of vehicle motor air intake duct simulation calculation | |
| Liu et al. | Investigation of multistage combustion inside a heavy-duty natural-gas spark-ignition engine using three-dimensional computational fluid dynamics simulations and the Wiebe-function combustion model | |
| Bozza et al. | The prediction of the performance and gasdynamic noise emitted by a medium-size spark-ignition engine by means of 1D and 3D analyses | |
| Kim et al. | A Simple Model of Cyclic Variation | |
| Falkenstein et al. | Experiments and large-eddy simulation for a flowbench configuration of the darmstadt optical engine geometry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180102 Termination date: 20190317 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |