CN104870780B - 爆震驻波发动机 - Google Patents

爆震驻波发动机 Download PDF

Info

Publication number
CN104870780B
CN104870780B CN201380063447.7A CN201380063447A CN104870780B CN 104870780 B CN104870780 B CN 104870780B CN 201380063447 A CN201380063447 A CN 201380063447A CN 104870780 B CN104870780 B CN 104870780B
Authority
CN
China
Prior art keywords
pinking
entrance
detonation engine
detonation
fuel
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
Application number
CN201380063447.7A
Other languages
English (en)
Other versions
CN104870780A (zh
Inventor
A.卡西莫夫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
King Abdullah University of Science and Technology KAUST
Original Assignee
King Abdullah University of Science and Technology KAUST
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by King Abdullah University of Science and Technology KAUST filed Critical King Abdullah University of Science and Technology KAUST
Publication of CN104870780A publication Critical patent/CN104870780A/zh
Application granted granted Critical
Publication of CN104870780B publication Critical patent/CN104870780B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23R—GENERATING COMBUSTION PRODUCTS OF HIGH PRESSURE OR HIGH VELOCITY, e.g. GAS-TURBINE COMBUSTION CHAMBERS
    • F23R7/00—Intermittent or explosive combustion chambers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C5/00—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion
    • F02C5/02—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion characterised by the arrangement of the combustion chamber in the chamber in the plant
    • F02C5/04—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion characterised by the arrangement of the combustion chamber in the chamber in the plant the combustion chambers being formed at least partly in the turbine rotor
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C5/00—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion
    • F02C5/02—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion characterised by the arrangement of the combustion chamber in the chamber in the plant
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C7/00—Features, components parts, details or accessories, not provided for in, or of interest apart form groups F02C1/00 - F02C6/00; Air intakes for jet-propulsion plants
    • F02C7/22—Fuel supply systems
    • F02C7/222—Fuel flow conduits, e.g. manifolds
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02K—JET-PROPULSION PLANTS
    • F02K7/00—Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof
    • F02K7/02—Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof the jet being intermittent, i.e. pulse-jet
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02K—JET-PROPULSION PLANTS
    • F02K9/00—Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof
    • F02K9/42—Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof using liquid or gaseous propellants
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02C—GAS-TURBINE PLANTS; AIR INTAKES FOR JET-PROPULSION PLANTS; CONTROLLING FUEL SUPPLY IN AIR-BREATHING JET-PROPULSION PLANTS
    • F02C5/00—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion
    • F02C5/10—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion the working fluid forming a resonating or oscillating gas column, i.e. the combustion chambers having no positively actuated valves, e.g. using Helmholtz effect
    • F02C5/11—Gas-turbine plants characterised by the working fluid being generated by intermittent combustion the working fluid forming a resonating or oscillating gas column, i.e. the combustion chambers having no positively actuated valves, e.g. using Helmholtz effect using valveless combustion chambers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02K—JET-PROPULSION PLANTS
    • F02K7/00—Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof
    • F02K7/02—Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof the jet being intermittent, i.e. pulse-jet
    • F02K7/06—Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof the jet being intermittent, i.e. pulse-jet with combustion chambers having valves
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2220/00—Application
    • F05D2220/30—Application in turbines
    • F05D2220/32—Application in turbines in gas turbines
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2240/00—Components
    • F05D2240/35—Combustors or associated equipment
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23C—METHODS OR APPARATUS FOR COMBUSTION USING FLUID FUEL OR SOLID FUEL SUSPENDED IN  A CARRIER GAS OR AIR 
    • F23C15/00—Apparatus in which combustion takes place in pulses influenced by acoustic resonance in a gas mass
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23R—GENERATING COMBUSTION PRODUCTS OF HIGH PRESSURE OR HIGH VELOCITY, e.g. GAS-TURBINE COMBUSTION CHAMBERS
    • F23R2900/00—Special features of, or arrangements for continuous combustion chambers; Combustion processes therefor
    • F23R2900/03281—Intermittent fuel injection or supply with plunger pump or other means therefor

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluidized-Bed Combustion And Resonant Combustion (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Jet Pumps And Other Pumps (AREA)

Abstract

一种爆震发动机(100)，能够在圆柱形爆震区域(205)内起爆燃料与氧化剂的混合物以做功。爆震发动机可具有第一和第二入口(135、140)，其具有从储箱流体地连接到爆震发动机的端部。所述第一和第二入口可沿共同轴线(145)排列。所述入口可以连接到喷嘴(105、110)，而分离器(115)可以放置在喷嘴之间并沿着共同轴线。

Description

爆震驻波发动机
优先权要求
本申请要求2012年10月12日提交的美国临时申请编号61/712,972的优先权，其全部内容通过引用并入。
技术领域
本发明涉及一种爆震驻波发动机(standing detonation wave engine)。
背景技术
气相爆震波可以在各种几何结构中传播并被其影响。在圆柱形刚性管中，所述传播在大直径管中可以是多元胞(multi-cellular)爆震的形式，在中等直径管中可以是以旋转爆震的形式，和/或在非常小直径的管中可以是驰振爆震(galloping detonation)。在各种宽高比的矩形横截面通道中，类似的元胞模式或驰振模式也可以达到。在平行板之间，当板之间的间隙比所述板的横向跨度小得多时，因为在垂直于板的方向上的横波被抑制，可以实现二维元胞爆震。其他的构造是那些在超音速流中稳定的爆震，其可以与爆震推进和爆震发动机的问题有关。
爆震燃烧是燃烧燃料与空气的混合物以释放化学能的有效方法。通过用输出功除热输入而计算出的爆震燃烧的理论效率大约是49％。作为比较，诸如恒定体积燃烧或恒压燃烧的更传统工艺分别具有47％和27％的理论效率。爆震燃烧的效率提高归因于它独特的热释放过程，其中燃料-空气混合物的燃烧发生的比常规燃烧快数万倍，这依赖于传播的火焰峰(flame front)。虽然更有效率，爆震燃烧也可以是更加难以控制的。例如，困难在启动和维持爆震燃烧时出现。然而，在发动机控制技术中的最新进展允许克服这些困难。
发明内容
所提供的是爆震发动机和方法。爆震发动机可产生稳定的爆震驻波。由稳定的爆震波产生的加热燃烧产物的流动可以被转换成机械能。例如，所述爆震发动机可具有爆震波在其中被稳定的燃烧腔室。该腔室可部分地被例如涡轮叶片围绕，所述涡轮叶片可提供障碍物以限制爆震波的自由膨胀，从而有助于形成稳定的驻波。
在一个方面，爆震发动机可包括具有流体地连接到第一储箱的第一端部和流体地连接到爆震发动机的第二端部的第一入口、具有流体地连接到第二储箱的第一端部和与第一入口相对地流体连接到爆震发动机的第二端部的第二入口。所述第一和第二入口可与共同轴线对齐。该爆震发动机还可包括连接到所述第一入口的第一喷嘴、连接到所述第二入口的第二喷嘴，以及布置在所述第一入口的第二端部与所述第二入口的第二端部之间分离器。第一入口、第二入口和分离器可以沿共同轴线定位。
在一些实施例中，爆震发动机可包括构造为在所述爆震发动机中稳定爆震的障碍物。爆震发动机可包括可旋转地附接到所述爆震发动机的一个或多个涡轮。涡轮还可以包括构造为在所述爆震发动机中稳定爆震的叶片。所述第一储箱可被构造为接收氧化剂。所述第二储箱可被构造为接收燃料。第一喷嘴可以靠近第一入口的第二端部。第二喷嘴可以靠近第二入口的第二端部。第一喷嘴可布置在第一入口内。第二喷嘴可布置在第二入口内。
在另一方面，一种用于引爆燃料与氧化剂的混合物的方法可以包括通过第一入口引入氧化剂到爆震发动机，通过第二入口引入燃料到爆震发动机，其中第一和第二入口可与共同轴线对准，然后在距共同轴线一定距离处的圆柱爆震区域内引爆燃料与氧化剂的混合物。一种方法可以包括加速氧化剂通过第一喷嘴，并加速燃料通过第二喷嘴。一种方法可以包括靠着分离器的第一侧并且从共同轴线径向向外地引导所述氧化剂，以及靠着分离器的第二侧并且从共同轴线径向向外地引导所述燃料，所述第二侧与第一侧相对。
在一些实施例中，第一入口具有流体地连接到第一储箱的第一端部和流体地连接到爆震发动机的第二端部。第二入口具有流体地连接到第二储箱的第一端部和流体连接到爆震发动机的第二端部。所述第一和第二入口可与共同轴线对齐。加速氧化剂通过第一喷嘴可以产生超音速流，并且加速燃料通过第二喷嘴可产生超音速流。一种方法可以包括膨胀由引爆混合物产生的燃烧产物通过涡轮以做功(produce work)。一种方法可包括提供障碍物以抑制由引爆混合物产生的燃烧产物的膨胀。一种方法可以包括在爆震发动机中稳定混合物的爆震。障碍物可以构造在爆震发动机中以实现爆震的稳定化。障 碍物可以是涡轮叶片，其可以，例如，将燃烧产物的动能转化为机械能。
其他方面、实施例和特征将从下面的描述、附图和权利要求书中显而易见。
附图说明
图1是爆震发动机的横截面视图。
图2是示出示例性圆柱形爆震区域的爆震发动机的横截面视图。
图3是在爆震中以压力表示的模拟驻波。
图4是在爆震中以温度表示的模拟驻波。
图5描绘了在径向膨胀流中的驻定爆震(standing detonation)的几何形状。
图6描绘了绝热膨胀和爆震冲击波的位置-速度曲线。
图7a)至i)描绘了作为以马赫数给出的入流的函数的爆震半径。
图8描述了示例性稳定爆震剖面(detonation profiles)。
图9描绘了作为初始流的函数的示例性爆震半径。
图10描述了作为E和Q的函数的最小爆震半径。
图11描绘了以方波状稳定解开始的爆震崩溃(collapse)。
图12描绘了以方波状稳定解开始的爆震膨胀。
图13描绘了膨胀爆震的稳定化。
图14描绘了在障碍物的存在时的爆震波启动(initiation of a detonationwave)。
图15描绘了反应区的详细结构。
具体实施方式
在实验性发动机中实现爆震驻波可以是挑战性的，但使用计算模型可以显著地提高成功的可能性。例如，通过迭代方法，计算模型可以使研究者能够无需建造一系列原型而确定有效几何形状和边界条件。该计算模型应当考虑守恒定律、化学反应机制以及稳定性考虑。爆震波可以被建模为跟随着反应区的冲击波。为了简化计算并减少计算时间，爆震波可以被建模为气体分子量没有变化的理想流体。能量增加可以被假定为发生在单步化学反应中。为了进一步简化计算，热传导和粘度的作用可以被忽略。所得系统的二维反 应欧拉方程能够例如在并行计算机集群上有效地解决，以便测试爆震驻波的稳定存在性。
在爆震发动机中，反应混合物的超音速流可以离开德拉瓦尔喷嘴(de Lavalnozzle)，并能形成，例如，马赫盘。在马赫盘中的气体压缩可以启动下游化学反应，使得反应区形成在离所述盘一定距离处。这种流构造的一个早期实验研究被在[10]中报道，其中，作者能够实现驻定冲击-反应-区域复合体。即使这种构造类似于爆震，它也许更适合被称为冲击诱导的燃烧，而不是爆震，因为马赫盘能够在这样的流中稳定地存在，而与化学反应的存在无关。尽我们所知，反应区在马赫盘的存在和性质中起作用的程度似乎没有被很详细地探索。
在推动燃料和氧化剂通过例如会聚-扩散喷嘴进入测试通道的高压供应罐的帮助下，燃料-空气混合物的超音速流可以产生。对与所述流一起传播或逆所述流传播的爆震进行了研究，并对通道内的流边界层在爆震速度上的影响进行了研究。在[21]中，通过改变通道几何形状的稳定化使用氢气-空气的多步动力学模型而数值地考虑。他们的结论是，通道横截面的变化可以用来稳定在通道中的爆震波。
爆震也可以被钝体在超音速流中稳定。在这样的构造中，爆震可以被稳定在反应气体的超音速流中的楔形体上。化学反应可以在楔形体前端下游的一定距离处开始，并可以影响附着到楔形体的冲击波的结构。类似于在德拉瓦尔喷嘴中的马赫盘下游的反应的情况，在这种构造中的冲击波即使没有化学反应也可以存在，所述化学反应即冲击诱导的燃烧。
通过在轴向上固定其位置，但允许它在方位角方向上旋转，爆震可以被稳定在超音速流中。这样的构造可以实现连续旋转爆震。超音速空气流可以进入两个同轴圆柱之间的薄缝隙并与燃料混合。燃料可以从例如内圆柱注入。该混合物可以在旋转爆震波中连续地燃烧，所述旋转爆震波在圆周方向上传播。这种构造可以在爆震发动机中使用，例如，作为脉冲爆震发动机的替代。
如下面进一步讨论的，对根据单步阿兰尼乌斯(Arrhenius)反应模型的理想气体反应使用可压缩欧拉方程，可以发现在从中心源径向发出的超音速流中驻定的稳态爆震的存在。在稳定爆震解的存在和结构中，入流的情况，诸如混合物的停滞焓和入流的马赫数，可以在爆震波的稳定性中发挥作用。 取决于所述参数，人们或者可以没有稳态问题的解，或者可以有它的一个或两个解。在两个稳定解共存的情况下，所述解中的一个可以对应于相对较小的爆震半径，而另一个对应于较大的半径。这两种类型的解可以具有明显不同的结构。小半径的解可以具有方波状结构，所述结构具有约是稳定平面ZND解长度的100倍的大感应区(induction zone)。与此相反，大半径的解没有明显的感应区。
此外，例如通过对二维反应欧拉方程进行数值积分，稳态解的稳定性可以被确定并且爆震的动力学可以被计算。作为初始条件，无论是方波状还是常规的稳定解都在下面进行了考虑。可以得到在几十个稳定半反应时间的时间尺度上的快速崩解，以及在几百个相同时间单位的尺度上的缓慢膨胀。大半径的解可以导致非常缓慢的膨胀爆震，其中，元胞结构在任何显著膨胀发生之前迅速地形成。这些膨胀爆震可以通过将小的刚性障碍物放置在所述稳态爆震半径的下游而被稳定，使得稳定的元胞爆震建立在距离所述源中心的一定距离处。
通过将障碍物最初放置在气体的非反应绝热流中而通过障碍物的爆震启动(Detonation initiation)也可以被数值地确定。例如，围绕障碍物形成的弓形冲击波可以快速地转化为开始膨胀的爆震。用相同的障碍物，在启动情况下可以获得相同的最终爆震结构。
一种实施例可以包括可燃混合物的超音速流。超音速流可以从圆形源径向地流动，并在下游经历爆震燃烧。所述流可以在两个平行板之间被引导，外部的混合物源从中心提供入流。在源出口处的高速流条件可通过经由例如喷嘴的快速扩张产生。由于冲击波的存在依赖于化学反应的存在，在本实施例中的爆震可以是自持波。所述流的单纯绝热膨胀会在没有化学反应存在的情况下发生。
控制反应欧拉方程可以允许在距离源一定有限距离处驻定的有自持爆震的稳态解。这样的解的性质可以通过检查该问题的各种参数而探究，所述参数诸如混合物性质和入流条件。使用二维模拟，可以对爆震的非线性动力学进行研究，所述动力学例如它的稳定性。爆震通常是不稳定的，并且所述不稳定不仅以形成元胞的形式表现其本身，也以爆震峰整体径向收缩和膨胀的形式表现其本身。膨胀可以通过例如在流的下游放置几个刚性障碍物而防止。
在没有化学反应存在的情况下，从给定半径的圆形源发出的理想可燃气体的二维径向对称超音速流是绝热的，并且由于它也是超音速的，在膨胀过程中，流速和马赫数随着距离增加，而压力、温度和密度都随着距离下降。如将在下面示出的，这些特征可以从气体动力学方程建立。利用例如在ZND理论的框架内的分析，稳定的径向对称爆震波可以在何种条件下存在于绝热径向膨胀流中可被确定。利用例如有单步阿兰尼乌斯动力学的二维反应欧拉方程的数值解，可以显示出这样的稳定结构对于二维扰动的稳定性。
在径向膨胀流中的驻定爆震的示意性几何结构示于图5中。中心源具有半径r0，反应气体在由压力p0、密度ρ0和流速u0给定的初始流条件下从所述源发出。由于所述流在超音速膨胀期间加速，在一些距离rs处，所述流条件可以变成使得爆震驻波结构可以建立在rs的下游。在爆震冲击波前的状态为ρ1、p1、u1，而紧随所述冲击波之后的状态为ρ2、p2、u2。对于驻定爆震，它的速度D应该与u1相同。在冲击波之后，有一个在r＝r*处的声速轨迹，其中的流速等于局部声速u*＝c*。
如以下所讨论的，这样的驻定爆震结构可在广泛的条件下存在。此外，在相同入流条件下的多个解的共存可以成为可能。
如果从所述源发出的气体的温度足够低，并且由于它在膨胀过程中下降，来自所述源的流可视为绝热的并且反应可以忽略不计。爆震半径rs先验地是未知的，但可以例如通过将上游状态与蓝金-许贡纽(Rankine–Hugoniot)条件和爆震冲击波的下游流条件匹配而确定。这样的爆震结构的一个重要成分是在冲击之后的声速点的存在。因为紧随所述冲击之后的流速是亚音速并且距冲击足够远的产物流动是超音速的，由此得出的结论是，在距源一定距离r*>rs处，流速等于局部声速。
气态爆炸速度通常是每秒几公里的量级，并且为了将这样爆震保持在固定距离处，所述流的初始能量可做得足够高。建立驻定爆震必需的初始流总能量可以估计。假设没有摩擦损失，所述流的总焓，
是一个守恒量。这里，T是所述流的温度，u是流速，W是该混合物的摩尔质量，R是通用气体常数。在绝热膨胀过程中，流的温度可以降低，因此流 速可以增加。即使流的所有势能转换成动能，所述流不能比更快。在另一方面，爆震速度具有可以估计的下限。假设会聚爆震的半径比反应区的尺寸大得多，爆震曲率效应可以被忽略不计，而爆震速度可以近似为 其中T是在爆震冲击前面的环境温度。由此，爆震速度总是大于为了将爆震保持在固定距离处，环境流应该加速到爆震速度，因此，对于初始流能量的条件可以写为H>Q(γ2-1)。
反应欧拉方程
二维可压缩反应理想气体可以假设为通过反应欧拉方程系统描述，其由质量、动量和能量守恒方程以及化学反应方程组成。
其中，ω(p,ρ,λ)是反应速率，其可以假设为通过反应物→产物的简化模型描述，速率通过下式给出
其中，E是活化能，ρ是密度，p是压强，k是反应速率常数，以及λ是反应进展变量。反应物的质量分率是1-λ，而产物的质量分率是λ，其中λ＝0对应于新鲜混合物，而λ＝0对应于完全燃烧气体。状态方程通过下式给出
其中，γ是比热的恒定比率。(4)中的总能量然后被限定为
冲击条件
蓝金-许贡纽跃变条件(jump condition)是
其中，D是冲击速度的垂直分量，而un是流速的垂直分量。像往常一样，冲击自身是非反应性的，因此(11)平凡地满足。对于稳态爆震，爆震速度等于环境流速，即D＝u1，因此，爆震在实验室参考系中是静止的。对于圆形稳态解，所述流是垂直于冲击面的，因此对于理想气体，所述蓝金-许贡纽条件可以写为：
其中，ρ1、u1、p1、ρ2、u2、p2分别是在冲击前和冲击后的密度、速度和压强。从这些方程，冲击后的状态可以明确地根据冲击前状态写出。
如在下面所讨论的，稳态系统可以化简为对u和λ的两个ODEs(ordinarydifferential equations，常微分方程)，和对于总能量和质量的两个守恒定律。通过假定反应速率在上游流为零，可以使用仅针对流速的蓝金-许贡纽条件，它有很简单的形式：
其中，是在源处的总能量，其在稳态中沿着所述流是守恒的。
无量纲方程和参数的选择
在爆震传播到参数恒定的静态中的问题中，相对于恒定状态或相对于冲击后状态缩放参数是自然的。爆震能够在非均匀介质中驻定或传播，由此，尺度的最佳选择是不会立即显现的。可以选择参考压力Pa、参考密度ρa、温度Ta＝Wpa/ρa以及速度这些可以视为对应于1个大气压和300K。通常地，独立于这样的参考状态的具体选择，并用如下所述的长度和时间尺度的选择，控制方程可以保持其形式。
剩余的尺度是用于传播到上述参考状态中的平面爆震的长度尺度，为此可以选择半反应区长度l1/2，以及时间尺度t1/2＝l1/2/ua。对于一组给定的反应参数Q、E和γ，设定这些尺度等于通过积分固定速率常数，
其中，
并且，
由此，在下面的计算中，尺度和时间尺度可以通过k确定，其根据基于p＝1、ρ＝1和T＝1的无量纲上游状态的Q、E和γ的值变化。
稳态径向对称解
在稳定径向对称的情况下，运动方程变为：
这些方程可以简化到对于u和λ的ODEs，
而质量和焓是守恒量：
rρu＝M＝constant
(22)
这里是局部声速，r是径向坐标以及M是质量通量。在下面的计算中，方程20和21可重写为三个方程的自治系统，其中未知数通过τ参数化，
这个系统以及对于质量和焓(方程22和23)的两个守恒定律，连同在源和远场处的边界条件和冲击条件(方程15)可以完全地确定驻定爆震结构。在源与爆震冲击之间的环境状态是处于绝热膨胀中。稳态方程可以简化为一个代数方程，例如，对于流速：
其中，r0是源的半径，并且u0是在源处的流速。
方程(20)的形式，其中有在所述流中存在声速点的可能性，在那里u＝c，并且通过在同一点设定方程(20)的分子为零而正规化的可能性，可以给出解的存在条件。原则上，求解过程需要人们对于给定入流条件发现整个结构。然而，无论是冲击位置，还是声波点的位置都不能直接地算出。人们需要一个迭代过程以确定后者。
由于声速点是鞍点，数值上更强劲(numerically more robust)的可以是，在声速点附近找到作为泰勒级数展开的解，然后从它离开并继续通过常规数值方法积分。声速点位置可以明确地是未知的。在源处固定边界条件，方程系统(24)在声速点具有一个猜测参数。它可以是声速点的半径r*，爆震半 径rs，或反应进展变量在声速点的值λ*。为了满足在源处的边界条件，人们可以扫描这些参数中的一个，例如λ*。
为了说明所述算法，假设源具有半径r0，所述流在源处具有压强p0、密度ρ0，并且超音速流速u0>c0。一旦在源处的流是已知的，人们就可以计算出H和M，其在整个域中是守恒的。因为H是固定的，通过设定(24)(a)中的分子至零，并使用c*＝u*，人们可以写出用于声速点的半径r*关于λ*的表达式：
其中，反应速率和声速通过下列表达式仅取决于λ*：
λ*的值在这个算法中可以是猜测参数。通过例如在0与1之间扫描它，在源处的初始条件可以得到满足。如果λ*在这一过程中被发现，由于H是已知的，人们可以将方程系统(24)积分到混合物变得新鲜的点λ(τs)＝0。这一点可以限定反应区的开始，这是后冲击状态。因此，跃变的半径由rs＝r(τs)给出。接着，施加跃变条件(15)，人们可以很容易地通过计算方程式(25)获得冲击之前的状态和所有的绝热剖面。
稳态爆震解可以只对在源处初始参数的一定范围存在。从另一个角度来考虑扫描过程，H可被假定是固定的。然后，对在声速点的反应进展变量的值的范围0<λ*<1解方程(24)可以给出后冲击流速对爆震半径的依赖关系。这个函数是单调递减和凸形的。使用方程(25)和(15)可对绝热膨胀计算出相同的依赖性。后一曲线也是凸形和单调递减的，但是它包含一个附加参数，在源处的初始流速调整u0，绝热曲线可相对于所述爆震曲线移动。在初始流速的一些值下，所述曲线相交，这意味着解的存在。这些曲线可以有多达两个的交点，这对于在源处的特定流给出了多个稳态解的存在性。不存在交点意味着所述问题不存在稳定解。
在图6中显示了三个条件。绝热膨胀和爆震冲击波的位置-速度曲线被 示出，其中曲线和有两个交点(a)，一个交点(b)，以及没有交点(c)。混合物的参数选择为：γ＝1.2；E＝40；Q＝30；H＝1.3Hmin；并且r0＝50。在图6中，(a)M0＝4.40，(b)M0＝4.0，以及(c)M0＝4.67。
上述讨论的一个重要结果是，对于特定混合物的爆震半径可只取决于两个参数：所述流的总焓的值以及初始流速的值。在这些参数的某些值下，可能不存在解，而在其他值下，一个或两个解是可能的。
稳态解的存在性和结构
可以是重要的是，确定何时存在解，并且如果存在，爆震位于哪里。在源处的混合物特性和入流条件在稳态解的结构上可以发挥作用。
图7描绘了爆震半径作为以马赫数给出的入流的函数。对于各种E、Q和γ，爆震半径被给出为所述源流的马赫数M0的函数。源的半径为r0＝50，而停滞焓H＝1.3Hmin＝1.3Q(γ2–1)。图7示出了在源的固定半径下和停滞焓的固定值下，入流马赫数可以如何影响爆震半径。附图还显示了活化能E、放热Q和比热比率γ发挥的作用。相对小的半径和相对大的半径的存在在这些图中是典型的。在γ＝1.2时，渐增的Q需要更大的M0值，以实现稳定解。与此同时，当Q从10改变到30时，上半径从超过9000到约1000以一个数量级减小。这可以表明入流马赫数是足够大的；半径是反应区尺寸的约1000倍大的爆震是可行的。渐增的γ具有不仅减少爆震的上半径到反应区的同样1000倍，而且还显著地降低入流马赫数的有趣效果，所述马赫数对于γ＝1.2从3-5下降到1-3。活化能的效果可以是非单调的，即随着活化能的增加，爆震半径可以先增加后降低。一个有趣的例子由在E＝30下的图7(d)表现。即使在源马赫数为1时，也存在两个解，下半径约是rs＝100而上半径约是3000。
如果混合物焓H是小的，即接近于最小值Hmin＝Q(γ2–1)，则能量不足以将所述流加速到高速直至无穷。那么，rs(M0)曲线的顶分支不存在，但底分支仍然存在。如果在另一方面，H是非常高的，那么顶分支是可行的，但底分支消失，原因是对应的半径变得小于源半径。这是在图7的一些图中没有底分支的原因。
探索对应于图7中所示的上解和下解的稳态解的结构是感兴趣的。在图8中，所绘制的是对应于下解(左栏)以及对应于在相似流条件下的上解的p、u、T和1-M的曲线。对于两者，所述解是方波形状的。
在存在两个分支时，左栏对应于在rs–M曲线的下分支上的方波爆震。左栏的参数是：γ＝1.2、Q＝10、E＝30、r0＝50、ρ0＝1、p0＝1.40、u0＝1.30、M0＝1.0、H＝2.1Hmin。对于右栏，仅存在上分支。右栏的参数是：γ＝1.4、Q＝10、E＝30、r0＝50、ρ0＝1、p0＝2.70、u0＝3.90、M0＝2.0、H＝1.75Hmin。
左边所示的解的一个显著特点是具有相当长的反应区的方波状结构。在冲击之后直到薄的能量释放区域，可看到压力、温度、速度和马赫数的曲线表现出几乎恒定的状态。有延伸约30个长度单位的感应区。随后，所有的能量可以在几个长度单位的距离上释放。与此相反，上解的结构没有显示任何感应区，并且反应区是尖锐的，仅跨越几个长度单位。混合物的性质对于这两种情况是相同的，然而，入流条件是不同地选择的，以便将爆震置于距原点大约相同的距离处。因此，对于相同的混合物，入流条件可被修改，使得两个非常不同的驻定爆震可以在同一半径处存在。它们的稳定性是这样的爆震是否存在的决定性因素。
在前面的计算中，入流焓是固定的，而对混合物参数和入流马赫数的作用进行了探讨。对于给定混合物，入流焓和入流马赫数在爆震半径上的影响在下面进行讨论。结果显示在图9中。
对于停滞焓H的不同值，图9示出了作为初始流马赫数的函数的爆震半径。混合物的参数是：γ＝1.3、Q＝10、E＝30。源的半径是r0＝50。通过增加混合物的焓，上解可以接近爆震半径的下值。下解的半径也可以随着H的增加而降低。这一结果与流的焓越高，就能够越容易地将其加速到在所述流中建立稳定爆震所必需的速度的预期是一致的。
在给定的M0处存在爆震解的Q–E参数空间的区域可以被确定。在图10中显示的是对于两个不同的流入马赫数M0＝1和M0＝2在γ＝1.2和γ＝1.4处作为Q和E的函数的爆震最小半径。顶部的数字对应于γ＝1.2而底部的数字对应于γ＝1.4。左边的数字对应于M0＝1而右边的数字对应于M0＝2。源的半径为r0＝50并且H＝1.3Hmin。在这些图中的一个有趣特征是，对于较大的γ，解的存在范围广得多，并且最小半径比γ＝1.2的情况要小得多。
会聚爆震的半径可取决于在源处的流并且取决于混合物的参数。为了找到依赖关系，对于这些参数的范围，可以解完整的方程系统(24)。在解这个问题之前，可以做出一些分析性的估计。例如，对于驻定爆震，在某些点 处的流速应该等于爆震速度，即uCJ＝D。当爆震半径较大时，爆震速度可通过例如查普曼-约戈特(Chapman-Jouguet)公式估计，
其中，T是就在爆震冲击前的温度。只要焓是固定的，绝热膨胀的环境温度只是流速的函数，
并且因此，人们可以对在爆震冲击位置处的绝热流速写出以下方程：
此方程可关于uCJ求解，并且它的解与方程(25)一起给出了爆震半径对初始流速的依赖关系，即rs(u0)。这个函数单调衰减：在源处的流越快，它越早到达查普曼-约戈特速度DCJ。这个估计可以给出rs–M0曲线的顶部分支的基本准确的描述。
二维模拟
虽然对于所述问题参数的广阔范围存在稳态解，了解它们的稳定性可以是重要的。气态爆震对于多维扰动往往是不稳定的。然而，本构造的两个元件可以起到本文阐明的作用，即爆震冲击的非均匀上游流和爆震冲击曲率。
基于以下简单的论证，人们可得出的结论是，事实上，即使相对于纵向扰动，在我们的构造中的爆震应始终是不稳定的。考虑在离中心一定距离处驻定的稳定径向对称爆震，并且想象它朝向所述源向内扰动小的距离。那么，由于在扰动后上游流比扰动之前更慢，并且由于爆震趋于相对于上游流在恒定速度下传播，扰动冲击将继续向内移动，这意味着不稳定。在爆震冲击向外扰动的相反情况下，扰动冲击的上游流比扰动之前更快。因此，出于爆震趋于以相对于上游状态的恒定速度传播的同样原因，扰动冲击将继续扩大，这又意味着不稳定。
但是，这些简单的论证忽视了两个可以在爆震动力学起作用的重要效应：曲率在爆震速度上的影响和导致爆震元胞(detonation cells)形成的横向不稳定性。随着曲率增加，爆震速度可以降低。由此，当爆震冲击被扰动 时，例如向内，它的稳态速度会下降。这样的下降对于补偿上游流的下降速度可以是足够的。由此，一旦爆震向内扰动，新的稳定状态是可能的。如果曲率效应超过了上游流动速度下降的影响而占支配地位，扰动实际上可以减少，导致爆震的稳定性。
在爆震半径较大时，曲率效应可以是较弱的。重要的多的是，两维不稳定可以到来(can set in)，引起高度非平凡的多维动力学，其中元胞结构开始发挥主导作用。正如下面所讨论，爆震的二维演化可以开始于对应方波状和常规爆震结构的稳态解。所述径向对称解在所有情况下都是不稳定的，但是，对于两种稳定解，不稳定性的本质是不同的。崩溃解和膨胀解两者都可以找到，它们之间的一个重要区别由不稳定性的时间尺度给出，其在崩溃情况下远小于在膨胀情况下的。
膨胀爆震可以通过在距中心一定距离在声速点之后放置几个障碍物来稳定。另外，障碍物稳定爆震可以通过超音速流被障碍物阻碍而启动。障碍物可以引起弓形冲击，其中，爆震被启动，并且可以通过连接来自障碍物的独立爆震峰形成围绕中心源的单个峰而建立其自身。
对于二维模拟，可以利用解算器，诸如由泰勒、卡西莫夫和斯图瓦特(Taylor,Kasimov,and Stewart)开发的。解算器可以使用有限体积法，其中，流量可以例如通过五阶WENO算法计算，而时间积分可以例如通过三阶TVD Runge-Kutta方法完成。空间域可以离散为均匀笛卡尔网格，分辨率是每稳定爆震半反应长度至少20个网格点。所述CFL数可以是0.5。代码可以例如通过使用虚源反射法(ghost-cell method)设计为分布式并行体系。在源处的流入边界条件可被设定，并且在域的端部的流出条件可通过变量外推设定。障碍物可以假定为绝对刚体，并且它们的边界可以使用例如浸入边界方法处理。
稳态圆形爆震的不稳定性
为了理解稳态解的不稳定性，可以分析如图8和图11所示的两种情况。图11示出了作为图8中所示的方波状稳定解开始的爆震的崩溃。从(a)至(c)的快照的时间分别是：t＝1、10和40。域的尺寸为600×600，而网格点的数量为1280×1280，其对应于每半反应区64个点。爆震的初始半径大约是150。经过t＝40较短时间，半径已减少到约100。对附图的仔细观察显示出，事实上，所述波已经经历了二维不稳定，并且爆震元胞出现。然而， 所述元胞是较弱的，不会明显改变冲击波的圆形形状。冲击压强在崩溃过程中从t＝0时的约1.1增加至t＝40时的2.3。动力学保持了基本的径向对称。
与上面考虑的崩溃情况相反，图12描绘了膨胀爆震。作为非方波状稳定解(如图8所示)开始的爆震的膨胀在图12中示出。从(a)至(c)的快照的时间分别是：t＝0、150和400。域的尺寸为600×600，而网格点的数量为5120z 5120，其对应于每半反应区20个点。初始条件是图8(右栏)的。爆震的初始半径与崩溃情况几乎是一样的，但爆震反应区的结构是非常不同的，压强在先导冲击之后急剧下降并且没有可见的感应区。在这种特殊情况下的二维不稳定是相当强的，并迅速导致强烈的元胞爆震的形成。这种情况与先前的崩溃情况的最重要区别在于，膨胀是显著更慢的。爆震扩大到两倍的初始半径需要大约400个时间单位。这种缓慢的膨胀表明，通过在冲击下游的流中放置减慢反应产物的膨胀流的障碍物可能能够防止所述膨胀。
图13示出了图12的膨胀爆震通过三个半径为10的障碍物的稳定化。从(a)至(c)的快照的时间分别是：t＝10、50和700。域的尺寸为600×600，而网格点的数量为2500×2500，其对应于每半反应区10个点。
爆炸通过障碍物的稳定化
为了查看前一节中的膨胀爆震是否确实可以被稳定化，几个障碍物可以被放置在流的稳态声速轨迹的正下游。会出现许多可能性，例如，取决于障碍物的数量、尺寸和形状。然而，重要的是，不多的几个障碍物足以防止膨胀。在障碍物前形成的弓形冲击能减缓产物流，使得爆震冲击在源与障碍物之间的区域中能够保持稳定。所得到的爆震波的精确位置和形状可以取决于障碍物的选择和混合物的细节以及源条件。
作为一个例子，在图13(a)中，小的弓形冲击的发展可围绕障碍物形成，而在图13(c)中，大的三角状弓形冲击的结构可在反应区下游驻定，提供对驻定爆震的稳定支持。在这种特殊情况下，在相同的半径上使用等间隔的多于三个的障碍物会导致所述波的崩溃。所述解是崩溃还是稳定可以敏感地取决于在现象中发挥作用的众多参数。
爆震的启动
计算完全可以将稳态解作为二维模拟的初始条件。这样的计算可以提供对给定稳态解的稳定性特征的理解。不稳定的稳态解可以是难以实现的，因为，例如，这样的爆震是通过源启动的，在不稳定的情况下，这可以无需进 入稳态解而导致脉冲或元胞爆震。由此，在径向出流中的驻定爆震的启动可以是重要的。可以实现这样做的不同手段，例如，通过从来自源的非反应性超音速流在这样的流遇到刚性障碍物后的爆震启动。
上面所讨论的同样障碍物可以放置在最初是非反应性和绝热的流中。图14示出了通过与图13中的相同类型和位置的三个障碍物的爆震启动。图对应于时间：(a)-t＝10、(b)-t＝50、(c)-t＝100，(d)-t＝200、(e)-t＝500、和(f)-t＝1500。域尺寸为600×600，网格点的数量为2500x 2500。如图14(a)所示，弓形冲击可以在障碍物前面形成。快速启动的爆震可通过障碍物便利于波的稳定。然而，一段时间后，这些爆震波彼此分离形成独立结构(图14(b))。不久后，分离的爆震峰融合并变形(图14(c-e))成图13所示的相同最终结构。
图15示出图11-13的反应区的详细结构。图15(a)-崩溃t＝10-对应于图11(c)中的情况。图15的(b)-膨胀t＝150-对应于图12(b)中的情况。图15(c)-爆震被障碍物稳定t＝700-对应于图13(c)中的情况。
示例性爆震发动机100可以包括用于输送氧化剂到发动机的第一入口135，和用于输送燃料到发动机的第二入口140。第一入口135可以流体地连接到包含氧化剂的第一储箱(未示出)，而第二入口140可以流体地连接到包含燃料的第二储箱(未示出)。在一个例子中，如图1所示，第一和第二入口(135、140)可以被定位在发动机的相对侧上。因此，流体和氧化剂可沿着共同轴线145在相反方向上输送到发动机100。第一入口135可以包括第一喷嘴105，而第二入口140可包括第二喷嘴110。随着氧化剂流过第一喷嘴105，超音速流可以在发动机100内实现。类似地，随着燃料流过第二喷嘴110，超音速流可以在发动机100内实现。
示例性爆震发动机100可以包括位于第一与第二入口(135、140)之间并沿着所述共同轴线145的分离器115。分离器115的第一侧可以从共同轴线145径向向外扩散第一氧化剂流。同样地，分离器115的第二侧可以从共同轴线145径向向外扩散第二燃料流。随着所述第一和第二流从共同轴线145径向向外行进，所述流最终行进超出分离器115的外周。在这里，第一和第二流开始混合。流的混合在一定的时间和距离上发生，并且随着它从共同轴线145向外移动，所述混合流变得更加均匀。当充分混合已经发生时，所述混合物的爆震可以发生。
爆震可以在较宽的初始压强条件范围内启动，但是一般应保持一定的温度和化学成分范围，以避免淬熄(quenching)。在1000K的量级上的温度和30-40巴的压强在气态爆震的冲击中是典型的，而2500-3000K以及20-30巴在爆炸产物中是典型的。这些参数的优选范围可以取决于特定设计条件和所需结果。由于爆震发动机100的几何结构，对爆震合适的条件只存在于距共同轴线145的一定距离处。因此，如图2所示，爆震可以发生在具有圆柱形几何形状的爆震区域205内。爆震区域205将包括爆震驻波。
一个典型的爆震区域可具有跟随着反应区的薄的冲击，所述反应区的范围可以在很大程度上取决于混合物的成分和燃烧条件。反应区通常可以是从几分之一毫米到几个毫米或更高。如果设备的尺寸显著超出这些尺度，这样快速的燃烧在最小化几何和损耗因子对燃烧过程的影响方面可以是有利的。随着混合物通过爆震区域205，爆震发生，并且热随着混合物转换为燃烧产物而被释放。热的燃烧产物继续从爆震区域205和共同轴线145径向向外行进，然后穿过并旋转涡轮120，从而产生有用功。随着燃烧产物穿过该涡轮，能量可以从燃烧产物中提取出来，导致其温度和压强的降低。
涡轮120可被构造为独立于发动机自由旋转。在一个例子中，如图1所示，涡轮可分别安装在附接到第一和第二入口(135、140)的第一和第二轴承(125、130)上。可选地，涡轮120可以在一个轴承或多于两个的轴承上转动。在另一实例中，涡轮120可以是未集成到发动机100。例如，如在汽车应用中常见的，涡轮120可安装在从发动机100延伸的排气管中。虽然涡轮120被示出并描述，这不是限制性的；用于从有升高的温度和压强的燃烧产物中产生功的任何类似部件可以被作为代用品。
为了实现爆震驻波，氧化剂和燃料能够以等于所得爆震波速度的速度被输送到爆震区域205。作为一个例子，如果爆震波否则会以2千米/秒的速率在爆震发动机100内的温度和压力下传播，则燃料和氧化剂应以2千米/秒的速率输送到爆震区域205以得到在爆震区域205中的爆震驻波。在一个例子中，所述氧化剂可以在约2千米/秒的超声速度下离开第一喷嘴105，而燃料可以在约2千米/秒的超声速度下离开第二喷嘴110。
图3和图4是示出根据压力或温度表示的稳定爆震波的示例性计算模型。这些模型基于对以70％的氦稀释的化学计量氢-氧混合物的理想气体代表的计算，所述混合物是从直径为10(以在稳定平面爆震中的反应区厚度为 单位)的源以超声速发出的。所述气体最初无化学反应地膨胀和加速，直到它运行成爆震冲击，在这之后，它开始燃烧。如在两个图中示出的，燃烧过程是不均匀的，并导致了典型的多维反应区，其有随着主冲击传播的横向冲击波。图3以巴为单位示出了压强场。所述压强，最初在从源膨胀的过程中减小，经历了越过爆震冲击的显著上升，并随后在流的外部区域中随着燃烧后的气体径向地在高速下喷出而减小。图4以395K为单位示出了流的温度。
在示例性爆震发动机中，燃烧产物可以与旋转叶片相互作用。这种相互作用可以实现几件事情。它可以提供必要的推力以产生有用功。所述相互作用可以通过减少径向膨胀而稳定爆震反应区。爆震可以通过在爆震冲击之后的流中放置障碍物而稳定化。为了模拟旋转叶片的这种稳定效果，一些小的刚性圆形障碍物可以沿着一定半径的环放置。在温度图中，这些障碍物可以清楚地作为小白点看到，其中可见从障碍物散发的漩涡条纹，换句话说，在涡列的头部。在示例性实施例中，几何形状与图中不同的旋转叶片可以在发动机内实现这种稳定性。这样的叶片可以既提供用于将爆震动能转化为发动机的机械能的装置，又提供用于稳定爆震的装置。障碍物的位置可影响稳态圆形爆震的半径，但通过这些障碍物的存在，可以防止爆震波的无限膨胀。稳定的爆燃燃烧能够以在这种几何形状和其他几何形状实现。
一个或多个实施例的细节在附图和说明书中阐述。其他的特征、目的和优点从说明书、附图和权利要求书中显见。虽然已经描述了本发明的多个实施例，应理解的是，可以在不脱离本发明的精神和范围的情况下做出各种修改。还应当理解的是，附图不一定是按比例的，其呈现了本发明的多种特征和基本原理的稍微简化的图示。
参考文献
下列各项通过引用以其整体并入本文。
F.A.Bykovskii,S.A.Zhdan,and E.F.Vedernikov.Continuous spindetonations.Journal of Propulsion and Power,22(6):1204–1216,2006.
F.A.Bykovskii,V.V.Mitrofanov,and E.F.Vedernikov.Continuous detonationcombustion of fuel-air mixtures.Combustion,Explosion and Shock Waves,33(3):344–353,1997.
A.Dadone.Symmetry techniques for the numerical solution of the 2DEuler equations at impermeable boundaries.International Journal for NumericalMethods in Fluids,28(7):1093–1108,1998.
W.Uber den detonationvorgang in gasen.Annalen der Physik,43(6/7):421–428,1943.W.Fickett and W.C.Davis.Detonation:theory andexperiment.Dover Publications,2011.
S.Gottlieb and C.W.Shu.Total variation diminishing runge-kuttaschemes.Mathematics of Computation,67(221):73–85,1998.
K.Kailasanath.Review of propulsion applications of detonationwaves.AIAA journal,38(9):1698–1708,2000.
J.H.S.Lee.The Detonation Phenomenon.Cambridge University Press,2008.
X.-D.Liu,S.Osher,and T.Chan.Weighted essentially non-oscillatoryschemes.Journal of Computational Physics,115(1):200–212,1994.
J.A.Nicholls and E.K.Dabora.Recent results on standing detonationwaves.In Symposium(International)on Combustion,volume 8,pages 644–655.Elsevier,1961.
G.D.Roy,S.M.Frolov,A.A.Borisov,and D.W.Netzer.Pulse detonationpropulsion:challenges,current status,and future prospects.Progress in Energyand Combustion Science,30(6):545–672,2004.
P.M.Rubins and R.C.Bauer.Review of shock-induced supersoniccombustion research and hypersonic applications.Journal of Propulsion andPower,10(5):593–601,1994.
R.I.Soloukhin.Shock waves and detonations in gases.Mono BookCorporation,1966.
B.D.Taylor,A.R.Kasimov,and D.S.Stewart.Mode selection in unstabletwo-dimensional detonations.Combust.Theory Model,13(6):973–992,2009.
A.A.Vasil’ev,V.I.Zvegintsev,and D.G.Nalivaichenko.Detonation waves ina reactive supersonic flow.Combustion,Explosion and Shock Waves,42(5):568–581,2006.
B.V.Voitsekhovskii,V.V.Mitrofanov,and M.Y.Topchian.The Structure ofDetonation Front in Gases.Report FTD-MTD-64-527.Foreign Technology Division,Wright Patterson Air Force Base,OH(AD 633-821).,1966.
J.von Neumann.Theory of detonation waves.Office of ScientificResearch and Development,Report 549.Technical report,National DefenseResearch Committee Div.B,1942.
P.Wolański.Detonative propulsion.Proceedings of the CombustionInstitute,2012.
Y.B.Zel’dovich.On the theory of propagation of detonation in gaseoussystems.J.Exp.Theor.Phys.,10(5):542–569,1940.
F.Zhang,editor.Shock Waves Science and Technology Library,Vol.6:Detonation Dynamics,volume 6.Springer,2012.
T.A.Zhuravskaya and V.A.Levin.Investigation of certain techniques forstabilizing detonation waves in a supersonic flow.Fluid Dynamics,47(6):793–801,2012.

Claims (17)

1.一种爆震发动机，包括：
第一入口，其具有流体地连接到第一储箱的第一端部和流体地连接到爆震发动机的第二端部；
第二入口，其具有流体地连接到第二储箱的第一端部和与第一入口相对地流体连接到所述爆震发动机的第二端部，其中，所述第一和第二入口可与共同轴线对准；
连接到所述第一入口的第一喷嘴；
连接到所述第二入口的第二喷嘴；以及
放置在所述第一入口的第二端部与所述第二入口的第二端部之间并且沿共同轴线放置的分离器，
其中，所述爆震发动机还包括：可旋转地附接到所述爆震发动机的至少一个涡轮，所述至少一个涡轮包括构造为在所述爆震发动机中稳定爆震的叶片。
2.根据权利要求1所述的爆震发动机，还包括：构造为在所述爆震发动机中稳定爆震的障碍物。
3.根据权利要求1所述的爆震发动机，其中，所述第一储箱构造为接收氧化剂。
4.根据权利要求1所述的爆震发动机，其中，所述第二储箱构造为接收燃料。
5.根据权利要求1所述的爆震发动机，其中，所述第一喷嘴靠近所述第一入口的第二端部。
6.根据权利要求1所述的爆震发动机，其中，所述第二喷嘴靠近所述第二入口的第二端部。
7.根据权利要求1所述的爆震发动机，其中，所述第一喷嘴布置在所述第一入口内。
8.根据权利要求1所述的爆震发动机，其中，所述第二喷嘴布置在所述第二入口内。
9.一种用于引爆燃料与氧化剂的混合物的方法，所述方法包括：
将氧化剂通过第一入口引入到爆震发动机；
将燃料通过第二入口引入到爆震发动机，其中，所述第一和第二入口与共同轴线对准；
将所述氧化剂加速通过第一喷嘴；
将所述燃料加速通过第二喷嘴；
靠着分离器的第一侧并且从所述共同轴线径向向外地引导所述氧化剂；
靠着分离器的与所述第一侧相对的第二侧并且从所述共同轴线径向向外地引导所述燃料；以及
在距共同轴线一定距离处的圆柱形爆震区域内引爆燃料与氧化剂的混合物。
10.根据权利要求9所述的方法，其中，所述第一入口具有流体地连接到第一储箱的第一端部和流体地连接到所述爆震发动机的第二端部。
11.根据权利要求10所述的方法，其中，所述第二入口具有流体地连接到第二储箱的第一端部和流体地连接到所述爆震发动机的第二端部。
12.根据权利要求9或11所述的方法，其中，将氧化剂加速通过所述第一喷嘴产生超音速流，并且其中，将燃料加速通过所述第二喷嘴产生超音速流。
13.根据权利要求12所述的方法，还包括：将由引爆所述混合物产生的燃烧产物膨胀通过涡轮以做功。
14.根据权利要求9所述的方法，还包括：提供障碍物以抑制燃烧产物的膨胀。
15.根据权利要求9所述的方法，还包括：在所述爆震发动机中稳定混合物的爆震。
16.根据权利要求15所述的方法，其中，障碍物构造在所述爆震发动机中以实现爆震的稳定化。
17.根据权利要求16所述的方法，其中，所述障碍物是涡轮叶片。
CN201380063447.7A 2012-10-12 2013-10-11 爆震驻波发动机 Expired - Fee Related CN104870780B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201261712972P 2012-10-12 2012-10-12
US61/712,972 2012-10-12
PCT/IB2013/003257 WO2014114980A2 (en) 2012-10-12 2013-10-11 Standing detonation wave engine

Publications (2)

Publication Number Publication Date
CN104870780A CN104870780A (zh) 2015-08-26
CN104870780B true CN104870780B (zh) 2017-03-08

Family

ID=50896333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380063447.7A Expired - Fee Related CN104870780B (zh) 2012-10-12 2013-10-11 爆震驻波发动机

Country Status (6)

Country Link
US (1) US9841192B2 (zh)
EP (1) EP2906795B1 (zh)
JP (1) JP6207617B2 (zh)
CN (1) CN104870780B (zh)
RU (1) RU2616143C2 (zh)
WO (1) WO2014114980A2 (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11149954B2 (en) 2017-10-27 2021-10-19 General Electric Company Multi-can annular rotating detonation combustor
RU2689056C1 (ru) * 2018-09-04 2019-05-23 Общество с ограниченной ответственностью "Новые физические принципы" Способ и устройство для газодинамического разгона массивных тел до высокой скорости
CN112562793B (zh) * 2020-12-10 2023-01-03 北京理工大学 一种针对燃料爆震燃烧的两步反应模型计算方法
DE102021004141A1 (de) 2021-03-31 2022-10-06 Mathias Herrmann Angepasstes Verfahrenskonzept und Leistungskonzept für Triebwerke (z.B. Raketen, Überschallraketenbrennkammern / Düsentriebwerke), luftatmende Antriebe (z.B. Unterschall-Staustrahltriebwerke, Ramjets, Scramjets, Dualmode, Pulsejets, Detonationstriebwerke, Raketen-Staustrahl- Antriebe), Turbopumpen bzw. Düsen (z.B. Glockendüsen, Aerospikes)
CN113153570B (zh) * 2021-05-27 2023-09-22 中国航空发动机研究院 一种脉冲爆震管性能计算方法和装置
CN114427689B (zh) * 2022-01-20 2024-07-12 中国空气动力研究与发展中心空天技术研究所 一种可实现超声速流场观测的圆盘形旋转爆震燃烧室

Family Cites Families (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US364866A (en) * 1887-06-14 Reaction-wheel
US1260015A (en) * 1917-08-13 1918-03-19 Clifton F Noakes Turbine hydrogas-motor.
US2746529A (en) * 1950-11-17 1956-05-22 Stewart Warner Corp Resonant combustion device for heating engines and similar structures
US2942412A (en) * 1952-09-30 1960-06-28 Curtiss Wright Corp Pulse detonation jet propulsion
US2796735A (en) * 1956-08-20 1957-06-25 Jr Albert G Bodine Acoustic jet engine with flow deflection fluid pumping characteristics
US3061337A (en) * 1958-03-05 1962-10-30 Perfecting Service Company Rotary fluid connector with noncommunicating passageways
US3217491A (en) * 1965-02-08 1965-11-16 Thomas N Diehl Method of producing energy in a reaction engine
US3336754A (en) * 1966-03-21 1967-08-22 Oswald H Lange Continuous detonation reaction engine
CH540429A (de) * 1970-04-16 1973-08-15 Andersson Louis Einrichtung zur Erzeugung von Treibgas
SU890030A1 (ru) * 1980-01-04 1981-12-15 Сибирский металлургический институт им.Серго Орджоникидзе Камера пульсирующего горени с резонансной трубой
US4741154A (en) * 1982-03-26 1988-05-03 The United States Of America As Represented By The Secretary Of The Navy Rotary detonation engine
JPS6093211A (ja) * 1983-10-28 1985-05-25 Toshiba Corp パルス燃焼装置
US4903480A (en) * 1988-09-16 1990-02-27 General Electric Company Hypersonic scramjet engine fuel injector
JP2726487B2 (ja) * 1989-03-31 1998-03-11 株式会社東芝 パルスバーナ
US5123835A (en) * 1991-03-04 1992-06-23 The United States Of America As Represented By The United States Department Of Energy Pulse combustor with controllable oscillations
US5138831A (en) * 1991-03-07 1992-08-18 Cowan Sr Howard H Air cooled rotary combustion engine
US6062018A (en) * 1993-04-14 2000-05-16 Adroit Systems, Inc. Pulse detonation electrical power generation apparatus with water injection
US6464490B1 (en) * 1998-08-31 2002-10-15 Clean Energy Combustion Systems, Inc. Circular pulsating combustors
US6938588B2 (en) * 1999-11-12 2005-09-06 Sarcos Investments, Lc Controllable combustion method and device
GB0128811D0 (en) * 2001-12-03 2002-01-23 Robinson Nicholas P Propulsion device-sky-wave
US6584765B1 (en) * 2001-12-21 2003-07-01 United Technologies Corporation Pulse detonation engine having an aerodynamic valve
JP3981743B2 (ja) * 2002-04-19 2007-09-26 北海道ティー・エル・オー株式会社 定常デトネーション燃焼器および定常デトネーション波生成方法
JP4256820B2 (ja) * 2004-06-29 2009-04-22 三菱重工業株式会社 デトネーションエンジンおよびこれを備えた飛行体
US7475914B2 (en) * 2004-09-14 2009-01-13 Pospisil Mark G Rotary union for directing a heat exchange fluid
RU48368U1 (ru) * 2005-03-18 2005-10-10 Открытое акционерное общество Опытно-конструкторское бюро "Сокол" Пульсирующий воздушно-реактивный двигатель
JP4912659B2 (ja) * 2005-10-14 2012-04-11 国立大学法人 筑波大学 衝撃波連続発生装置
JP5166698B2 (ja) 2006-01-18 2013-03-21 帝人化成株式会社 押出成形用シートの原料として用いるペレットの製造方法
US20090266047A1 (en) * 2007-11-15 2009-10-29 General Electric Company Multi-tube, can-annular pulse detonation combustor based engine with tangentially and longitudinally angled pulse detonation combustors
GB0822676D0 (en) * 2008-12-12 2009-01-21 Rolls Royce Plc A gas turbine engine
GB2468515B (en) * 2009-03-09 2011-02-02 Salahuddin Rasheedsedqi Saeed Clean aerospace engine
WO2011155248A1 (ja) * 2010-06-10 2011-12-15 学校法人早稲田大学 エンジン
US20120192630A1 (en) * 2011-01-28 2012-08-02 General Electric Company Pulse Detonation Turbine Engine Using Turbine Shaft Speed for Monitoring Combustor Tube Operation
JP2013044455A (ja) * 2011-08-23 2013-03-04 Yoshiyuki Nakata 内燃機関
JP2014171406A (ja) 2013-03-06 2014-09-22 Chugoku Electric Power Co Inc:The 不要物回収装置

Also Published As

Publication number Publication date
CN104870780A (zh) 2015-08-26
WO2014114980A2 (en) 2014-07-31
WO2014114980A3 (en) 2015-01-08
RU2015117637A (ru) 2016-12-10
EP2906795A2 (en) 2015-08-19
US20150285506A1 (en) 2015-10-08
RU2616143C2 (ru) 2017-04-12
EP2906795B1 (en) 2019-05-01
US9841192B2 (en) 2017-12-12
JP6207617B2 (ja) 2017-10-04
JP2015531456A (ja) 2015-11-02

Similar Documents

Publication Publication Date Title
CN104870780B (zh) 爆震驻波发动机
Jiang et al. Criteria for hypersonic airbreathing propulsion and its experimental verification
Gaillard et al. Numerical simulation of a rotating detonation with a realistic injector designed for separate supply of gaseous hydrogen and oxygen
Evans et al. Influence of chemical kinetics and unmixedness on burning in supersonic hydrogen flames
Davidenko et al. Numerical study of the continuous detonation wave rocket engine
Meng et al. Numerical investigation of the effect of inlet mass flow rates on H2/air non-premixed rotating detonation wave
Huang et al. Flow-field analysis of a typical hydrogen-fueled dual-mode scramjet combustor
Shi et al. Numerical investigation on the forced oscillation of shock train in hypersonic inlet with translating cowl
Candon et al. Numerical analysis and design optimization of supersonic after-burning with strut fuel injectors for scramjet engines
Teasley et al. A review towards the design optimization of high performance additively manufactured rotating detonation rocket engine injectors
Cai et al. Adaptive mesh refinement–based numerical simulation of detonation initiation in supersonic combustible mixtures using a hot jet
Fotia et al. Experimental study of performance scaling in rotating detonation engines operated on hydrogen and gaseous hydrocarbon fuel
Wu et al. Hysteretic behaviors in counter-flow thrust vector control
Cai et al. Detonation simulations in supersonic combustible mixtures with nonuniform species
Luo et al. Powder fuel transport process and mixing characteristics in cavity-based supersonic combustor with different injection schemes
Saeid et al. Effect of diffusion time on the mechanism of deflagration to detonation transition in an inhomogeneous mixture of hydrogen-air
Ying et al. Initiation characteristics of oblique detonation waves from a finite wedge under argon dilution
Vinogradov et al. A review of fuel pre-injection in supersonic, chemically reacting flows
Radhakrishnan et al. A detailed modelling on spray atomisation and combustion of LOX/GCH4 in variable area pintle injector rocket engines
Cai et al. Parametric investigation of secondary injection in post-chamber on combustion performance for hybrid rocket motor
Shi et al. Numerical investigation on behaviors of shock train in a hypersonic inlet with translating cowl
Liang et al. Comprehensive modeling of a liquid rocket combustion chamber
Sosa et al. Exploration of nozzle flow circumferential attenuation and efficient expansion for rotating detonation rocket engines
Brophy et al. Detonation initiation improvements using swept-ramp obstacles
Sun et al. Prediction of detonation initiation in an air-breathing pulse detonation combustor utilizing a simplified ethylene-air reaction mechanism

Legal Events

Date Code Title Description
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170308

CF01 Termination of patent right due to non-payment of annual fee