應(yīng)用等離子體激勵(lì)器抑制串列雙圓柱噪聲的數(shù)值仿真研究

施傲 曹琦 燕群 劉興強(qiáng) 趙越

摘 要：針對(duì)等離子體激勵(lì)下的串列雙圓柱繞流噪聲抑制問(wèn)題，通過(guò)將等離子體體積力模型、脫落渦模擬、聲比擬理論等技術(shù)相結(jié)合的數(shù)值模擬方法，研究不同來(lái)流速度下等離子體激勵(lì)器安裝位置對(duì)雙圓柱分離流形態(tài)控制與遠(yuǎn)場(chǎng)噪聲抑制效果的影響。結(jié)果表明，當(dāng)所施加的等離子體激勵(lì)位于圓柱流動(dòng)分離點(diǎn)附近時(shí)，控制措施可有效減小分離渦尺度和湍流強(qiáng)度，并顯著降低遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)的總聲壓級(jí)。隨著來(lái)流速度增大，等離子體激勵(lì)器的降噪效果增強(qiáng)，同時(shí)最優(yōu)安裝位置前移。當(dāng)來(lái)流速度達(dá)到55m/s時(shí)獲得最優(yōu)降噪效果，其遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)聲壓級(jí)頻譜峰值和總聲壓級(jí)分別降低11.5dB和8.3dB。而隨著來(lái)流速度的進(jìn)一步增大，等離子體激勵(lì)器的降噪效果逐漸減弱。所得結(jié)果對(duì)于等離子體流動(dòng)控制抑制串列圓柱噪聲的實(shí)際應(yīng)用有一定指導(dǎo)意義。

關(guān)鍵詞：串列圓柱； 等離子體； 主動(dòng)流動(dòng)控制； 數(shù)值模擬； 氣動(dòng)噪聲

中圖分類(lèi)號(hào)：V211.3 文獻(xiàn)標(biāo)識(shí)碼：A DOI：10.19452/j.issn1007-5453.2023.04.003

隨著近年來(lái)國(guó)內(nèi)外民用航空運(yùn)輸業(yè)的高速發(fā)展，其所帶來(lái)的噪聲污染越發(fā)嚴(yán)重[1]。由起落架機(jī)構(gòu)復(fù)雜外形導(dǎo)致的氣動(dòng)噪聲作為大型客機(jī)起降階段的主要機(jī)體噪聲源之一，引起了工業(yè)和學(xué)術(shù)界的高度重視，起落架外流中存在大分離流動(dòng)和渦系干擾等復(fù)雜現(xiàn)象，具有較大的流動(dòng)尺度變化和較強(qiáng)的非定常特性。為方便開(kāi)展流動(dòng)和噪聲機(jī)理分析，常使用簡(jiǎn)化的串列雙圓柱模型[2]來(lái)表征起落架基本結(jié)構(gòu)用以研究其流動(dòng)特性?；诖?，已有諸多研究人員采用不同的研究手段，圍繞起落架降噪設(shè)計(jì)這一目標(biāo)開(kāi)展了主、被動(dòng)流動(dòng)控制措施的噪聲抑制效果相關(guān)研究工作。

在被動(dòng)控制降噪方面，Chen Xu等[3]通過(guò)對(duì)圓柱繞流下游區(qū)域添加多孔介質(zhì)的數(shù)值模擬，發(fā)現(xiàn)下游流動(dòng)由湍流轉(zhuǎn)變?yōu)閷恿?，?duì)控制其噪聲水平極為有效。F.V.Hutcheson等[4]通過(guò)在多圓柱表面添加O形環(huán)狀凸起結(jié)構(gòu)，實(shí)現(xiàn)了顯著降噪效果。在此基礎(chǔ)上，王大慶等[5]更細(xì)致地研究了多種不同參數(shù)的表面凸起結(jié)構(gòu)對(duì)雙圓柱干涉噪聲的影響規(guī)律。相比于被動(dòng)控制方法，主動(dòng)控制措施雖然因結(jié)構(gòu)復(fù)雜而難以在工程中得到廣泛應(yīng)用，但其可根據(jù)流動(dòng)環(huán)境進(jìn)行即時(shí)調(diào)節(jié)的優(yōu)秀性能仍廣泛地獲得了研究人員的青睞。 D.Angland等[6]通過(guò)試驗(yàn)驗(yàn)證了對(duì)圓柱采用吹氣措施可以有效地抑制剪切層厚度增長(zhǎng)，進(jìn)而顯著降低寬頻噪聲。李潮隆等[7]在圓柱繞流問(wèn)題中加入合成雙射流控制措施，發(fā)現(xiàn)該方法可通過(guò)改善虛擬迎流型面并加速分離渦再附，有效地抑制流動(dòng)分離。

在眾多流動(dòng)控制方法中，等離子體技術(shù)通過(guò)在電磁場(chǎng)力的作用下對(duì)流場(chǎng)施加可控?cái)_動(dòng)，具有結(jié)構(gòu)簡(jiǎn)單、控制位置靈活、能耗低、響應(yīng)快等優(yōu)點(diǎn)[8]，逐漸成為了近20年來(lái)流動(dòng)控制研究領(lǐng)域的熱點(diǎn)。J.Roth等[9]通過(guò)在低速風(fēng)洞中對(duì)平板流動(dòng)施加等離子體控制，展示了等離子體激勵(lì)對(duì)流動(dòng)的誘導(dǎo)作用。T.McLaughlin等[10]通過(guò)圓柱繞流控制試驗(yàn)表明，等離子體激勵(lì)器可有效調(diào)控分離渦的脫落頻率和展向相關(guān)性，從而影響分離流動(dòng)形態(tài)。J.Huang等[11]通過(guò)試驗(yàn)分析了等離子體激勵(lì)器對(duì)渦輪葉片分離流動(dòng)的抑制效果，發(fā)現(xiàn)在流動(dòng)分離點(diǎn)附近位置加載激勵(lì)的效果最佳。M.Forte等[12]研究了等離子體激勵(lì)器的數(shù)量、間距與電壓幅值對(duì)平板流動(dòng)的影響，結(jié)果表明，多個(gè)激勵(lì)器串聯(lián)能有效提高流動(dòng)誘導(dǎo)速度。

在等離子體激勵(lì)控制圓柱流動(dòng)方面，F(xiàn).O.Thomas等[13]通過(guò)單介質(zhì)阻擋放電（SDBD）等離子體激勵(lì)器，顯著降低圓柱尾流的湍流和噪聲水平。D.Kim等[14]在此基礎(chǔ)上采用大渦模擬方法和體積力模型進(jìn)行了數(shù)值仿真，結(jié)果和試驗(yàn)吻合良好，驗(yàn)證了該數(shù)值方法的準(zhǔn)確性。A.V.Kozlov等[15]將展向布置等離子體激勵(lì)器和流向渦流發(fā)生器（PSVG）應(yīng)用于雙圓柱繞流問(wèn)題，試驗(yàn)表明，兩種方法在不同雷諾數(shù)下可有效降低下游圓柱表面的壓力脈動(dòng)水平。L.A.Al-Sadawi等[16]研究了介質(zhì)阻擋放電（DBD）等離子體激勵(lì)器對(duì)單圓柱和串聯(lián)圓柱結(jié)構(gòu)氣動(dòng)噪聲的影響，結(jié)果表明上游圓柱的渦脫落是聲頻譜中峰值噪聲的主要來(lái)源，且可通過(guò)同時(shí)對(duì)上、下游圓柱使用等離子體激勵(lì)來(lái)進(jìn)一步提高控制效果。國(guó)內(nèi)的蘇長(zhǎng)兵等[17]通過(guò)基于等離子體控制的單圓柱繞流試驗(yàn)驗(yàn)證了該措施可有效抑制鈍體分離流，且其效果隨雷諾數(shù)的增大而降低。王建明等[18]通過(guò)試驗(yàn)研究了等離子體激勵(lì)參數(shù)對(duì)單圓柱繞流控制效果的影響。結(jié)果表明，激勵(lì)頻率對(duì)控制效果的影響較大，且當(dāng)激勵(lì)位于流動(dòng)分離點(diǎn)附近區(qū)域時(shí)控制效果較優(yōu)。李紅麗等[19]通過(guò)試驗(yàn)驗(yàn)證了等離子體激勵(lì)抑制串列雙圓柱氣動(dòng)噪聲的效果。賴(lài)晨光等[20]通過(guò)數(shù)值仿真研究了來(lái)流速度和等離子體激勵(lì)器安裝位置對(duì)單圓柱流動(dòng)控制效果的影響。結(jié)果表明，其控制效果在激勵(lì)安裝于下游分離區(qū)表面時(shí)更好，且隨來(lái)流速度增大而減弱。

目前已開(kāi)展的圓柱繞流等離子體流動(dòng)控制相關(guān)研究主要集中于流動(dòng)機(jī)理與控制效果影響分析，很少涉及等離子體對(duì)遠(yuǎn)場(chǎng)噪聲的影響趨勢(shì)分析。國(guó)內(nèi)鮮有利用數(shù)值模擬方法開(kāi)展等離子體激勵(lì)對(duì)串列雙圓柱繞流問(wèn)題的噪聲控制效果研究?；诖?，本文以第一屆機(jī)體噪聲計(jì)算研討會(huì)標(biāo)模算例（BANC-1）中雙圓柱模型為研究對(duì)象[2]，通過(guò)數(shù)值模擬方法探究來(lái)流速度與安裝位置對(duì)等離子體激勵(lì)的流動(dòng)控制效果和噪聲抑制效果的影響趨勢(shì)。

1 等離子體激勵(lì)模型

本文使用目前應(yīng)用較為廣泛的介質(zhì)阻擋放電等離子體激勵(lì)器，其電極之間有起鎮(zhèn)流作用的電介質(zhì)層，可以使放電趨于穩(wěn)定。該激勵(lì)器由交流電源、絕緣電介質(zhì)與電極組成。電極分為上下兩電極，上層電極安裝于絕緣材料表面，下方電極內(nèi)嵌于絕緣介質(zhì)中。如圖1所示。

等離子體激勵(lì)器工作時(shí)，通過(guò)強(qiáng)電場(chǎng)作用力將激勵(lì)器作用區(qū)域內(nèi)的空氣擊穿、電離。被電離的帶電粒子在空間不均勻電場(chǎng)的作用下定向運(yùn)動(dòng)，形成高速的等離子風(fēng)，并通過(guò)引射效應(yīng)誘導(dǎo)周?chē)鷼怏w運(yùn)動(dòng)。

目前常用的等離子體激勵(lì)器建模仿真方法[21]有唯象學(xué)仿真法、集總電路仿真法、物理原理多場(chǎng)耦仿真法和粒子群-蒙特卡羅模擬法。由于串列雙圓柱繞流問(wèn)題中存在復(fù)雜的流動(dòng)分離與渦系運(yùn)動(dòng)，為兼顧計(jì)算資源與時(shí)間消耗以及對(duì)仿真結(jié)果的可接受性，本文采用Shyy提出的體積力模型[22]。該模型不關(guān)注放電的具體過(guò)程和所涉及的化學(xué)反應(yīng)，而是基于試驗(yàn)數(shù)據(jù)與理論研究結(jié)果，將激勵(lì)效果簡(jiǎn)化為施加在一定區(qū)域的體積力。之后作為Navier-Stokes（N-S）方程的源項(xiàng)，通過(guò)求解N-S方程模擬等離子體激勵(lì)控制對(duì)流場(chǎng)流動(dòng)的影響。

該模型將等離子體的激勵(lì)作用影響范圍簡(jiǎn)化為如圖2所示的三角形OAB。該作用區(qū)域內(nèi)電場(chǎng)強(qiáng)度呈線(xiàn)性分布，在原點(diǎn)O處具有最大電場(chǎng)強(qiáng)度E0，在OA、OB和AB邊上的電場(chǎng)強(qiáng)度為截?cái)嚯妶?chǎng)強(qiáng)度Eb。當(dāng)電場(chǎng)強(qiáng)度小于截?cái)鄰?qiáng)度時(shí)，電場(chǎng)力不足以形成等離子層，因而可以在其外部區(qū)域不考慮等離子體作用。

圖3給出了數(shù)值模擬靜止空氣中等離子體激勵(lì)所形成誘導(dǎo)流場(chǎng)的流線(xiàn)分布。流動(dòng)在等離子體作用下向暴露電極末端附近位置聚集，并以射流形式向下游發(fā)展，射流速度可達(dá)約1m/s。其流動(dòng)形式與試驗(yàn)結(jié)果一致。

圖4給出了x = 0.005m站位處的流動(dòng)速度沿壁面法向分布曲線(xiàn)。等離子體加速效果隨壁面距離減小而增大，在距壁面約0.003m處流體速度達(dá)到最大值。隨著壁面距離近一步減小，流體速度因壁面黏性而逐漸減小到零。數(shù)值模擬與試驗(yàn)結(jié)果吻合良好，驗(yàn)證了本文所使用等離子體激勵(lì)數(shù)值模擬方法的準(zhǔn)確性。

2 數(shù)值方法驗(yàn)證

2.1 計(jì)算模型與計(jì)算設(shè)置

本文采用BANC-1中的串列雙圓柱繞流標(biāo)準(zhǔn)模型，圓柱直徑D = 0.05715m，兩圓柱的中心距離L = 3.7D。來(lái)流速度v = 44m/s，馬赫數(shù)Ma = 0.1274，基于圓柱直徑的雷諾數(shù)Re = 1.66×10。對(duì)比試驗(yàn)數(shù)據(jù)來(lái)自美國(guó)國(guó)家航空航天局（NASA）的QFF模型試驗(yàn)[2]。模型設(shè)置和角度定義如圖5所示。

流場(chǎng)計(jì)算采用有限體積N-S方程模擬軟件Fluent。計(jì)算域?yàn)?.5m×6m的矩形，上游圓柱中心與計(jì)算域入口和出口的距離分別為3m和6.5m，與上、下邊界的距離為3m。邊界條件采用速度入口、壓力出口，圓柱壁面為無(wú)滑移的固壁?？臻g離散采用多塊結(jié)構(gòu)網(wǎng)格所轉(zhuǎn)化的非結(jié)構(gòu)網(wǎng)格，如圖6所示。為充分捕捉不同尺度的渦系結(jié)構(gòu)，對(duì)圓柱周?chē)拔擦鲄^(qū)域進(jìn)行局部網(wǎng)格加密處理。近壁處第一層網(wǎng)格厚度為1.4×10-5m。網(wǎng)格總數(shù)在7萬(wàn)個(gè)左右。

定常計(jì)算收斂后，采用SA-DDES（spalart allmarasdelayed detached eddy simulation）模型進(jìn)行非定常流動(dòng)模擬，使用可穿透積分面的FW-H（ffowcs williams-hawkings）方程計(jì)算遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)的聲學(xué)信息。積分面的選擇區(qū)域包括了圓柱壁面以及部分尾跡區(qū)域，如圖7所示。計(jì)算時(shí)間步長(zhǎng)Δt = 2.5×10-4s，采集數(shù)據(jù)時(shí)間為0.175s。

2.2 結(jié)果對(duì)比

將計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，其前后圓柱的表面時(shí)均壓力系數(shù)Cp分布如圖8所示?？梢钥吹?，前圓柱的表面時(shí)均壓力系數(shù)與試驗(yàn)結(jié)果吻合良好。而后圓柱表面時(shí)均壓力系數(shù)的負(fù)壓峰值和角度與試驗(yàn)結(jié)果存在一定的偏移，即數(shù)值模擬得到的下游圓柱邊界層分離較早，尾流分離區(qū)域更大。根據(jù)參考文獻(xiàn)[2]中的研究結(jié)果，多位學(xué)者使用不同的數(shù)值方法均一定程度上出現(xiàn)了類(lèi)似誤差。其原因可能是在圓柱間距較敏感的情況下，前圓柱尾流形態(tài)的微小誤差將導(dǎo)致其對(duì)后圓柱流動(dòng)產(chǎn)生較大影響。圖9給出了前后圓柱表面壓力系數(shù)均方根（RMS）分布對(duì)比。圖中數(shù)值模擬與試驗(yàn)結(jié)果貼合較好。

圖10給出了位于（9.11D，32.49D）處的遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)聲壓級(jí)頻譜特性對(duì)比?？梢园l(fā)現(xiàn)，數(shù)值預(yù)測(cè)結(jié)果的峰值聲壓級(jí)與試驗(yàn)結(jié)果基本一致，誤差僅為1.46%。峰值頻率偏移5Hz，具有較高的準(zhǔn)確性。數(shù)值計(jì)算結(jié)果在高于600Hz的頻段明顯低于試驗(yàn)數(shù)據(jù)。這類(lèi)高頻誤差通常是由離散網(wǎng)格的尺度截?cái)嘈?yīng)所導(dǎo)致的，不影響整體計(jì)算精度。綜上所述，本文所采用的數(shù)值模擬方法具有較好的預(yù)測(cè)精度。

3 等離子體激勵(lì)對(duì)流動(dòng)分離與降噪效果的影響

為了探究來(lái)流速度和等離子體激勵(lì)位置對(duì)雙圓柱繞流的流動(dòng)與降噪效果的影響趨勢(shì)，在研究中保持等離子體激勵(lì)器參數(shù)不變。由于等離子體誘導(dǎo)的速度大小有限，本文研究的來(lái)流速度為20～100m/s。根據(jù)J.Huang等[11]的試驗(yàn)結(jié)果，等離子體激勵(lì)施加在分離點(diǎn)附近時(shí)的抑制分離效果最好。因此，本文的等離子體激勵(lì)加載位置在85°～94°范圍內(nèi)。

3.1 無(wú)控制措施的雙圓柱繞流分析

圖11為35m/s、44m/s、55m/s和85m/s來(lái)流速度的雙圓柱繞流流線(xiàn)圖。可以看出，氣流在受到前圓柱阻擋后沿著圓柱表面向下游發(fā)展，并發(fā)生了明顯的氣流分離現(xiàn)象，流體流經(jīng)下游圓柱在90°附近再次發(fā)生流動(dòng)分離。來(lái)流速度越高，前圓柱尾部的分離區(qū)域越大，但是前后圓柱分離點(diǎn)的位置變化不顯著。

圖12為各來(lái)流速度下的前圓柱的壓力系數(shù)均方根值曲線(xiàn)圖?？梢钥闯觯€(xiàn)的峰值位置（以左峰值為例）逐漸由91°增大到95°，即分離位置隨著來(lái)流速度的增大而向后移動(dòng)。其原因是流速增大時(shí)，相當(dāng)于給氣體增加了一定的動(dòng)能，從而增強(qiáng)了邊界層抵抗逆壓梯度的能力，導(dǎo)致分離點(diǎn)后移。

圖13給出了各來(lái)流速度下（9.11D，32.49D）處監(jiān)測(cè)點(diǎn)的聲壓級(jí)頻譜特性。可以發(fā)現(xiàn)，各來(lái)流速度下雙圓柱遠(yuǎn)場(chǎng)聲壓級(jí)曲線(xiàn)的趨勢(shì)基本一致。隨著來(lái)流速度的增大，峰值噪聲從92.8dB增大到112.8dB，且峰值頻率也從75.7Hz平移到203.8Hz，寬頻噪聲也有約3dB的降低。其原因?yàn)閬?lái)流速度越高，前圓柱繞流流速越大，湍流能量變大，遠(yuǎn)場(chǎng)噪聲強(qiáng)度增大，單位時(shí)間內(nèi)前圓柱脫落渦撞擊后圓柱次數(shù)增加，即渦脫落頻率增加，使得峰值頻率變大。

3.2 激勵(lì)位置對(duì)控制效果的影響分析

圖14為來(lái)流速度為0時(shí)等離子體在90°位置的平均速度云圖。其中，等離子體加載位置是指等離子體作用區(qū)域的左邊界。由圖可發(fā)現(xiàn)在等離子體體積力作用下，流動(dòng)沿壁面向下游方向加速，形成類(lèi)似射流的流動(dòng)形式。值得注意的是，等離子體誘導(dǎo)的流動(dòng)加速效果在體積力作用區(qū)域的中部達(dá)到最強(qiáng)。該處的誘導(dǎo)流速為20m/s，距離等離子體作用區(qū)起始位置約3°夾角的距離。

圖15給出了來(lái)流速度為35m/s時(shí)，無(wú)等離子體激勵(lì)以及等離子體激勵(lì)施加位置為70°、90°和110°的流線(xiàn)圖。對(duì)比發(fā)現(xiàn)，施加等離子體激勵(lì)后，前圓柱分離點(diǎn)向后推遲，流動(dòng)分離得到明顯抑制。其原因可能是等離子體激勵(lì)誘導(dǎo)流體加速，增加近壁面附近流動(dòng)能量，抵抗壓差力與黏性力的能力增強(qiáng)，從而抑制流動(dòng)的分離。由圖發(fā)現(xiàn)，在90°位置施加等離子體激勵(lì)的分離區(qū)更小，抑制流動(dòng)分離的效果更好。其原因分析為當(dāng)施加激勵(lì)位置在流動(dòng)分離點(diǎn)之前時(shí)，激勵(lì)作用區(qū)域流動(dòng)的能量較高，受到等離子體的影響較弱。當(dāng)施加激勵(lì)位置在流動(dòng)分離點(diǎn)之后時(shí)，由于流動(dòng)已經(jīng)分離，等離子體激勵(lì)影響范圍內(nèi)的流體密度變小，導(dǎo)致其能量增加較少，抵抗逆壓梯度的能力變?nèi)?。這也符合之前的研究者得出的將激勵(lì)位置安裝在分離點(diǎn)附近的結(jié)論[9]。

圖16給出了無(wú)等離子體激勵(lì)與等離子體激勵(lì)分別施加在70°、90°和110°位置的遠(yuǎn)場(chǎng)聲壓級(jí)頻譜特性圖。相比于無(wú)等離子體激勵(lì)情況，施加等離子體激勵(lì)后，遠(yuǎn)場(chǎng)聲壓級(jí)曲線(xiàn)圖出現(xiàn)了兩個(gè)峰值。其中，第一個(gè)峰值是由前圓柱渦脫落撞擊后圓柱產(chǎn)生的，第二個(gè)峰值是前圓柱渦脫落產(chǎn)生的[24]。結(jié)合圖17所示的瞬時(shí)渦量圖，可以發(fā)現(xiàn)在等離子體激勵(lì)作用效果下，前圓柱剪切流強(qiáng)度降低，渦脫落撞擊強(qiáng)度降低，從而曲線(xiàn)中第一個(gè)峰值降低。第二個(gè)峰值出現(xiàn)的原因是等離子體激勵(lì)誘導(dǎo)流體加速，增大了前圓柱的卡門(mén)渦街的總能量，湍流能量增大，遠(yuǎn)場(chǎng)聲壓級(jí)增大。

為了探究等離子體施加在分離點(diǎn)附近位置的效果，圖18給出了等離子體施加在85°～94°位置的遠(yuǎn)場(chǎng)監(jiān)測(cè)點(diǎn)聲壓級(jí)頻譜的峰值降低量。從圖中可以看到，在分離點(diǎn)附近施加等離子體激勵(lì)對(duì)遠(yuǎn)場(chǎng)峰值噪聲的降低都有一定的效果。其中，當(dāng)?shù)入x子體安裝在92°位置時(shí)，雙圓柱遠(yuǎn)場(chǎng)峰值噪聲的降低量最大。該安裝角度與流動(dòng)分離點(diǎn)94.7°之間的夾角為2.7°。這與前文分析等離子體在作用區(qū)域起始位置3°夾角的距離后可使流體達(dá)到最大誘導(dǎo)速度的結(jié)果相符。

3.3 來(lái)流速度對(duì)控制效果的影響分析

圖19給出了各個(gè)速度下等離子體激勵(lì)的最優(yōu)角度圖?？梢钥闯鲭S著來(lái)流速度的增加，最佳的等離子體激勵(lì)角度減小。原因分析為隨著來(lái)流速度的增大，在等離子體影響區(qū)域內(nèi)，流體微元在更靠后的區(qū)域達(dá)到了加速的最佳效果。圖20給出了最優(yōu)角度等離子體在各來(lái)流速度的峰值噪聲降低量。由圖可以看出，在低來(lái)流速度時(shí)，等離子體激勵(lì)的降噪效果隨來(lái)流速度增大而增強(qiáng)。來(lái)流速度在55m/s時(shí)，當(dāng)前激勵(lì)參數(shù)下的等離子體的降噪效果達(dá)到最好。隨著來(lái)流速度的進(jìn)一步增大，等離子體的噪聲抑制效果逐漸減弱。其原因可能是高來(lái)流速度時(shí)，流體繞流能量大，等離子體激勵(lì)對(duì)流動(dòng)狀態(tài)的影響較小。在來(lái)流速度較低時(shí)，等離子體激勵(lì)誘導(dǎo)產(chǎn)生的速度與來(lái)流速度相當(dāng)，遠(yuǎn)場(chǎng)的噪聲水平可能會(huì)因?yàn)榈入x子體激勵(lì)誘導(dǎo)速度過(guò)大而增大，使遠(yuǎn)場(chǎng)噪聲的降低量變小。

4 結(jié)論

本文采用基于體積力模型的數(shù)值模擬方法，對(duì)比分析了串列雙圓柱繞流問(wèn)題中來(lái)流速度、安裝位置等關(guān)鍵因素對(duì)等離子體激勵(lì)流動(dòng)控制和噪聲抑制效果的影響。得到以下結(jié)論：

（1） 等離子體體積力模型能夠準(zhǔn)確反映等離子體對(duì)流動(dòng)的誘導(dǎo)加速效果。流動(dòng)速度剖面與試驗(yàn)數(shù)據(jù)吻合度較高。

（2） 對(duì)于低速串列雙圓柱繞流問(wèn)題，在上游圓柱對(duì)稱(chēng)安裝等離子體激勵(lì)能夠有效地抑制流動(dòng)分離，并降低遠(yuǎn)場(chǎng)噪聲水平。在來(lái)流速度為55m/s、等離子體安裝角度為88°時(shí)，遠(yuǎn)場(chǎng)噪聲總聲壓級(jí)降低量最多可達(dá)8.3dB。

（3） 相同來(lái)流速度下，在分離點(diǎn)附近位置施加等離子體激勵(lì)對(duì)流動(dòng)分離具有更明顯的抑制效果，也可獲得更低的遠(yuǎn)場(chǎng)噪聲水平。

（4）等離子體激勵(lì)的流動(dòng)控制和降噪效果在低速條件下隨來(lái)流速度增大而增強(qiáng)，其最優(yōu)安裝角度隨來(lái)流速度增大而減小。當(dāng)來(lái)流速度達(dá)到55m/s并繼續(xù)增大時(shí)，等離子體激勵(lì)對(duì)流動(dòng)分離和噪聲的抑制能力明顯減弱。

受限于計(jì)算耗費(fèi)，本文只圍繞二維串列雙圓柱繞流問(wèn)題開(kāi)展了相關(guān)數(shù)值模擬研究，未能考慮更多不同等離子體激勵(lì)電壓條件。未來(lái)將在本文研究所獲得的規(guī)律的基礎(chǔ)上，面向三維等離子體激勵(lì)實(shí)際應(yīng)用問(wèn)題，開(kāi)展更翔實(shí)的研究。

參考文獻(xiàn)

[1]閆國(guó)華，李成晨. 民機(jī)機(jī)體邊線(xiàn)噪聲預(yù)測(cè)方法研究[J]. 航空科學(xué)技術(shù)，2021，32（3）：6-13. Yan Guohua， Li Chengchen. Research on prediction method of airframe lateral noise of civil aircraft[J]. Aeronautical Science & Technology， 2021，32（3）：6-13. （in Chinese）

[2]Lockard D. Summary of the tandem cylinder solutions from the benchmark problems for airframe noise computations-I workshop[C]//49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition. Orlando， Florida， 2011： 353.

[3]Xu C， Mao Y， Hu Z. Control of cylinder wake flow and noise through a downstream porous treatment[J]. Aeronautical Science & Technology， 2019，88： 233-243.

[4]Hutcheson F V， Brooks T F. Noise radiation from single and multiple rod configurations[J]. International Journal of Aeroacoustics， 2012，11（3-4）： 291-333.

[5]王大慶，仝帆，馮和英，等.節(jié)徑凸起結(jié)構(gòu)對(duì)串列雙圓柱干涉噪聲的影響研究[J]. 航空動(dòng)力學(xué)報(bào)，2022（5）：1-12. Wang Daqing， Tong Fan， Feng Heying， et al. Influence of the pitch diameter protruding structure on the interference noise of tandem double cylinders[J]. Journal of Aerospace Power， 2022（5）：1-12. （in Chinese）

[6]Angland D， Zhang X， Goodyer M. The use of blowing flow control to reduce bluff body interaction noise[C]//16th AIAA/ CEAS Aeroacoustics Conference. Stockholm， Sweden： American Institute of Aeronautics and Astronautics. Stockholm， Sweden， 2010.

[7]李潮隆，夏智勛，羅振兵，等.合成雙射流控制水下圓柱繞流流動(dòng)分離數(shù)值模擬研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào)，2020，38（2）： 254-259+267. Li Chaolong， Xia Zhixun， Luo Zhenbing， et al. Numerical study on flow separation of underwater cylinder with dual synthetic jets actuator[J]. Acta Aerodynamica Sinica， 2020， 38（2）： 254-259+267. （in Chinese）

[8]張野平，侯銀珠，汪發(fā)亮. 飛行器繞流介質(zhì)阻擋放電等離子體流動(dòng)控制技術(shù)綜述[J]. 航空科學(xué)技術(shù)，2016，27（6）：5-10. Zhang Yeping， Hou Yinzhu， Wang Faliang. Brief introduction of DBD plasma flow control in aircraft design[J]. Aeronautical Science & Technology， 2016，27（6）：5-10. （in Chinese）

[9]Roth J， Sherman D， Wilkinson S. Boundary layer flow control with a one atmosphere uniform glow discharge surface plasma[C]//36th AIAA Aerospace Sciences Meeting and Exhibit. Reno， NV， U.S.A， 1998.

[10]McLaughlin T， Munska M， Vaeth J， et al. Plasma-based actuators for cylinder wake vortex control[C]//2nd AIAA Flow Control Conference. Portland， Oregon， 2004.

[11]Huang J， Corke T C， Thomas F O. Unsteady plasma actuators for separation control of low-pressure turbine blades[J]. AIAA Journal， 2006， 44（7）： 1477-1487.

[12]Forte M， Jolibois J， Pons J， et al. Optimization of a dielectric barrier discharge actuator by stationary and non-stationary measurements of the induced flow velocity： Application to airflow control[J]. Experiments in Fluids， 2007， 43（6）： 917-928.

[13]Thomas F O， Kozlov A， Corke T C. Plasma actuators for cylinder flow control and noise reduction[J]. AIAA Journal， 2008， 46（8）： 1921-1931.

[14]Kim D， Wang M. Large-eddy simulation of flow over a circular cylinder with plasma-based control[C]//47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Orlando， Florida， 2009.

[15]Kozlov A V， Thomas F O. Plasma flow control of cylinders in a tandem configuration[J]. AIAA Journal， 2011， 49（10）： 2183-2193.

[16]Al-Sadawi L A， Chong T P. Circular cylinder wake and noise control using DBD plasma actuator[C]//25th AIAA/CEAS Aeroacoustics Conference. Delft， The Netherlands， 2019.

[17]蘇長(zhǎng)兵，宋慧敏，李應(yīng)紅，等.基于等離子體激勵(lì)的圓柱繞流控制實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)流體力學(xué)，2006（4）： 45-48. Su Changbing， Song Huimin， Li Yinghong， et al. Experiments of the flow field structure control around a circular cylinder based on plasma actuation[J]. Journal of Experiments in Fluid Mechanics， 2006（4）： 45-48. （in Chinese）

[18]王建明，江海亮，明曉杰，等.等離子體激勵(lì)參數(shù)對(duì)圓柱繞流影響的風(fēng)洞實(shí)驗(yàn)研究[J]. 科學(xué)技術(shù)與工程，2017，17（35）： 149-155. Wang Jianming， Jiang Hailiang， Ming Xiaojie， et al. Influence of plasma actuation parameters on the flow around a circular cylinder by a wind tunnel experiment[J]. Journal of Research in Science and Engineering， 2017， 17（35）： 149-155. （in Chinese）

[19]李紅麗，劉興強(qiáng)，延浩.等離子體抑制雙柱繞流噪聲實(shí)驗(yàn)研究[J]. 航空工程進(jìn)展，2016，7（3）： 294-300. Li Hongli， Liu Xingqiang， Yan Hao. Experimental investigation of noise control around tandem cylinder based on plasma actuation[J]. Advances in Aeronautical Science and Engineering，2016，7（3）：294-300. （in Chinese）

[20]賴(lài)晨光，付航，凌志偉，等.等離子體激勵(lì)控制圓柱繞流的大渦模擬研究[J]. 重慶理工大學(xué)學(xué)報(bào)（自然科學(xué)），2021，35（8）： 81-89. Lai Chenguang， Fu Hang， Ling Zhiwei， et al. Large eddy simulation of plasma actuation to control flow around a cylinder[J]. Journal of Chongqing University of Technology（Natural Science）， 2021，35（8）： 81-89. （in Chinese）

[21]李應(yīng)紅，梁華，吳云，等.等離子體氣動(dòng)激勵(lì)建模仿真綜述[J].空軍工程大學(xué)學(xué)報(bào)（自然科學(xué)版），2008（5）： 1-5. Li Yinghong， Liang Hua， Wu Yun， et al. Research status quo of modeling and simulation technology of plasma aerodynamic actuation[J]. Journal of Air Force Engineering University（Natural Science Edition）， 2008（5）： 1-5. （in Chinese）

[22]Shyy W， Jayaraman B， Andersson A. Modeling of glow discharge-induced fluid dynamics[J]. Journal of applied physics， American Institute of Physics， 2002， 92（11）： 6434-6443.

[23]Rizzetta D P， Visal M R. Large-eddy simulation of plasmabased turbulent boundary-layer separation control[J]. AIAA Journal， 2010， 48（12）： 2793-2810.

[24]Eltaweel A， Wang M， Kim D， et al. Numerical investigation of tandem-cylinder noise reduction using plasma-based flow control[J]. Journal of Fluid Mechanics， 2014， 756： 422-451.

Numerical Simulation Study on Plasma Noise Control of Cylinders in a Tandem Configuration

Shi Ao， Cao Qi， Yan Qun， Liu Xingqiang， Zhao Yue Aircraft Strength Research Institute， Xi’an 710065，China

Abstract： Focusing on tandem cylinder noise suppression via plasma actuator， the influence of flow velocity and plasma applying positions on flow separation control and far field noise reduction is studied by combining plasma body-force model， detached eddy simulation and acoustic analogy theory. The results show that when plasma actuator is applied near the separation position of front cylinder， the turbulence intensity and far-field total sound pressure level can be effectively reduced. As flow velocity increases， the effect of plasma noise control is enhanced， and the optimal position of plasma actuator moves forward. When the flow velocity comes to 55m/s， a significant noise suppression is observed. The tonal peak of sound pressure level spectra and total sound pressure level decreases by 11.5dB and 8.3dB， respectively. With further increase of flow velocity， the effect of plasma noise control becomes weakened.

Key Words： tandem cylinder； plasma； active flow control； numerical simulation； aerodynamic noise