基于等離子流動控制的翼型減阻技術

王元元 周鑄 陳明環

摘要：為研究基于等離子流動控制的減阻技術，基于Langtry-Menter轉捩模型提出邊界層轉捩數值模擬技術。該技術可有效結合轉捩模型與湍流模型，用標準模型驗證其精確性，為采用等離子流動控制抑制邊界層分離和轉捩研究提供數值模擬平臺。采用基于現象學模型的等離子流動控制數值模擬技術，對流動分離以及邊界層轉捩抑制進行數值模擬，為基于等離子流動控制的翼型減阻技術提供參考。

關鍵詞：等離子流動控制;流動分離;邊界層;轉捩;現象學模型

中圖分類號：TP391.99;V211.412

文獻標志碼：B

文章編號：1006-0871（2021）03-0067-06

DOI：10.13340/j.cae.2021.03.012

Abstract：To?study?the?drag?reduction?technology?based?on?plasma?flow?control，?the?numerical?simulation?technology?of?boundary?layer?transition?is?proposed?based?on?Langtry-Menter?transition?model.?The?technique?can?effectively?combine?the?transition?model?and?turbulence?model，?and?its?accuracy?is?verified?by?standard?model.?It?provides?a?numerical?simulation?platform?for?the?suppression?of?boundary?layer?separation?and?transition?by?plasma?flow?control.?The?numerical?simulation?technology?of?plasma?flow?control?based?on?phenomenological?model?is?adopted?to?numerically?simulate?the?flow?separation?and?boundary?layer?transition?suppression.?It?provides?a?reference?for?airfoil?drag?reduction?technology?based?on?plasma?flow?control.

Key?words：plasma?flow?control;flow?separation;boundary?layer;transition;phenomenological?model

0?引?言

基于等離子流動控制的飛行器氣動減阻技術是一種新概念減阻技術，典型的等離子激勵控制方式包括介質阻擋放電、電弧放電、激光離子等，其中介質阻擋放電和電弧放電等對高速流場中激波流動控制機理方面的研究較多[1-5]，可以驗證這些方法對某一范圍空氣流動中的分離、升阻特性和噪聲控制的有效性，揭示動量傳遞以及熱傳導效應機理，尤其在抑制分離方面取得豐碩的成果。但是，在等離子激勵位置對邊界層轉捩抑制效果方面的系統性研究仍然缺乏。

本文面向分離、轉捩等問題和等離子控制技術，建立以N-S為主控方程的計算流體力學模型，對高壓納秒級脈沖介質阻擋放電等離子流動的控制進行數值模擬，研究大型客機的超臨界翼型低速流動減阻控制，驗證該方法對分離抑制的有效性，并初步揭示等離子激勵位置對轉捩位置的抑制規律。

1?邊界層轉捩數值模擬

當M=0.1、Re=1.0×106時，S809翼型自由轉捩計算的升力曲線、阻力特性曲線和力矩曲線數值模擬結果與風洞試驗數據的對比分別見圖1～3，上下表面轉捩點預測與風洞試驗的結果對比見圖4。轉捩位置預測較為準確，可為本文開展邊界層轉捩抑制研究提供可靠的數值依據。

2?納秒級脈沖介質阻擋放電等離子體激勵模型

第一定律模型能夠準確模擬真實空氣中的放電反應，通過求解電荷連續方程、動量方程和電勢方程可模擬等離子體微尺度變化過程。然而，空氣為多組分氣體，放電反應非常復雜，電荷連續方程和動量方程可以多到幾十個，求解難度和計算量也相當大。

對基于現象學的模型進行改進完善，可以很好地描述等離子激勵器的工作原理，可明顯簡化與其他學科耦合數值模擬的復雜程度。

研究對象選擇納秒級高壓脈沖介質阻擋放電等離子流場。對于納秒級脈沖放電，其放電周期處于10-8～10-9?s的量級，與電荷分布周期基本處于一個量級，流動周期約比等離子作用周期慢2個量級，因此可以認為電荷重布是即時連續的，建立激勵器模型可以采用準定常假設。

通過求解麥克斯韋方程模擬等離子體激勵器工作過程的計算量較為龐大，基于電荷重布時間、放電周期以及相應時間之間在量級上的差異，可以推導電勢拉普拉斯方程和電荷誘導電勢泊松方程[9]，即

為進一步進行電場力分布求解，給出電勢和電荷分布邊界條件，分別見圖5和6。

對上述方程采用LUSGS[10]隱式時間進行求解，空間離散格式采用Roe格式[11]，湍流模型采用SST?k-ω模型[12]。

3?基于等離子流動控制的分離抑制

以NACA0015翼型失速流動控制為研究對象，來流速度為34?m/s，迎角15°，將等離子激勵裝置置于上表面1%弦長位置，等離子激勵下的電場強度分布見圖7，控制前的流場分布和渦量場分布分別見圖8和9，等離子激勵器開啟后的流場分布和渦量場分布分別見圖10和11。由此可以看出，等離子激勵作用下，流場分離區域大小得到充分控制，將分離點的位置從前緣推遲到70%弦長位置，等渦量影響區大幅減小，分離控制效果明顯。

4?基于等離子流動控制的層流減阻

以自主設計的超臨界翼型邊界層轉捩流動控制為研究對象，來流速度為20?m/s，迎角為2°，電壓為12?kV，將等離子激勵產生的徹體力引入N-S方程的動量方程和能量方程，激勵器置于不同位置（10%、20%、30%、40%、50%弦長）對邊界層轉捩的不同控制效果（湍動能、摩阻等影響）見圖12～21。由此可以看出，原始翼型的轉捩點大致位于60%弦長處，隨著等離子激勵器位置從10%向后移動到50%弦長位置時，轉捩抑制效果呈現出先變好后變差的變化過程，即當激勵器距離轉捩位置前方30%弦長時（x/c=0.3），控制效果最好，推遲轉捩位置近10%。

5?結?論

（1）基于轉捩模型的邊界層轉捩數值模擬技術，以S809翼型低速繞流的邊界層轉捩進行模擬，通過與風洞試驗數據對比驗證該方法的精度。

（2）基于現象學模型的等離子流動控制數值模擬技術，對NACA0015分離流場進行等離子控制模擬，結果將分離點位置后移70%，控制效果明顯。

（3）研究等離子流動控制技術對邊界層轉捩推遲的影響，以自主設計的超臨界翼型為研究對象，初步揭示低速狀態下激勵器位置對轉捩位置的影響，即激勵器距離轉捩位于前方合適位置時（太近或者太遠都不合適）轉捩抑制效果最好。針對文中案例的翼型，當激勵器距離轉捩位置前方30%弦長時控制效果最好，推遲轉捩位置近10%。

參考文獻：

[1]?POST?M?L，?CORKE?T?C.?Separation?control?on?high?angle?of?attack?airfoil?using?plasma?actuators[J].?AIAA?Journal，?2004，?42（11）：2117-2184.?DOI：10.2514/1.2929.

[2]?MADHAN?R，?RHEL?J，?ROTH?J?R.?Aerodynamic?flow?control?using?paraelectric?and?peristatic?acceleration?of?a?one?atmosphere?uniform?glow?discharge?plasma[C/OL].?（2004-12-31）[2021-08-10].?https：//www.researchgate.net/publication/261415338_Aerodynamic_flow_control_using_paraelectric_and_peristaltic_acceleration_of_a_One_Atmosphere_Uniform_Glow_Discharge_Plasma.?DOI：10.1109/PLASMA.2004.1339868.

[3]?ENLOE?C?L，?MCLAUGHLIN?T?E，?VANDYKEN?R?D，?et?al.?Mechanisms?and?responses?of?a?single?dielectric?barrier?plasma?actuator：Plasma?morphology[J].?AIAA?Journal，?2004，?42（3）：589-594.?DOI：10.2514/1.2305.

[4]?POST?M?L.?Plasma?actuators?for?separation?control?on?stationary?and?oscillating?airfoils[D].?South?Bend：University?of?Notre?Dame，?2004.

[5]?ENLOE?C?L，?MCLAUGHLIN?T?E，?VANDYKEN?R?D，?et?al.?Mechanisms?and?responses?of?a?single?dielectric?barrier?plasma?actuator：Geometric?effects[J].?AIAA?Journal，?2004，?42（3）：595-604.?DOI：10.2514/1.3884.

[6]?MENTER?F?R，?LANGTRY?R?B，?LIKKI?S?R，?et?al.?A?correlation?based?transition?model?using?local?variables?-?part?1：Model?formulation[J].?Journal?of?Turbomachinery，?2006，?128（3）：413-422.?DOI：10.1115/1.2184352.

[7]?LANGTRY?R?B，?MENTER?F?R，?LIKKI?S?R，?et?al.?A?correlation-based?transition?model?using?local?variables：Part?Ⅱ-?Test?cases?and?industrial?applications[C]//?Proceedings?of?ASME?Turbo?Expo：Power?for?Land，?Sea，?and?Air.?Vienna：International?Gas?Turbine?Institute，?2004.?DOI：10.1115/GT2004-53454.

[8]?LANGTRY?R?B，?MENTER?F?R.?Correlation-based?transition?modeling?for?unstructured?parallelized?computational?fluid?dynamics?codes[J].?AIAA?Journal，?2009，?47（12）：2894-2906.?DOI：10.2514/1.42362.

[9]?ORLOV?D，?COEKE?T.?Numerical?simulation?of?aerodynamic?plasma?actuator?effects[C]//?Proceeding?of?43rd?AIAA?Aerospace?Sciences?Meeting?and?Exhibit.?Reno：AIAA，?2005.?DOI：10.2514/6.2005-1083.

[10]?YOON?S，?JAMESON?A.?Lower-upper?Symmetric-Gauss-Seidel?method?for?Euler?and?Navier-Stokes?equations[J].?AIAA?Journal，?1988，?26（9）：1025-1026.?DOI：10.2514/3.10007.

[11]?ROE?P?L.?Approximate?Riemann?solver，?parameter?vectors，?and?difference?schemes[J].?Journal?of?Computational?Physics，?1981，?43（2）：357-372.?DOI：10.1016/0021-9991（81）90128-5.

[12]?MENTER?F?R.?Two-equation?eddy-viscosity?turbulence?models?for?engineering?applications[J].?AIAA?Journal，?1994，?32（8）：269-289.?DOI：10.2514/3.12149.

（編輯?武曉英）