CN110135060A - 方形多孔热防护材料的优化设计方法 - Google Patents
方形多孔热防护材料的优化设计方法 Download PDFInfo
- Publication number
- CN110135060A CN110135060A CN201910399186.5A CN201910399186A CN110135060A CN 110135060 A CN110135060 A CN 110135060A CN 201910399186 A CN201910399186 A CN 201910399186A CN 110135060 A CN110135060 A CN 110135060A
- Authority
- CN
- China
- Prior art keywords
- rectangular
- sound
- protective materials
- wave
- thermally protective
- 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.)
- Pending
Links
- 239000000463 material Substances 0.000 title claims abstract description 81
- 230000001681 protective effect Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000013461 design Methods 0.000 title claims abstract description 29
- 238000005457 optimization Methods 0.000 claims abstract description 30
- 230000007704 transition Effects 0.000 claims abstract description 17
- 239000012530 fluid Substances 0.000 claims abstract description 8
- 239000011148 porous material Substances 0.000 claims abstract description 8
- 239000002245 particle Substances 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 7
- 230000000737 periodic effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000002829 reductive effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract description 4
- 238000011161 development Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公开的方形多孔热防护材料的优化设计方法，涉及抑制边界层内Mack第二模态不稳定波的方形多孔热防护材料的优化设计方法，属于航空航天领域。本发明实现方法为：定义方形微孔形状参数与背景空气介质参数。建立声场模型，确定声场模型中的声压反射系数。声场模型包括入射波模型、反射波模型和孔内波模型。通过数值方法优化使0,0阶的声压反射系数|R00|最小，确定方形微孔几何参数，得到满足热防护目的的多孔热防护材料。将优化得到的满足热防护的目的的多孔热防护材料应用于热防护领域，使优化的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速来流中的热载荷保持在层流区的热载荷转态，最终实现热防护的目的。
Description
技术领域
本发明涉及一种新型的方形多孔热防护材料的优化设计方法，尤其涉及抑制边界层内Mack第二模态不稳定波的方形多孔热防护材料的优化设计方法，属于航空航天领域。
背景技术
随着的航空航天事业的蓬勃发展，其对高速飞行器热防护提出了越来越高的要求。热防护的主要技术思路之一是通过抑制边界层转捩来减少飞行器的热载荷，这是因为边界层流态对高超飞行器气动热环境的影响显著，湍流区气动加热远大于层流区，而热流峰值一般出现在流动转捩区，高超声速时可达层流热载荷的3倍，即对于热载荷有：层流区＜湍流区＜转捩区。因此，飞行器热防护对如何有效抑制转捩有着迫切的需求，它也是目前国内外高超声速领域的研究热点。
对高超声速飞行器转捩抑制的研究，主要集中在如何更为有效地抑制边界层内Mack第二模态的发展。微型多孔覆盖层是在壁面加工一系列远小于扰动波长的缝/孔，在不影响主流的前提下，通过孔内粘性来耗散吸收扰动波的能量，从而达到抑制Mack第二模态发展的目的。这项技术已经在高超风洞实验中证实了能够显著延迟转捩，并作为技术储备授权于美国Boeing公司。对多孔覆盖层的应用，以及更进一步的参数优化设计，首先需要完成多孔材料的微结构-声场-流场耦合理论模型的研究，取得其在预设的来流条件下最小的吸声系数，以抑制Mack第二模态的发展，并获得良好的热防护效果。
2001年Fedorov等人提出了细长孔的声学特性，给出单孔对流场和声场的作用模型：v′＝Ap′。其中v′为壁面法向脉动速度，p′为壁面声场压力，A为导纳。2005年Kozlov等人给出了不同截面形状的微孔导纳表达式。需要指出的是，上述理论模型都基于单孔声学特性推导，并未计及开孔处实体壁面引起的高阶模态散射，也未计及孔与孔之间的干扰，从而降低了模型的预测精度。
发明内容
本发明公开的方形多孔热防护材料的优化设计方法要解决的技术问题是：通过所述方形多孔热防护材料的优化设计方法对方形多孔热防护材料进行优化设计，进而有效抑制方形多孔热防护材料表面边界层转捩的发生，使设计的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速时的热载荷保持在层流区的热载荷转态，通过方形多孔热防护材料实现热防护的目的。
本发明目的是通过下述技术方案实现的。
本发明公开的方形多孔热防护材料的优化设计方法，定义方形微孔形状参数与背景空气介质参数。建立声场模型，并确定声场模型中的声压反射系数。所述声场模型包括入射波模型、反射波模型和孔内波模型。通过数值方法优化使0,0阶的声压反射系数|R00|最小，确定方形微孔几何参数，得到满足热防护目的的多孔热防护材料。将优化得到的满足热防护的目的的多孔热防护材料应用于热防护领域，使优化的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速来流中的热载荷保持在层流区的热载荷转态，通过方形多孔热防护材料最终实现热防护的目的。
本发明公开的方形多孔热防护材料的优化设计方法，包括如下步骤：
步骤一：定义方形微孔形状参数与背景空气介质参数。
所述方形多孔热防护材料由周期规律分布的多个方形微孔组成，定义孔深方向为y方向，与孔深方向垂直的材料表面平面定义为xz平面，xz平面上以s为周期规律地分布着边长为2b、孔深为h的方形微孔。无量纲几何参数孔隙率φ由上述参数表示为φ＝4b2/s2，无量纲几何参数宽深比为Ar由上述参数表示为Ar＝2b/h。背景介质的密度为ρ0，声波在此介质中的传播速度为c0。
步骤二：建立声场模型，并确定声场模型中的声压反射系数。所述声场模型包括入射波模型、反射波模型和孔内波模型。
步骤2.1：建立入射波模型。
背景介质中存在与时间相关的微幅扰动e-jωt，则在声场中传播的任意一束平面入射波表达为
其中，pi为入射压力，vy,i为质点的y方向速度；为垂直动量；为平行于平面的动量；波数k0＝ω/c0；ω为角频率；
步骤2.2：建立反射波模型。
对于m,n阶反射波，其压力场和y方向上的质点速度分别为
其中 Rmn为m,n阶反射系数，特殊的0,0阶反射，即为镜面反射。
步骤2.3：建立孔内波模型。
声波在方形微孔中的传播，会因为热传导和粘性的作用而产生损耗。而由于方形微孔的边长远小于波长，即2b＜＜λacs，所以在孔内波长的限制下，声波以基础模态波为主。因此方形微孔中的声压和质点y方向速度表示为
式中的有效密度ρh、压缩系数Ch、波数kh都是与角频率ω相关的复数参数，并分别有
ρh(ω)＝ρ0/Ψv，
其中气体绝热指数γ＝Cp/Cv，Cp为定压比热，Cv为定容比热；
其中μ为动力粘度，κ为导热系数。
步骤2.4：通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定声场模型中的声压反射系数表达式，并通过声压反射系数进一步确定0,0阶声压反射系数R00表达式。
平面声波入射到方形微孔表面产生的声场环境与阻抗特性通过平面波展开法获得，即通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定第m,n阶衍射波的声压反射系数如下。
式中，δmn,00为Kronecker函数，当m＝n＝0时，δmn,00＝1；当m,n不都为零时，δmn,00＝0。Smn为m,n阶衍射模态与孔内主导模态的重叠积分
Mack第二模态不稳定波简化为垂直入射，且孔周期远小于入射波波长，即s＜＜λacs，则有效声导率A表示为
其中R00为0,0阶的声压反射系数，即镜面声压反射系数。
步骤三：通过数值方法优化使0,0阶的声压反射系数|R00|最小，确定方形微孔几何参数，得到满足热防护目的的多孔热防护材料。
在给定的来流条件下，通过数值方法优化使方形多孔热防护材料的0,0阶声压反射系数|R00|最小，确定方形微孔几何参数最佳的孔隙率φ、宽深比Ar，即完成方形多孔热防护材料的优化设计。
具体实现方法为：优化目标为在优化区间内使0,0阶的声压反射系数|R00|最小，约束条件为：将方形微孔的几何参数限定在材料实际工程加工中的能够实现范围内，且入射波频率不超出精度极限
作为优选，所述优化区间选为下述三个优化区间分别为loop 1：0.2≤φ≤0.8，loop 2：0.12≤Ar≤0.6，loop 3：所述三个优化区间目的是用于满足约束条件。
通过使0,0阶的声压反射系数|R00|最小的数值方法优化得到的方形微孔几何参数最佳孔隙率φ、宽深比Ar，通过优化后的方形微孔几何参数最佳孔隙率φ、宽深比Ar实现的多孔热防护材料，能够有效抑制壁面边界层内Mack第二模态的发展，进而有效抑制方形多孔热防护材料表面边界层转捩的发生，使设计的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速来流中的热载荷保持在层流区的热载荷转态，通过方形多孔热防护材料最终实现热防护的目的。
还包括步骤四：将步骤三优化得到的满足热防护的目的的多孔热防护材料应用于热防护领域，使优化的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速来流中的热载荷保持在层流区的热载荷转态，通过方形多孔热防护材料最终实现热防护的目的。
有益效果：
1、本发明公开的方形多孔热防护材料的优化设计方法，建立的反射波模型中计及高阶模态衍射波，能够提高声压反射系数|R00|的预测精度。
2、本发明公开的方形多孔热防护材料的优化设计方法，在约束条件中将方形微孔的几何参数限定在材料实际工程加工中能够实现的范围内，进而保证方形多孔热防护材料易于在材料实际工程加工。
3、本发明公开的方形多孔热防护材料的优化设计方法，在给定的来流条件下用数值方法优化确定方形微孔几何参数，使方形多孔热防护材料的0,0阶声压反射系数|R00|最小，实现有效抑制壁面边界层内的Mack第二模态的发展和转捩的发生，进而有效抑制扰动强度和脉动压力幅值，最终达到热防护的目的。
4、本发明公开的方形多孔热防护材料的优化设计方法，方形多孔热防护材料技术简单、易于加工，能够在加工过程中节约成本；方形多孔热防护材料结构简单、质量轻，在保证热防护效果的同时为飞行器减轻结构重量。
附图说明
图1是本发明公开的方形多孔热防护材料的优化设计方法的流程图；
图2是方形多孔热防护材料的物理模型图；
图3是φ＝0.3，Ar＝0.5的方形多孔热防护材料的理论模型预测反射系数对比图；
图4是φ＝0.3，Ar＝0.16的方形多孔热防护材料的理论模型预测反射系数对比图；
图5是φ＝1，Ar＝0.5的方形多孔热防护材料的理论模型预测反射系数对比图；
图6是φ＝1，Ar＝0.16的方形多孔热防护材料的理论模型预测反射系数对比图；
图7是平板流动的压力脉动云图；
图8是平板流动壁面的压力幅值对比图。
具体实施方式
为了更好的说明本发明的目的和优点，下面结合附图和实例对发明内容做进一步说明。
实施例1：
如图1所示，本实施例公开的方形多孔热防护材料的优化设计方法，具体实现步骤如下：
步骤一：定义方形微孔形状参数与背景空气介质参数。
所述方形多孔热防护材料如图2是由周期规律分布的多个方形微孔组成，定义孔深方向为y方向，与孔深方向垂直的材料表面平面定义为xz平面，xz平面上以s为周期规律地分布着边长为2b、孔深为h的方形微孔。无量纲几何参数孔隙率φ由上述参数表示为φ＝4b2/s2，无量纲几何参数宽深比为Ar由上述参数表示为Ar＝2b/h。背景介质空气的密度为ρ0，声波在此介质中的传播速度为c0。
步骤二：建立声场模型，并确定声场模型中的声压反射系数。所述声场模型包括入射波模型、反射波模型和孔内波模型。
步骤2.1：建立入射波模型。
背景介质中存在与时间相关的微幅扰动e-jωt，则在声场中传播的任意一束平面入射波表达为
其中，pi为入射压力，vy,i为质点的y方向速度；为垂直动量；为平行于平面的动量；波数k0＝ω/c0；ω为角频率；
步骤2.2：建立反射波模型。
对于m,n阶反射波，其压力场和y方向上的质点速度分别为
其中 Rmn为m,n阶反射系数，特殊的0,0阶反射，即为镜面反射。
步骤2.3：建立孔内波模型。
声波在方形微孔中的传播，会因为热传导和粘性的作用而产生损耗。而由于方形微孔的边长远小于波长，即2b＜＜λacs，所以在孔内波长的限制下，声波以基础模态波为主。因此方形微孔中的声压和质点y方向速度可以表示为
式中的有效密度ρh、压缩系数Ch、波数kh都是与角频率ω相关的复数参数，并分别有
ρh(ω)＝ρ0/Ψv，
其中气体绝热指数γ＝Cp/Cv，Cp为定压比热，Cv为定容比热；
其中μ为动力粘度，κ为导热系数。
步骤2.4：通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定声场模型中的声压反射系数表达式，并通过声压反射系数进一步确定0,0阶声压反射系数R00表达式。
平面声波入射到方形微孔表面产生的声场环境与阻抗特性可通过平面波展开法获得，即通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定第m,n阶衍射波的声压反射系数如下。
式中，δmn,00为Kronecker函数，当m＝n＝0时，δmn,00＝1；当m,n不都为零时，δmn,00＝0。Smn为m,n阶衍射模态与孔内主导模态的重叠积分。
Mack第二模态不稳定波简化为垂直入射，且孔周期远小于入射波波长，即s＜＜λacs，则有效声导率A表示为
其中R00为0,0阶的声压反射系数，即镜面声压反射系数。
步骤三：验证反射系数R00预测的准确性
图3至图6展示了在温度均匀的静止环境中，单频声波传入多孔材料表面的反射系数对比图。图中包含了本专利所提出模型，Fedorov提出的单孔无耦合模型，和有限元求解器(COMSOL )，共三组计算数据。图3额外提供了一组COMSOL网格加密仿真数据以验证网格无关性，该组数据在COMSOL仿真中的网格数量是另外一组的8.6倍。方形多孔材料参数选取孔隙率φ＝0.16,0.5，宽深比Ar＝0.3,1；将入射波频率f无量纲化后有无量纲入射频率facs＝f H/c0。可以发现，本专利提出的模型与COMSOL结果吻合较好，而Fedorov提出的模型由于未计及高阶衍射模态，预测的反射系数在高频范围内出现漂移。
步骤四：通过数值方法优化使0,0阶的声压反射系数|R00|最小，确定方形微孔几何参数，得到满足热防护的目的的多孔热防护材料。
以某给定的高超声速流动环境为例进行说明。在这个给定的流动环境下，背景介质空气密度ρ0＝0.005588337kg/m3，温度T＝293K，动力粘度μ＝1.8127×10-5N·s/m2，使用数值方法对方形微孔的几何参数进行优化，目的为在优化区间内使方形多孔热防护材料的0,0阶声压反射系数|R00|最小。优化区间为：loop 1：0.2≤φ≤0.8，loop 2：0.12≤Ar≤0.6，loop 3：以上范围参考了实际中的工程可实现性以及模型的精度极限：为模型的共振频率，是facs的频率上限。多孔微结构的几何参数优化结果如下表1所示，得到最小|R00|＝2.8×10-5。
表1某给定的来流条件下，方形微孔的几何参数优化结果
步骤五：验证优化后的方形多孔热防护材料抑制Mack第二模态的效果。
将上述表1中经过几何参数优化设计的方形多孔热防护材料安装在长度为Lref＝0.2米的平板的后半部，而平板的前半部则未安装热防护材料，即前半部为光滑平板，并以此为平板流动验证算例，假定Mach 6.0高超声速来流，使用数值模拟获得平板的流场脉动压力分布。对比图7(a)和图7(b)，可以发现方形多孔热防护材料能够有效抑制扰动强度，实现抑制Mack第二模态的目的。如图8所示，将普通光滑材料与方形多孔热防护材料的壁面脉动压力幅值进行对比，可以观察到热防护材料(即平板后半部)对脉动压力幅值和转捩有着非常明显的抑制作用。
通过图7和图8可以综合判断，方形多孔热防护材料有效抑制壁面边界层内的Mack第二模态的发展，即有效抑制转捩的发生，令方形多孔热防护材料有效规避了热流峰值一般出现在的流动转捩区，使高超声速来流中的热载荷保持如在层流区的热载荷转态，最终达到热防护的目的。
以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.方形多孔热防护材料的优化设计方法，其特征在于：包括如下步骤：
步骤一：定义方形微孔形状参数与背景空气介质参数；
步骤二：建立声场模型，并确定声场模型中的声压反射系数；所述声场模型包括入射波模型、反射波模型和孔内波模型；
步骤三：通过数值方法优化使0,0阶的声压反射系数|R00|最小，确定方形微孔几何参数，得到满足热防护目的的多孔热防护材料。
2.如权利要求1所述的方形多孔热防护材料的优化设计方法，其特征在于：还包括步骤四：将步骤三优化得到的满足热防护的目的的多孔热防护材料应用于热防护领域，使优化的方形多孔热防护材料有效规避流动转捩区的热流峰值，使高超声速来流中的热载荷保持在层流区的热载荷转态，通过方形多孔热防护材料最终实现热防护的目的。
3.如权利要求1或2所述的方形多孔热防护材料的优化设计方法，其特征在于：步骤一实现方法为，
所述方形多孔热防护材料由周期规律分布的多个方形微孔组成，定义孔深方向为y方向，与孔深方向垂直的材料表面平面定义为xz平面，xz平面上以s为周期规律地分布着边长为2b、孔深为h的方形微孔；无量纲几何参数孔隙率φ由上述参数表示为φ＝4b2/s2，无量纲几何参数宽深比为Ar由上述参数表示为Ar＝2b/h；背景介质的密度为ρ0，声波在此介质中的传播速度为c0。
4.如权利要求3所述的方形多孔热防护材料的优化设计方法，其特征在于：步骤二实现方法为，
步骤2.1：建立入射波模型；
背景介质中存在与时间相关的微幅扰动e-jωt，则在声场中传播的任意一束平面入射波表达为
其中，pi为入射压力，vy,i为质点的y方向速度；为垂直动量；为平行于平面的动量；波数k0＝ω/c0；ω为角频率；
步骤2.2：建立反射波模型；
对于m,n阶反射波，其压力场和y方向上的质点速度分别为
其中m,n＝-∞,...,-1,0,1,...,+∞；Rmn为m,n阶反射系数，特殊的0,0阶反射，即为镜面反射；
步骤2.3：建立孔内波模型；
声波在方形微孔中的传播，会因为热传导和粘性的作用而产生损耗；而由于方形微孔的边长远小于波长，即2b＜＜λacs，所以在孔内波长的限制下，声波以基础模态波为主；因此方形微孔中的声压和质点y方向速度表示为
式中的有效密度ρh、压缩系数Ch、波数kh都是与角频率ω相关的复数参数，并分别有
ρh(ω)＝ρ0/Ψv，
其中气体绝热指数γ＝Cp/Cv，Cp为定压比热，Cv为定容比热；
其中μ为动力粘度，κ为导热系数；
步骤2.4：通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定声场模型中的声压反射系数表达式，并通过声压反射系数进一步确定0,0阶声压反射系数R00表达式；
平面声波入射到方形微孔表面产生的声场环境与阻抗特性通过平面波展开法获得，即通过建立的入射波模型、反射波模型和孔内波模型声场模型，推导确定第m,n阶衍射波的声压反射系数如下；
式中，δmn,00为Kronecker函数，当m＝n＝0时，δmn,00＝1；当m,n不都为零时，δmn,00＝0；Smn为m,n阶衍射模态与孔内主导模态的重叠积分
Mack第二模态不稳定波简化为垂直入射，且孔周期远小于入射波波长，即s＜＜λacs，则有效声导率A表示为
其中R00为0,0阶的声压反射系数，即镜面声压反射系数；
5.如权利要求4所述的方形多孔热防护材料的优化设计方法，其特征在于：步骤三实现方法为，
在给定的来流条件下，通过数值方法优化使方形多孔热防护材料的0,0阶声压反射系数|R00|最小，确定方形微孔几何参数最佳的孔隙率φ、宽深比Ar，即完成方形多孔热防护材料的优化设计。
6.如权利要求5所述的方形多孔热防护材料的优化设计方法，其特征在于：步骤三具体实现方法为，
优化目标为在优化区间内使0,0阶的声压反射系数|R00|最小，约束条件为：将方形微孔的几何参数限定在材料实际工程加工中的能够实现范围内，且入射波频率不超出精度极限
7.如权利要求6所述的方形多孔热防护材料的优化设计方法，其特征在于：所述优化区间选为下述三个优化区间分别为loop 1：0.2≤φ≤0.8，loop 2：0.12≤Ar≤0.6，loop 3：所述三个优化区间目的是用于满足约束条件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910399186.5A CN110135060A (zh) | 2019-05-14 | 2019-05-14 | 方形多孔热防护材料的优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910399186.5A CN110135060A (zh) | 2019-05-14 | 2019-05-14 | 方形多孔热防护材料的优化设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110135060A true CN110135060A (zh) | 2019-08-16 |
Family
ID=67573941
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910399186.5A Pending CN110135060A (zh) | 2019-05-14 | 2019-05-14 | 方形多孔热防护材料的优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110135060A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651897A (zh) * | 2020-06-18 | 2020-09-11 | 北京理工大学 | 一种能够抑制高超声速边界层转捩的大尺寸缝隙微结构 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
| CN104198584A (zh) * | 2014-09-18 | 2014-12-10 | 合肥工业大学 | 一种获得圆形开口声传递率及声传递损失的方法 |
-
2019
- 2019-05-14 CN CN201910399186.5A patent/CN110135060A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104198032A (zh) * | 2014-08-14 | 2014-12-10 | 合肥工业大学 | 一种矩形开口声传递率及声传递损失计算方法 |
| CN104198584A (zh) * | 2014-09-18 | 2014-12-10 | 合肥工业大学 | 一种获得圆形开口声传递率及声传递损失的方法 |
Non-Patent Citations (6)
| Title |
|---|
| ALEXANDER V. FEDOROV 等: "Stabilization of Hypersonic Boundary Layers by Porous Coatings", 《AIAA JOURNAL》 * |
| R. ZHAO 等: "Theoretical Modeling and Optimization of Porous Coating for Hypersonic Laminar Flow Control", 《AIAA JOURNAL》 * |
| TUO LIU 等: "Inherent losses induced absorptive acoustic rainbow trapping with a gradient metasurface", 《JOURNAL OF APPLIED PHYSICS》 * |
| VIOLA WARTEMANN 等: "Numerical Investigation of Hypersonic Boundary-Layer Stabilization by Porous Surfaces", 《AIAA JOURNAL》 * |
| 涂国华 等: "多孔表面抑制第二模态失稳的最优开孔率和孔半径分析", 《空气动力学学报》 * |
| 赵瑞 等: "声学超表面抑制Mack第2模态机理与优化设计", 《气体物理》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651897A (zh) * | 2020-06-18 | 2020-09-11 | 北京理工大学 | 一种能够抑制高超声速边界层转捩的大尺寸缝隙微结构 |
| CN111651897B (zh) * | 2020-06-18 | 2022-10-11 | 北京理工大学 | 一种能够抑制高超声速边界层转捩的大尺寸缝隙微结构 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chuvakhov et al. | Spontaneous radiation of sound by instability of a highly cooled hypersonic boundary layer | |
| Morris et al. | Prediction of broadband shock-associated noise using Reynolds-averaged Navier-Stokes computational fluid dynamics | |
| CN110135062A (zh) | 圆形多孔热防护材料的优化设计方法 | |
| Chen et al. | CFD numerical simulation of exhaust muffler | |
| West et al. | Jet noise prediction using a permeable FW-H solver | |
| CN110135060A (zh) | 方形多孔热防护材料的优化设计方法 | |
| Chakrabarti et al. | The dynamics of azimuthal modes in rectangular jets | |
| Liu et al. | Structured porous blunt trailing edge with uniform and non-uniform parameters for vortex shedding noise reduction | |
| He et al. | Numerical investigation on interference and absorption of electromagnetic waves in the plasma-covered cavity using FDTD method | |
| Fedorov et al. | Instability of hypersonic boundary layer on a wall with resonating micro-cavities | |
| Carton de Wiart et al. | A discontinuous Galerkin method for implicit LES of moderate Reynolds number flows | |
| Crawley et al. | Control of a supersonic rectangular jet using plasma actuators | |
| Knisely et al. | The supersonic mode and the role of wall temperature in hypersonic boundary layers with thermochemical nonequilibrium effects | |
| Ewert et al. | A CAA based approach to tone haystacking | |
| Wong et al. | A parabolised stability equation based broadband shock-associated noise model | |
| Kalyan et al. | Towards a suitable turbulence length and temporal scale model for broadband shock associated noise | |
| Huet et al. | Simulation of flow control with microjets for subsonic jet noise reduction | |
| Brès et al. | Alternate designs of ultrasonic absorptive coatings for hypersonic boundary layer control | |
| Abbasi | Effects of blowing location on aeroacoustics of the flow over a circular cylinder | |
| Zhang et al. | Flow characteristics of an air chamber in a rail tunnel | |
| Gang et al. | Numerical study on noise reduction of wind turbine blade vortex generator | |
| Özdemir et al. | An Analysis of Geometric Parameters‟ Effects on Flow Characteristic of a Reactive Muffler | |
| PHILPOT | The role of rotor blade blockage in the propagation of fan noise interaction tones | |
| Huet | On the use of plasma synthetic jets for the control of jet flow and noise | |
| Lakshmi Narasimha Prasad et al. | Flow-field and Acoustics of Over-expanded Rectangular Jets Subjected to LAFPA Based Control |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhao Rui Inventor after: Zhang Xinxin Inventor after: Wei Haogong Inventor after: Wen Zhiyong Inventor before: Zhao Rui Inventor before: Zhang Xinxin |
|
| CB03 | Change of inventor or designer information | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190816 |
|
| WD01 | Invention patent application deemed withdrawn after publication |