• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Unsteady detonation with thermodynamic nonequilibrium effect based on the kinetic theory

    2023-09-28 06:22:10XianliSuandChuandongLin
    Communications in Theoretical Physics 2023年7期

    Xianli Su and Chuandong Lin

    Sino-French Institute of Nuclear Engineering and Technology,Sun Yat-sen University,Zhuhai 519082,China

    Abstract

    In this paper,unsteady detonation is simulated and investigated from the viewpoint of kinetic theory.The deviations of the velocity distribution function from the equilibrium state are studied in the evolution of detonation.It has been discovered that the characteristics of the deviation around the detonation wave are significantly different from those in the post-wave region.Besides,the kinetic moments of the reaction term have been simulated,verified and analyzed in detail.In addition,the reaction manifestation is defined to describe the global effects of kinetic moments due to chemical reactions.It is interesting to find that there are three types of periodic oscillations of the reaction manifestation during the evolution of the unsteady detonation.Via the fast Fourier transform,it can be seen that the reaction manifestation is mainly composed of several signal frequencies.Moreover,the impact of rate constants of the two-step reaction scheme on the reaction manifestation is studied,and the influence of chemical heat is investigated as well.

    Keywords: unsteady detonation,chemical reaction,kinetic theory,discrete Boltzmann method

    1.Introduction

    Detonation is a kind of supersonic combustion where a chemical reaction is coupled with shock dynamics [1,2].Detonation is the key factor causing catastrophic damage in various explosion accidents,such as gas explosions [3].Besides,detonation has been widely used in aerospace,defense,mining,and demolition applications [4-6].In view of its importance,detonation has been researched by experimental [7-10],theoretical [11,12],and numerical methods[13-17].Especially,with the rapid development of computer technology and computational fluid dynamics,numerical simulation has become an indispensable method to study detonation [2,18].However,the evolutionary process of detonation involves complex physicochemical phenomena[11,12].During the evolution of the detonation,the chemical reaction and hydrodynamic processes interact with each other,which produces various interfacial and mechanical structures,and contains a wealth of hydrodynamic and thermodynamic nonequilibrium effects [19].Therefore,the accurate prediction and effective control of detonation are still open issues.

    In recent years,the kinetic methods based on the Boltzmann equation have become a popular mesoscopic simulation approach in various fields [20-23].In particular,the discrete Boltzmann method(DBM)[24,25]is suitable for supersonic waves with essential nonequilibrium effects,such as the steady or unsteady detonation.A finite number of velocities is employed to discretize the Boltzmann equation in the modeling of DBM.Thanks to its finite discrete velocities and distribution functions,the DBM owns a high computational efficiency [24,26].Moreover,the hydrodynamic and thermodynamic nonequilibrium effects beyond the macroscopic governing equations can be obtained and quantified as the higher-order kinetic moments in the DBM[27].Therefore,as a mesoscopic kinetic model,the DBM not only inherits the function of the Boltzmann equation to recover the macroscopic transport equations in the continuous limit,but also has the ability to capture the thermodynamic nonequilibrium effects beyond the macroscopic equations.The DBM has been widely utilized to obtain some deeper insights into nonequilibrium behaviours in various complex fluids and has excellent applications in the analysis of nonequilibrium,nonlinear,transient complex flow systems [22-27,29].

    Recently,the DBM has achieved remarkable success in simulating combustion and detonation [24,25,30-38].The pioneering DBM for detonation [30] was presented by Yan et al who adopted the Lee-Tarver model [39] to control the chemical reaction.In 2014,Lin et al [31] proposed a novel polar coordinate DBM and investigated typical implosion and explosion processes.In 2015,Xu et al presented a twodimensional multiple-relaxation-time model [32] and studied the hydrodynamic and thermodynamic nonequilibrium in the process of steady detonation.In 2016,a double-distributionfunction DBM for combustion and detonation was constructed by Lin et al [33],where the chemical reactant is described by one distribution function,and the chemical product by the other distribution function.Then,Lin et al reported a multi-component DBM for premixed,nonpremixed,or partially premixed nonequilibrium reactive flows [35].In 2019,Lin et al developed a multiple-relaxation-time DBM for unsteady reactive flows[24],where the specific heat ratio and Prandtl number are adjustable.In 2020,Lin et al simulated unsteady detonation with essential hydrodynamic and thermodynamic nonequilibrium effects,and studied the impact of the perturbation amplitude,wavelength,and chemical heat release on the physical field of unsteady detonation with nonequilibrium effects[25].Next,a three-dimensional DBM for the simulation of detonation was developed by Ji et al[36],and they captured the features of the three-dimensional detonation via the DBM[37].Very recently,Su et al recovered the particle velocity distribution function of reactive flows from the gas kinetic theory and investigated the higher-order kinetic moments and nonequilibrium quantities around the detonation wave by using the DBM [38].

    In the present study,we carry out more in-depth and comprehensive research on unsteady detonation via the DBM.In section 2,the details of the modeling construction are presented.Then,the DBM is employed to simulate the unsteady detonation.The deviation of the velocity distribution function from the equilibrium state and the kinetic moments of the reaction term are studied in section 3.Reaction manifestations in moment space are investigated from a kinetic perspective in section 4.Finally,section 5 provides a summary.

    2.Discrete Boltzmann model

    The DBM is based on the following equation,

    where τ denotes the relaxation time,t the time,fiandfieqrepresent the discrete distribution function and its equilibrium counterpart,respectively.Moreover,videnotes the discrete velocity with i=1,2,3,…,N,and N=16 is the total number of discrete velocities.Here a two-dimensional sixteen-velocity model is utilized,see figure 1.

    In especial,the original expression of equilibrium distribution function feqis as follows [24]

    where n,u,and T represent the particle number density,flow velocity and temperature,respectively.Besides,D=2 denotes the dimension,I counts extra degrees of freedom due to vibration and/or rotation,and η is used to describe the corresponding vibrational and/or rotational energies.

    On the right-hand side of equation (1),Riindicates the reaction term describing the change rate of the distribution function due to chemical reactions.In this paper,a two-step reaction scheme is employed to control the chemical reaction process [40],

    The first step in equation (3) denotes a thermally neutral induction zone or ignition process,and the second step in equation (4) describes the rapid energy release after the branched-chain thermal explosion.ξ′ denotes the change rate of the reaction progress variable in a thermally neutral induction period,whileλ′ is the change rate of the mass fraction of the chemical product.Besides,Tsis the temperature after the preshock wave.EIand ERrepresent the activation energies.KIand KRindicate the rate constants.

    Under the assumption that the time scale of molecular collisions < the time scale of chemical reactions < the characteristic time scale of flows,the original reaction term can be obtained [38],

    where Q indicates the chemical heat release per unit mass of chemical reactant.

    In this work,in order to recover the Navier-Stokes (NS)equations,the following moment relations should be satisfied,

    with the particle velocities Ψ=1,v,(v·v+η2),vv,(v·v+η2)v,vvv,(v·v+η2)vv,and the corresponding discrete velocities,Ψi=1,vi,(vi·vi+η2i),vivi,(vi·vi+ηi2)vi,vivivi,(vi·vi+ηi2)vivi.Substituting equations (2) into (6) leads to an explicit expression:

    where M is a 16×16 square matrix bridging the velocity and moment space,consisting of the functions of the discrete velocities.In addition,feqandHfeqare matrix expressions of the equilibrium distribution function in velocity and moment space,respectively.The details are in appendix A.

    Similarly,the reaction terms follow the relations,

    The elements of Ψ and Ψiare the same as those in equation (6).Substituting equations (5) into (8) results in,HR=MR,where R and HRcorrespond to the matrices of reaction term in velocity and moment space,respectively.Hence,

    where M?1is the inverse of M.Mathematically,the matrix HRis expressed byHR=(HR1HR2…HR16)T.Physically,HR1represents the change rate of mass density due to chemical reactions.HR2and HR3describe the corresponding change rate of momentum in the x and y directions,respectively.HR4indicates the change rate of energy because of the chemical reactions.HR5,HR6and HR7describe the change rate of nonorganised momentum fluxes due to chemical reactions.HR8,HR9,HR10,HR11,HR12,and HR13are related to the change rate of nonorganised energy fluxes due to chemical reactions.HR14,HR15,and HR16correspond to the change rate of the fluxes of nonorganised energy fluxes due to chemical reactions.

    3.Investigation of unsteady detonation

    In this section,the DBM is employed to simulate the onedimensional unsteady detonation.The initial configuration,obtained from the Hugoniot relation,takes the form

    where the subscript L indicates 0 ≤x ≤0.01,and R indicates 0.01

    The pressure,temperature,density and horizontal velocity are illustrated in figure 2(a).Obviously,as the detonation wave propagates from left to right,the macroscopic physical quantities in the medium first rise drastically and then fall sharply to form a high peak.Then there are regular fluctuations after the peak.The three dashed lines in figure 2(a)denote the locations x=1.1724,x=1.1799,and x=1.2518,respectively.

    Figure 2(b) exhibits the nonorganized energy in the x direction.The inset shows the corresponding profiles within the spatial range 1.245

    Next,let us investigate the velocity distribution function of unsteady detonation.In [38],via the Chapman-Enskog analysis,the local particle velocity distribution function of reactive flow is derived from the gas kinetic theory.In short,on the NS level,the local particle velocity distribution function reads

    where feqis the equilibrium particle velocity distribution function,and fneqis the deviation of the velocity distribution function from the equilibrium state.Besides,fneqcan be written as [38],

    From equations(10)and(11),it can be seen that both spatial gradients and reaction term make contributions to the deviation of the local velocity distribution function from the equilibrium state,while the velocity distribution function is composed of the equilibrium particle velocity distribution function and the corresponding deviation counterpart.

    Next,we conduct an in-depth study of fneq.Since the macroscopic quantities show regular oscillations in figure 2,we first investigate the deviations at a trough (x=1.1724) and a peak (x=1.1799) of the pressure in the region after the detonation front,and the position of the highest peak(x=1.2518)of pressure around the detonation wave is also researched.To this end,figure 3 depicts the deviation of the particle velocity distribution function from the equilibrium statefneq=fneq(vx,vy)in the left column and corresponding contours in the velocity space(vx,vy)in the right column.To be specifci,fgiures 3(a)and(b)exhibit the deviation fneqand its corresponding contour at the valley of the pressure.Figures 3(c) and (d) display the deviation fneqand the relevant contour at the peak of the pressure.Besides,the deviation fneqand its contour at the position of the highest peak of pressure around the detonation wave are drawn in figures 3(e) and (f),respectively.

    As shown in figures 3(a) and (c),there are both positive and negative deviations.Figure 3(a) depicts that along the vxdirection,a deep negative valley first appears,and then there is a positive peak,afterwards it decreases to form a shallow negative valley.On the contrary,there is an opposite trend in figure 3(c).Firstly,a small positive peak appears,then drops to a negative trough,and later it increases to a high positive peak.Meanwhile,it is interesting to note that the curves of figures 3(b) and (d) look quite similar to each other.Both structures are composed of three parts,and the corresponding shapes of the three parts are the same.However,the difference between figures 3(b)and(d)mainly lies in the values of each contour line,especially whether the values are positive or negative.Specifically,from left to right,the values of three parts of figure 3(b) are first negative,then positive,and finally negative,which is consistent with the trend of figure 3(a).Similarly,the values of the three parts in figure 3(d) are consistent with the trend of figure 3(c),which contains positive,negative,and positive parts.Moreover,From the three-dimensional view,figures 3(e) looks like a bulging positive circle surrounding a negative valley.Besides,it can be found from the contour that figure 3(f) is symmetric in the vxand vydirections.

    In sum,it can be found in figure 3 that the characteristics of the deviation near the wavefront of the detonation are significantly different from those in the post-wave region.Besides,the nonequilibrium intensity around the detonation wave is much greater than that in the post-wave region.Moreover,in the post-wave region of the unsteady detonation wave,the deviation of the particle velocity distribution function from the equilibrium state changes with the macroscopic quantities and nonorganized energy.

    It should be pointed out that the DBM is a mesoscopic kinetic method,so the chemical reaction of unsteady detonation can be studied from the kinetic moment viewpoint.Figure 4 depicts the sixteen moments of reaction term HR.

    Figure 4(a) exhibits the profiles of moments of reaction term HR1to HR4.The symbols denote the DBM results,and the solid lines represent the theoretical solutions.Clearly,the simulation results are in nice agreement with the theoretical solutions.Besides,it can be seen that during the evolution of unsteady detonation wave,HR1,HR2and HR3keep zero due to the conservation of mass and moment.HR4first increases sharply and then decreases smoothly,and forms a towering asymmetric peak.The reason is that there is violent chemical heat released around the detonation wave.

    Figure 4(b) displays the moments of reaction term HR5,HR6and HR7.It is obvious that HR6is zero,while HR5and HR7first grow and then decline,forming a positive crest.Besides,the trends of HR5and HR7are exactly the same,which is also consistent with the theoretical solutionsHR5=ρT′ andHR7=ρT′.

    The patterns of moments of reaction term HR8,HR9,HR10,HR11,HR12,and HR13are shown in figure 4(c).Clearly,HR9,HR11,and HR13keep zero,while HR8,HR10,and HR12all increase first and then drop sharply.Moreover,the peak of HR8is higher and wider than HR10,and the peak of HR12is thinner and lower than HR10.The evolution in figure 4(c) is in accordance with the theoretical solutionsHR8=(D+I+2)ρuxT′,HR9=(D+I+2)ρuyT′,HR10=3ρuxT′,HR11=ρuyT′,HR12=ρuxT′,andHR13=3ρuyT′,respectively.

    The rest change rate of the fluxes of nonorganized energy flux due to chemical reactions HR14to HR16is plotted in figure 4(d).Similarly,HR15keeps zero as the detonation wave propagates from left to right,while HR14and HR16show positive peaks.In addition,the peak of HR14is taller and wider than HR16.The theoretical solutionsHR14=ρ[ 2T(D+I+2)+(D+I+5)ux2+uy2]T′,HR15=ρuxuy(D+I+4)T′,andHR16=ρ[2T(D+I+2)+ux2+(D+I+5)uy2]T′ explain the variation patterns of figure 4(d).

    4.Reaction manifestation in moment space

    In order to have a better understanding of chemical reactions of unsteady detonation from the perspective of kinetic moments,the reaction manifestation GRis defined to describe the global effects due to the chemical reaction in moment space

    where the integral is over the whole computational domain.It should be mentioned that it is a coarse-grained way to investigate physical characteristics and reaction manifestations by using equation(12),which denotes the magnitude of the quantity in the moment space(HR1,HR2,...,HR16).That is to say,the mathematical expression of equation (12) consists of 16 independent variables.

    Figure 5 shows the evolution of the reaction manifestation in moment space.Figure 5(a)shows the whole process of reaction manifestation,it can be seen that when t<0.1,the reaction manifestation changes very drastically,in especial,an extremely high peak was observed around t=0.04.After t>0.15,the reaction manifestation begins to show regular oscillations.In figure 5(b),the red dashed lines represent the envelope curves.From the upper envelop curve,it can be seen that the maximum of reaction manifestation displays a periodic small change when t>0.15.One cycle of the first type of periodic oscillation is shown by the red rectangular frame,and the first cycle is about T1=0.040 28.Then,we further narrow the time scale of the observed reaction manifestation.Figure 5(c)illustrates the second type of cycle of the reaction manifestation around T2=0.006 69,i.e.the pattern of change from a valley to a peak,and one of the cycles is displayed by the green rectangular frame.Besides,after partially zooming in figure 5(c),we see that the data points which are plotted by hollow squares also seem to evolve with a much smaller period from figure 5(d),i.e.they alternate in groups of about seven points.Similarly,one cycle of the third period is about T3=0.000 012 signed by the black rectangular frame.

    As shown in figure 5,there are three types of periodic changes in the reaction manifestation in the moment space of unsteady detonation during the post-stabilization evolution.Therefore,the fast Fourier transform is introduced to deal with the reaction manifestation in figure 5(b).Figure 6 illustrates the amplitude versus frequency f(HZ).It is apparent that there are many peaks in various frequencies.As the frequency increases,the amplitude decreases gradually.The first eight peaks are more distinct,which denotes that the reaction manifestation is composed mainly of the first several groups' signal of different frequencies.

    In addition,the specific frequency and corresponding amplitude are shown in table 1.The first line denotes the number of the first eight frequencies,and the second and third lines give the values of the corresponding frequency and amplitude,respectively.It can be found from the table that the frequencies are all roughly multiples of f=150 (HZ).

    The data in table 1 are fitted to further explore the relations among the number,frequency,and amplitude.Figure 7(a)depicts the values of frequency versus the number of frequency.The squares denote the numerical data,and the solid line represents the fitting function f=150n,where f and n stand for the frequency and number,respectively.It is evident that the original data agree well with the fitting solutions in figure 7(a),which means that the reaction manifestation may consist of a superposition of frequencies of different multiples of f=150 (HZ).Figure 7(b) displays the amplitude (A) versus frequency.In the same way,the circles represent the numerical data,and the line is for the fitting functionA=0.01637 ×exp(-f/210.9455).From the ftiting result,it can be seen that there is an exponentially decreasing relationship between the frequency and amplitude.The amplitude decreases exponentially as the frequency increases.

    Next,let us study the effects of different physical factors on the reaction manifestation in the detonation process.We focus on the changes in the maximum and amplitude of the reaction manifestation.In fact,the bifurcation boundary between steady and unsteady detonations depends on the ratio of the length of the heat release layer to that of the induction zone layer [40,41].Generally,the length of the induction zone can be adjusted by adjusting KI,meanwhile,the length of the heat release layer can be determined by KR.The two rate constants play an important role in the evolution of the detonation.Therefore,we first consider unsteady detonations with the rate constant KItunable,and the other parameters fixed the same as the above simulation.

    Figures 8(a)and(b)exhibit the maximum and amplitude of the reaction manifestation versus the rate constant KI,respectively.The symbols denote the simulation results and the lines are for the fitting results.The fitting equations in figures 8(a)and (b) areYMax=0.11305-2.2505 ×10-4×KIand A=0.10767 ?2.41021×10?4×KI,respectively.Obviously,with the increasing rate constant KI,the maximum and amplitude of the reaction manifestation decrease.Figure 9 illustrates the physical quantities around the detonation wave with two different rate constants KI.As shown in figures 9(a)-(d),with the increasing KI,the amplitudes of the oscillation of density,pressure,horizontal velocity and temperature around the detonation wave become smaller,which means that the chemical reactions become less violent.In other words,the reaction manifestation becomes weaker.

    Table 1.The frequencies and corresponding amplitudes of the reaction manifestation.

    Next,we simulate the unsteady detonation wave with various values of rate constant KRfrom KR=7000 to KR=13000.Figure 10 illustrates the simulation results and fitting solutions which are displayed by the symbols and solid lines,respectively.The fitting equations areYMax=0.02175+2.3068 ×10-6×KRand A=0.00818+2.64065×10?6×KR,respectively.Obviously,the values of maximum and amplitude grow linearly with the increasing KR.Figures 11(a)-(d) display the evolutions of density,pressure,horizontal velocity and temperature with different rate constants KR,respectively.It is clear in figure 11 that with the increase of KR,the macroscopic physical quantities change more violently,which means that the intensity of chemical reactions is enhanced.Consequently,the global effect caused by chemical reactions on the moment space is enhanced,and the reaction manifestation becomes stronger.

    Figure 1.Sketch of the two-dimensional sixteen-velocity model.

    Figure 2.Profiles of physical quantities (a) and nonorganized energy (b) of the unsteady detonation wave.

    Figure 4.Kinetic moments of reaction term around the detonation wave.(a)HR1 to HR4,(b)HR5 to HR7,(c)HR8 to HR13,and(d)HR14 to HR16.The symbols show the numerical results,and the solid lines represent the theoretical solutions.

    Figure 5.Evolution of the reaction manifestation in moment space.(a) 0 ≤t ≤0.4,(b) 0.15 ≤t ≤0.35,(c) 0.24 ≤t ≤0.26,and (d)0.25 ≤t ≤0.2501.

    Figure 6.The amplitude versus frequency from the fast Fourier transform performed on the reaction manifestation.

    Figure 7.The fitting results of the reaction manifestation.(a) the relation between the number and frequency,and (b) amplitude versus frequency.

    Figure 8.The maximum (a) and amplitude (b) of the reaction manifestation versus rate constant KI.

    Figure 9.The physical quantities with different rate constants KI,(a) density,(b) pressure,(c) horizontal velocity,and (d) temperature.

    Figure 10.The maximum (a) and amplitude (b) of the reaction manifestation versus rate constant KR.

    Figure 11.The physical quantities with different rate constants KR,(a) density,(b) pressure,(c) horizontal velocity,and (d) temperature.

    Finally,the simulations of unsteady detonation waves are carried out for various values of chemical heat release per unit mass from Q=3 to 8.Figures 12(a) and (b) show the simulation and fitting results.Symbols represent the simulation results,and the solid lines denote the fitting equationsYMax=?0.12108+0.04272×Q and A=?0.10301+0.03554×Q,respectively.It can be found that as the chemical heat increases,the values of maximum and amplitude increase in a linear form.Physically,for larger chemical heat,the chemical reaction occurs more violently,the physical quantities change more sharply around the detonation wave,and the reaction manifestation becomes larger.

    5.Conclusions

    In this work,we have investigated the characteristics of unsteady detonation from the kinetic viewpoint.The unsteady detonations are simulated via the DBM,which has the capability of capturing detonation waves with nonequilibrium effects.As the detonation wave propagates forwards,the macroscopic physical quantities in the medium first rise drastically and then fall sharply to form a peak.There are regular oscillations after the detonation peak.Besides,the deviation of the velocity distribution function from the local equilibrium state of unsteady detonation oscillates with the macroscopic quantities.Especially,the deviations show opposite patterns at the peak and trough of pressure after the detonation wave.Moreover,the characteristics of the deviation at the peak of pressure are significantly different from those in the post-wave region.In addition,the nonequilibrium intensity around the detonation wave is much greater than that in the post-wave region.Furthermore,the DBM is adopted to investigate the chemical reaction of unsteady detonation,and the change rate of the kinetic moments due to chemical reactions is consistent with their theoretical solution.

    In order to have a deeper understanding of chemical reactions from the kinetic perspective,the reaction manifestation is defined to describe the global effects.There are three types of periodic oscillations in the reaction manifestation in moment space during the evolution of the unsteady detonation,and the corresponding cycles are about T1=0.04 028,T2=0.006 69 and T3=0.000 012,respectively.Via the fast Fourier transform,it can be seen that the reaction manifestation is mainly composed of several signal frequencies.Besides,it has been demonstrated that the reaction manifestation consists of a superposition of frequencies of different multiples of f=150 (HZ),and the corresponding amplitude decreases exponentially as the frequency increases.Finally,the influences of rate constants and the chemical heat on the reaction manifestation have been investigated.From the numerical simulations and corresponding fitting equations,the following points can be obtained: (I) with the increasing rate constant KI,the maximum and amplitude of the reaction manifestation decrease.(II) The values of maximum and amplitude of the reaction manifestation grow linearly with the increasing rate constant KR.(III) As the chemical heat increases,the values of maximum and amplitude increase in a linear form.

    Acknowledgments

    This work is supported by the National Natural Science Foundation of China (under Grant No.51806116),and Guangdong Basic and Applied Basic Research Foundation(under Grant No.2022A1515012116).

    Appendix A

    The matrixHfeqtakes the form,

    The matrix HRis expressed by,

    whose elements are HR1=0,HR2=0,HR3=0,HR4=ρ(D+I)T′,HR5=ρT′,HR6=0,HR7=ρT′,HR8=(D+I+2)ρuxT′,HR9=(D+I+2)ρuyT′,HR10=3ρuxT′,HR11=ρuyT′,HR12=ρuxT′,HR13=3ρuyT′,HR14=ρ[ 2T(D+I+2)+(D+I+5)ux2+uy2]T′,HR15=ρuxuy(D+I+4)T′,andHR16=ρ[ 2T(D+I+2)+ux2+(D+I+5)uy2]T′.

    The matrix M is specified as

    containing the blocks,Mi=(Mi1Mi2…Mi16),with elementsM1i=1,M2i=vix,M3i=viy,M4i=vi2+ηi2,M5i=vi2x,M6i=vixviy,M7i=v2iy,M8i=

    Appendix B

    Most of the previous research on the mechanism of unstable detonation is studied using somehow CFD method,which is quite different from DBM.Roughly speaking,a DBM is approximately equivalent to a continuous fluid model plus a coarse-grained model of other relevant thermodynamic nonequilibrium effects.Therefore,DBM has the ability to capture unsteady detonation with nonequilibrium effects,including diffusion and thermal conduction.In this sector,the DBM results are compared with a CFD simulation.The same initial configurations are used in the CFD simulation.

    Figure 13 depicts the pressure profiles of unsteady detonation at various times.The red lines represent the DBM results,and the blue lines plot the CFD solutions via an Euler solver.It is clear that the profiles simulated by the DBM are smoother than those using the Euler solver,and the amplitudes simulated by DBM are smaller compared to Euler solver solutions.Physically,the increasing diffusion smooths the density gradient,and the increasing thermal conduction smooths the temperature gradient.Consequently,the pressure amplitude is smaller when considering nonequilibrium effects,including the diffusion and thermal conduction [28].The simulation via Euler solver lacks essential thermodynamics of nonequilibrium effects.In addition,the detonation velocities calculated by the DBM and Euler solver are 3.166 and 3.178,respectively.The corresponding theoretical solution is 3.208.Therefore,the results are generally satisfactory.

    Figure 13.Comparisons between Euler solver solutions and DBM results of the unsteady detonation at t=0.36(a),t=0.38(b),and t=0.4(c),respectively.

    国产乱人伦免费视频| 9191精品国产免费久久| 在线观看免费午夜福利视频| 制服丝袜大香蕉在线| 国产欧美日韩一区二区三| 午夜福利在线观看免费完整高清在 | 最近最新中文字幕大全电影3| 亚洲欧美日韩高清专用| 国产精品女同一区二区软件 | 国产精品久久视频播放| 国产伦精品一区二区三区视频9 | 久久精品91蜜桃| 精品久久久久久久人妻蜜臀av| 两人在一起打扑克的视频| 精品久久久久久久末码| 欧美日韩中文字幕国产精品一区二区三区| 18禁黄网站禁片免费观看直播| 日韩中文字幕欧美一区二区| 91九色精品人成在线观看| 国产精品永久免费网站| 热99re8久久精品国产| 欧美激情在线99| 黄片小视频在线播放| 精品人妻1区二区| 手机成人av网站| 美女被艹到高潮喷水动态| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区精品| 51午夜福利影视在线观看| 亚洲精品色激情综合| 国产男靠女视频免费网站| 欧美xxxx黑人xx丫x性爽| 亚洲精品粉嫩美女一区| 国产精品香港三级国产av潘金莲| 我的老师免费观看完整版| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人a区在线观看| 国产v大片淫在线免费观看| 成年人黄色毛片网站| 亚洲av第一区精品v没综合| 九色成人免费人妻av| 国产精品香港三级国产av潘金莲| 国产精品综合久久久久久久免费| 久久久久国产精品人妻aⅴ院| 国产av麻豆久久久久久久| 亚洲人成网站在线播| 美女 人体艺术 gogo| 精品久久久久久久末码| 欧美性感艳星| 波多野结衣高清作品| 最新中文字幕久久久久| 日韩免费av在线播放| 欧美最新免费一区二区三区 | 美女 人体艺术 gogo| 国产成人影院久久av| 我要搜黄色片| 亚洲精品国产精品久久久不卡| 久久精品国产亚洲av香蕉五月| 久久人人精品亚洲av| 两个人的视频大全免费| 亚洲色图av天堂| 亚洲成av人片在线播放无| 国产精品一区二区免费欧美| 精品日产1卡2卡| aaaaa片日本免费| 久久午夜亚洲精品久久| 丝袜美腿在线中文| 欧美性猛交黑人性爽| 97超视频在线观看视频| eeuss影院久久| 久久精品人妻少妇| 亚洲中文字幕日韩| 国内精品一区二区在线观看| 麻豆一二三区av精品| 男女那种视频在线观看| 色综合欧美亚洲国产小说| 精品无人区乱码1区二区| 给我免费播放毛片高清在线观看| 国产精品久久久久久人妻精品电影| 欧美性猛交╳xxx乱大交人| 免费高清视频大片| 一进一出好大好爽视频| 亚洲成人中文字幕在线播放| 好男人电影高清在线观看| 国产精品久久久久久久电影 | 男女之事视频高清在线观看| 麻豆国产97在线/欧美| 精品乱码久久久久久99久播| 在线观看一区二区三区| 国产男靠女视频免费网站| 精品电影一区二区在线| 精品国产亚洲在线| 亚洲精品一卡2卡三卡4卡5卡| 国产探花极品一区二区| 亚洲国产日韩欧美精品在线观看 | 免费电影在线观看免费观看| 精品一区二区三区人妻视频| 亚洲中文日韩欧美视频| 久久国产精品人妻蜜桃| 亚洲人与动物交配视频| 亚洲av一区综合| 此物有八面人人有两片| 精品国产超薄肉色丝袜足j| 国产精品1区2区在线观看.| 无限看片的www在线观看| 看片在线看免费视频| 亚洲欧美日韩高清在线视频| 午夜久久久久精精品| 熟女少妇亚洲综合色aaa.| 蜜桃亚洲精品一区二区三区| 美女 人体艺术 gogo| 日韩有码中文字幕| tocl精华| 禁无遮挡网站| 免费看十八禁软件| 色精品久久人妻99蜜桃| 俺也久久电影网| 一区二区三区免费毛片| 国产精品精品国产色婷婷| 亚洲精品色激情综合| 午夜福利欧美成人| 国产不卡一卡二| 夜夜看夜夜爽夜夜摸| 少妇高潮的动态图| 黄色丝袜av网址大全| 国产中年淑女户外野战色| 男女视频在线观看网站免费| 床上黄色一级片| 亚洲美女视频黄频| 久久中文看片网| 久久久久久久久大av| 日韩 欧美 亚洲 中文字幕| ponron亚洲| 国产午夜精品论理片| 国产精品 欧美亚洲| av女优亚洲男人天堂| 国产欧美日韩一区二区精品| 制服丝袜大香蕉在线| 精品国产美女av久久久久小说| 精品人妻偷拍中文字幕| 黄色日韩在线| 18+在线观看网站| 亚洲天堂国产精品一区在线| 久久精品国产99精品国产亚洲性色| 最新在线观看一区二区三区| 在线看三级毛片| 欧美3d第一页| 欧美中文日本在线观看视频| 国产97色在线日韩免费| e午夜精品久久久久久久| 在线免费观看不下载黄p国产 | 两个人的视频大全免费| 人妻夜夜爽99麻豆av| av在线蜜桃| 午夜福利高清视频| 亚洲 国产 在线| 精品人妻1区二区| 精品久久久久久,| 美女免费视频网站| 精华霜和精华液先用哪个| 美女免费视频网站| 一区福利在线观看| 国产蜜桃级精品一区二区三区| 国产精品久久久人人做人人爽| 一进一出好大好爽视频| 国产精品免费一区二区三区在线| 日韩欧美三级三区| 变态另类丝袜制服| 午夜福利免费观看在线| 听说在线观看完整版免费高清| 深爱激情五月婷婷| 欧美成人性av电影在线观看| 69av精品久久久久久| 欧美黑人巨大hd| 男人的好看免费观看在线视频| 国产亚洲av嫩草精品影院| 99精品久久久久人妻精品| bbb黄色大片| 国产亚洲欧美98| 人妻久久中文字幕网| 久久久久精品国产欧美久久久| 亚洲人成网站高清观看| 日日夜夜操网爽| 国产亚洲精品av在线| 亚洲va日本ⅴa欧美va伊人久久| 白带黄色成豆腐渣| 小说图片视频综合网站| 国内精品美女久久久久久| 国产爱豆传媒在线观看| 国产av在哪里看| 一级a爱片免费观看的视频| 男女下面进入的视频免费午夜| 国产精品一及| 少妇的逼好多水| 极品教师在线免费播放| 男女之事视频高清在线观看| 夜夜爽天天搞| АⅤ资源中文在线天堂| 中文字幕熟女人妻在线| 欧美激情在线99| 给我免费播放毛片高清在线观看| 51国产日韩欧美| 国产视频内射| 成人18禁在线播放| 免费人成在线观看视频色| 国产亚洲精品一区二区www| 国产乱人伦免费视频| 国产一区二区三区视频了| 精品人妻1区二区| 国产主播在线观看一区二区| 国产成年人精品一区二区| 最近视频中文字幕2019在线8| 99国产精品一区二区蜜桃av| 俺也久久电影网| 国内精品久久久久久久电影| 精品久久久久久成人av| 一级毛片女人18水好多| 国产精品99久久久久久久久| 女同久久另类99精品国产91| 亚洲 国产 在线| 欧美bdsm另类| 99视频精品全部免费 在线| 国产视频一区二区在线看| 男人舔奶头视频| 男人舔女人下体高潮全视频| 免费av毛片视频| 亚洲av电影在线进入| 69av精品久久久久久| 久9热在线精品视频| www国产在线视频色| 制服人妻中文乱码| 9191精品国产免费久久| 在线观看66精品国产| 中文字幕av成人在线电影| 一a级毛片在线观看| 变态另类丝袜制服| 欧美激情在线99| 亚洲人成电影免费在线| 成人亚洲精品av一区二区| 看免费av毛片| 欧美日韩瑟瑟在线播放| 1000部很黄的大片| 亚洲成人中文字幕在线播放| 天堂动漫精品| 最近最新中文字幕大全免费视频| 国产成人影院久久av| 日本黄色片子视频| xxx96com| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲 欧美 日韩 在线 免费| 日韩欧美精品v在线| 天堂动漫精品| 成人午夜高清在线视频| 成人性生交大片免费视频hd| 亚洲avbb在线观看| 免费一级毛片在线播放高清视频| 亚洲色图av天堂| 欧美又色又爽又黄视频| 99久久精品热视频| 国产精品自产拍在线观看55亚洲| 日韩欧美精品免费久久 | 狂野欧美白嫩少妇大欣赏| 色尼玛亚洲综合影院| 搡老熟女国产l中国老女人| 精品国产亚洲在线| 欧美在线黄色| 亚洲av日韩精品久久久久久密| 免费看十八禁软件| 亚洲一区高清亚洲精品| 免费在线观看日本一区| 欧美日韩精品网址| 婷婷精品国产亚洲av在线| 美女黄网站色视频| 两个人视频免费观看高清| 无人区码免费观看不卡| 国产在视频线在精品| 久久精品国产亚洲av涩爱 | e午夜精品久久久久久久| 久久久久久久亚洲中文字幕 | 色老头精品视频在线观看| 97碰自拍视频| 国产成年人精品一区二区| 色播亚洲综合网| 国产69精品久久久久777片| 免费人成在线观看视频色| 女同久久另类99精品国产91| 一级a爱片免费观看的视频| 久久精品综合一区二区三区| 国模一区二区三区四区视频| 成年女人毛片免费观看观看9| 国产探花在线观看一区二区| 久久国产精品影院| 精华霜和精华液先用哪个| 噜噜噜噜噜久久久久久91| 国产欧美日韩一区二区精品| www.熟女人妻精品国产| 日韩国内少妇激情av| 人人妻,人人澡人人爽秒播| 天堂动漫精品| 成人性生交大片免费视频hd| 美女高潮的动态| 午夜精品一区二区三区免费看| 国产毛片a区久久久久| 1000部很黄的大片| 最新中文字幕久久久久| 国产精品三级大全| 午夜福利在线观看吧| 亚洲成a人片在线一区二区| 久久伊人香网站| 国产探花在线观看一区二区| 国产伦在线观看视频一区| 一级黄色大片毛片| 国产精品久久久久久亚洲av鲁大| 精品国产三级普通话版| 亚洲aⅴ乱码一区二区在线播放| 国产97色在线日韩免费| 99精品在免费线老司机午夜| 国产精品亚洲一级av第二区| 亚洲 国产 在线| 国产一区在线观看成人免费| 一进一出好大好爽视频| 90打野战视频偷拍视频| 国产精品久久久久久人妻精品电影| 很黄的视频免费| 国产av一区在线观看免费| 天堂av国产一区二区熟女人妻| 激情在线观看视频在线高清| 国产精品综合久久久久久久免费| 欧美bdsm另类| 亚洲无线在线观看| 国模一区二区三区四区视频| 少妇的丰满在线观看| 小蜜桃在线观看免费完整版高清| 亚洲成人久久性| 亚洲国产精品sss在线观看| 嫩草影院精品99| 一进一出抽搐动态| 99精品在免费线老司机午夜| 动漫黄色视频在线观看| 色吧在线观看| 最后的刺客免费高清国语| 欧美中文综合在线视频| 婷婷亚洲欧美| www日本黄色视频网| 精品免费久久久久久久清纯| 国产成人欧美在线观看| 国产成+人综合+亚洲专区| 国产野战对白在线观看| 露出奶头的视频| 免费在线观看亚洲国产| 日韩欧美在线乱码| 午夜福利在线在线| 中文字幕熟女人妻在线| 国产视频内射| 久久精品亚洲精品国产色婷小说| 国产av在哪里看| 国产精品女同一区二区软件 | 国产真实伦视频高清在线观看 | 成人国产综合亚洲| 日韩 欧美 亚洲 中文字幕| 成人午夜高清在线视频| 校园春色视频在线观看| 成人国产综合亚洲| 草草在线视频免费看| 成人无遮挡网站| 最新美女视频免费是黄的| 亚洲国产高清在线一区二区三| 亚洲,欧美精品.| 精品欧美国产一区二区三| 免费人成在线观看视频色| 国内少妇人妻偷人精品xxx网站| 麻豆国产97在线/欧美| bbb黄色大片| 国产亚洲精品av在线| 亚洲自拍偷在线| 久久久久久久精品吃奶| svipshipincom国产片| 九九在线视频观看精品| 国产高清视频在线观看网站| 亚洲最大成人手机在线| 成人一区二区视频在线观看| 国产乱人视频| 欧美性猛交╳xxx乱大交人| 成年女人毛片免费观看观看9| 久久久久国产精品人妻aⅴ院| 日本 av在线| 一本久久中文字幕| 脱女人内裤的视频| 日韩欧美国产一区二区入口| 亚洲性夜色夜夜综合| 99久国产av精品| 国产精品1区2区在线观看.| 97超视频在线观看视频| 成人av在线播放网站| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 久久6这里有精品| www.熟女人妻精品国产| 午夜福利在线观看免费完整高清在 | 老司机午夜十八禁免费视频| 熟妇人妻久久中文字幕3abv| 国产主播在线观看一区二区| 黄片大片在线免费观看| 怎么达到女性高潮| 国产精品影院久久| 国产高清有码在线观看视频| www.熟女人妻精品国产| 三级国产精品欧美在线观看| 国产亚洲欧美在线一区二区| 色视频www国产| 国产精品99久久99久久久不卡| 97超视频在线观看视频| 18禁黄网站禁片午夜丰满| 又黄又粗又硬又大视频| 黄色女人牲交| 国产淫片久久久久久久久 | 欧美丝袜亚洲另类 | 久久亚洲真实| 亚洲在线观看片| 久久久久亚洲av毛片大全| 国产高清有码在线观看视频| 制服人妻中文乱码| 亚洲 国产 在线| 国产高清有码在线观看视频| 久久国产精品人妻蜜桃| 日韩欧美 国产精品| 久久99热这里只有精品18| 18禁裸乳无遮挡免费网站照片| 精品一区二区三区av网在线观看| 伊人久久大香线蕉亚洲五| 色吧在线观看| 亚洲五月婷婷丁香| 两个人看的免费小视频| 国内精品一区二区在线观看| 国产精品久久久久久精品电影| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲美女久久久| 一进一出好大好爽视频| 亚洲人与动物交配视频| 欧美日韩瑟瑟在线播放| 男女床上黄色一级片免费看| 日本撒尿小便嘘嘘汇集6| 成年版毛片免费区| 中文字幕av成人在线电影| 亚洲av第一区精品v没综合| 最近最新中文字幕大全电影3| 久久精品人妻少妇| 91久久精品国产一区二区成人 | 在线天堂最新版资源| 日韩精品中文字幕看吧| 一进一出好大好爽视频| 美女 人体艺术 gogo| 网址你懂的国产日韩在线| 国产av不卡久久| 国产精品久久久久久久久免 | 国产综合懂色| av天堂在线播放| 久久久久国内视频| 18+在线观看网站| 国产精品久久久人人做人人爽| 国产麻豆成人av免费视频| 精品福利观看| 亚洲国产欧美人成| 内地一区二区视频在线| 亚洲专区中文字幕在线| 老鸭窝网址在线观看| 一级毛片高清免费大全| 热99在线观看视频| 黄色片一级片一级黄色片| 国产精品98久久久久久宅男小说| 国产国拍精品亚洲av在线观看 | 日韩国内少妇激情av| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 国产高潮美女av| 色播亚洲综合网| 午夜视频国产福利| 黄色女人牲交| 夜夜夜夜夜久久久久| 亚洲精品在线美女| 性色av乱码一区二区三区2| 一级毛片女人18水好多| 一a级毛片在线观看| 51国产日韩欧美| 精品不卡国产一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 一级作爱视频免费观看| 黄色成人免费大全| 欧美一区二区亚洲| 日韩大尺度精品在线看网址| 一a级毛片在线观看| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 91在线精品国自产拍蜜月 | 亚洲国产欧洲综合997久久,| 午夜福利成人在线免费观看| 成人国产一区最新在线观看| 18美女黄网站色大片免费观看| 欧美乱码精品一区二区三区| 黄片小视频在线播放| 久久草成人影院| 日韩精品青青久久久久久| 久久性视频一级片| 欧美不卡视频在线免费观看| 乱人视频在线观看| 国内精品久久久久精免费| 婷婷六月久久综合丁香| 欧美黄色淫秽网站| 精品电影一区二区在线| 嫁个100分男人电影在线观看| 757午夜福利合集在线观看| 我要搜黄色片| 窝窝影院91人妻| 久久精品国产自在天天线| 女警被强在线播放| 亚洲中文字幕日韩| 国产高清有码在线观看视频| av专区在线播放| 国产精品久久久人人做人人爽| 在线看三级毛片| 欧美日韩黄片免| 91久久精品国产一区二区成人 | 欧美日韩黄片免| 免费观看人在逋| 久久国产精品人妻蜜桃| 亚洲乱码一区二区免费版| 俺也久久电影网| 国产精品亚洲美女久久久| 国产精品免费一区二区三区在线| 日本撒尿小便嘘嘘汇集6| 精品不卡国产一区二区三区| 少妇裸体淫交视频免费看高清| 好看av亚洲va欧美ⅴa在| 99在线人妻在线中文字幕| 丰满人妻熟妇乱又伦精品不卡| 搡老岳熟女国产| 亚洲成a人片在线一区二区| 国内少妇人妻偷人精品xxx网站| 国产亚洲精品久久久com| 亚洲色图av天堂| 制服丝袜大香蕉在线| 久久精品国产综合久久久| 女生性感内裤真人,穿戴方法视频| 亚洲在线自拍视频| 男女做爰动态图高潮gif福利片| 亚洲第一欧美日韩一区二区三区| 成人国产一区最新在线观看| 亚洲精品影视一区二区三区av| 天天一区二区日本电影三级| 琪琪午夜伦伦电影理论片6080| 国产色爽女视频免费观看| 久久久成人免费电影| 老司机福利观看| 亚洲国产精品sss在线观看| 最新在线观看一区二区三区| 国产精品影院久久| 亚洲成av人片免费观看| 白带黄色成豆腐渣| 久久6这里有精品| 亚洲va日本ⅴa欧美va伊人久久| 超碰av人人做人人爽久久 | 精品人妻偷拍中文字幕| 亚洲欧美日韩高清专用| 成年女人看的毛片在线观看| a在线观看视频网站| 日本成人三级电影网站| 欧美成人a在线观看| 免费无遮挡裸体视频| 色播亚洲综合网| 精品福利观看| 色播亚洲综合网| 激情在线观看视频在线高清| 亚洲av熟女| 欧美精品啪啪一区二区三区| 精品免费久久久久久久清纯| 特级一级黄色大片| 天美传媒精品一区二区| 欧美绝顶高潮抽搐喷水| 一夜夜www| 欧美日韩瑟瑟在线播放| av视频在线观看入口| 白带黄色成豆腐渣| 天天躁日日操中文字幕| av天堂中文字幕网| 日韩欧美精品免费久久 | 日韩中文字幕欧美一区二区| 嫁个100分男人电影在线观看| 亚洲国产高清在线一区二区三| 午夜福利视频1000在线观看| 国产极品精品免费视频能看的| 国产成人影院久久av| 成人特级黄色片久久久久久久| 美女cb高潮喷水在线观看| 国产99白浆流出| 亚洲国产色片| 色哟哟哟哟哟哟| 国产精品久久视频播放| 国产精品99久久久久久久久| 九九热线精品视视频播放| 免费观看精品视频网站| 久久久久久久久久黄片| 精品人妻一区二区三区麻豆 | 成人18禁在线播放| eeuss影院久久| 亚洲精品456在线播放app | 一个人免费在线观看的高清视频| 夜夜躁狠狠躁天天躁| 亚洲 欧美 日韩 在线 免费| 夜夜爽天天搞| av福利片在线观看| 十八禁人妻一区二区| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲精品久久久com|