留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种宽角扫描紧耦合阵列天线单元设计

陆娇君 吴鸿超

旷圆圆, 卢艳. 双模瑞利-泰勒不稳定性的预热烧蚀效应研究[J]. 强激光与粒子束, 2022, 34: 082203. doi: 10.11884/HPLPB202234.220133
引用本文: 陆娇君, 吴鸿超. 一种宽角扫描紧耦合阵列天线单元设计[J]. 强激光与粒子束, 2019, 31: 113202. doi: 10.11884/HPLPB201931.190158
Kuang Yuanyuan, Lu Yan. Study on preheating ablative effects of two-mode Rayleigh-Taylor instability[J]. High Power Laser and Particle Beams, 2022, 34: 082203. doi: 10.11884/HPLPB202234.220133
Citation: Lu Jiaojun, Wu Hongchao. Design of a wide-angle scanning tightly coupled dipole array element[J]. High Power Laser and Particle Beams, 2019, 31: 113202. doi: 10.11884/HPLPB201931.190158

一种宽角扫描紧耦合阵列天线单元设计

doi: 10.11884/HPLPB201931.190158
基金项目: 

“十三五”装备预研共用技术项目 41413060404

详细信息
    作者简介:

    陆娇君(1994—), 女,硕士研究生,从事宽带宽角扫描相控阵研究; saraeling@163.com

    通讯作者:

    吴鸿超(1976—), 男,高工,从事超宽带天线、电磁计算研究; babaywind@sina.com

  • 中图分类号: TN822

Design of a wide-angle scanning tightly coupled dipole array element

  • 摘要: 提出一种新型紧耦合偶极子阵列天线单元,并结合等效电路对天线参数进行分析与优化。引入频率选择表面置于天线口径上方,取代传统的电介质板,用以改善宽角扫描时阻抗变换。巴伦采用微带线到共面平行双线过渡结构,实现平衡馈电及阻抗变换。仿真结果表明,该阵列能够实现3倍频(2~6 GHz) E面80°、H面45°的波束扫描,且在扫描范围内有源驻波比均小于3。仿真得到的阵列法向交叉极化隔离度保持在25 dB,由于阻性FSS的损耗,天线增益有所下降。该天线结构简单紧凑,易于加工制作,实现了紧耦合阵列的小型化。
  • 当轻流体加速重流体时,流体交界面处会发生瑞利-泰勒不稳定性(RTI)[1-4]。惯性约束聚变(ICF)[5-8]内爆过程中,靶丸外表面高温高压低密度等离子体加速压缩内层高密度流体时,流体交界面会发生烧蚀瑞利-泰勒不稳定性(ARTI)[9-12],并且由于激光强度分布不均、靶丸表面粗糙程度等因素,ARTI的发展会被进一步放大,从而破坏内爆壳层的对称性,导致靶丸[5]破裂,甚至影响点火热斑的形成。因此,ARTI的研究对聚变靶的内爆压缩及点火烧蚀具有重要意义。通常利用数值模拟方法研究ARTI,首先在无扰动下对一维流场进行稳定性分析,得到一个加速度基本保持不变的平衡流场,然后在烧蚀面处引入物理量的小扰动,追踪流场发展获得进一步的不稳定性分析。以往的单模扰动[13-18]研究表明,ARTI的演化大致可以分为线性阶段与非线性阶段,其增长会受到烧蚀面附近物质的对流[7]、烧蚀面密度梯度[6]和预热[19-21]等因素的影响,其中密度梯度效应对ARTI的线性增长和非线性阶段模耦合发展具有显著影响[22]

    由于引起界面小扰动的因素很多,因此在烧蚀面附近可能引发多模扰动,从而激发许多新的谐波[23-25],这些耦合谐波会进一步影响ARTI的增长。本文研究的主要内容是双模扰动下带有预热的ARTI发展问题,这里的预热是指较低能的电子预热,所以能量是沉积在烧蚀面附近,因此如果预热程度越高,意味着烧蚀面密度梯度越大,烧蚀面宽则烧蚀程度强。由于预热会引起烧蚀面宽度(密度梯度效应)变化,数值模拟中设置强、中等和弱三种不同预热条件,通过一维稳定性分析,首先得到三种不同预热条件下的稳态流场,并在此基础上引入波长不同的两个速度扰动[26]。研究中将数值模拟分为两组:第一组是长波与短波分别在强、中和弱三种预热条件下的耦合;第二组是短波与短波分别在强、中和弱三种预热条件下的耦合。通过数值计算给出高次谐波幅值随时间演化的曲线和ARTI演化一段时间后的气泡尖钉结构图,比较分析两种不同耦合方式下预热强度对谐波演化的影响可以发现,随着预热强度的增大,谐波增长趋势变缓。

    我们采用的物理模型为200μm厚的CH平面靶,初始密度为1g/cm3,激光功率线性上升至4 ns后达到峰值强度1014W/cm3然后维持峰值不变,临界密度为0.03g/cm3。将等离子体视为理想气体,考虑电子热传导,守恒型控制方程如下[27]

    ρt+ρv=0
    (1)
    (ρv)t+(ρvv)+p=0
    (2)
    t(12ρvv+ε)+[(12ρvv+ε+p)v]=(κT)+S
    (3)

    式中:ρvTε=cvρT分别为密度、速度、温度和内能,其中cv=Γ/(γh1)是定容比热,γh=5/3是CH材料中气体绝热指数;κ为电子热传导系数;S为激光能量。

    对于电子热传导,我们采用的模型为

    κ=κSHf(T)h
    (4)

    其中κSH=kST5/2为经典电子热传导系数,h取1,f(T)=c/T3/2+b/T+1为预热函数[19]。数值模拟中可以通过改变bc的值来改变预热程度。我们选取强预热(strong preheating,SP)、中等预热(moderate preheating,MP)和弱预热(weak preheating,WP)所对应的bc的值,如表1所示[28-30]

    表  1  三种预热条件的参数设置
    Table  1.  Parameter settings of strong, medium and weak preheating conditions
    casebc
    strong preheating (SP)8.61.6
    moderate preheating (MP)20.4
    weak preheating (WP)0.860.24
    下载: 导出CSV 
    | 显示表格

    一维计算到8 ns左右,烧蚀面附近流场加速度基本保持不变,达到稳态。8 ns时三种预热情况下的一维平衡流场如图1所示。以最小密度梯度定标长度Lm=min[|ρ(dρ/dx)1|]来表征烧蚀面的宽度,SP,MP和WP三种情况下的一维平衡流场对应的Lm分别为1.830.410.23μm,这表明三种模型中,预热程度越大,烧蚀面越宽。

    图  1  不同预热情况下的一维平衡流场的密度、速度和温度分布
    Figure  1.  Density, velocity and temperature distribution of one-dimensional equilibrium flow field under different preheating conditions

    设一维平衡流场以加速度g作匀加速运动,在加速运动的非惯性系坐标系下,二维形式的控制方程可以写为

    ρt+ρux+ρvy=0
    (5)
    ρut+(ρu2+p)x+ρuvx=ρg
    (6)
    ρvt+ρuvx+(ρv2+p)y=0
    (7)
    ρεt+(ρε+p)ux+(ρε+p)vy=xκTx+y(κTy)+ρug
    (8)

    式中:ux方向速度;vy方向速度。该物理模型中烧蚀面距离激光吸收位置(临界密度附近)较远,因此可以忽略激光能量项S,能量靠电子热传导传输。引入一个初始速度扰动

    u=iuicos(kiy+ϕi)exp(ki|xx0|)
    (9)
    v=ivisin(kiy+ϕi)exp(ki|xx0|)
    (10)

    其中ki=2π/λ表示波数,ϕi表示相位,本文不考虑相位变化,即ϕi=0。通过数值计算得到各个波长下对应的线性扰动增长率如图2所示,WP情况下线性增长率最大,MP次之,SP最小。WP,MP和SP情况下对应的截止波长分别为λWP=4μmλMP=4.2μmλSP=5μm,这表明随着预热增强,截止波长在变大。三种预热情况下,扰动波长λ=18μm左右具有最大线性增长率。

    图  2  强、中和弱预热情况下的线性增长率
    Figure  2.  Linear growth rate under strong, medium and weak preheating

    接下来为了研究预热烧蚀效应对不同模耦合发展的影响,我们计算两组对比数值模拟。第一组扰动选在截止波长附近λ1=10μm的短波与远离截止波长λ2=40μm的长波,第二组在截止波长附近λ1=10μmλ2=12μm的两个短波。不同波长初始幅值相同,为满足扰动周期性特点,y方向计算域为120μm

    随着不稳定性的增长进入弱非线性阶段,两个初始扰动模的相互作用开始变得明显,将产生和频谐波与差频谐波及介于两者之间的其他谐波[31]

    长波与短波耦合所得到谐波幅值η随时间的演化如图3所示。图3(a)为SP情况,图3(b)为MP情况,图3(c)为WP情况。比较发现,在三种预热条件下都产生了波长为40μm20μm的幅值涨幅较大的长波模,以及波长为1310μm的幅值涨幅较小的短波模。在5 ns左右,SP情况下长波模幅值迅速增长,短波模幅值增长被抑制,随后长波模一直占据上风,到7.5 ns左右,40μm20μm的长波模幅值达到第一个极值点30μm左右,随后开始以震荡形式衰减,短波模幅值同时震荡衰减。在MP与WP情况下,5 ns左右长波模迅速增长,短波模幅值增长被抑制,长波模一直占据上风,到7.5 ns左右40μm长波模幅值达到30μm左右,而20μm的长波模幅值达到了70μm左右,随后长波模与短波模开始迅速衰减。这表明长波与短波耦合情况,在三种不同预热程度下,长波模在弱非线性阶段占主导地位,而短波模受到抑制。在预热程度相对较少的MP与WP情况下,20μm的长波模幅值达到的峰值要高很多,这表明预热可以抑制谐波的增长。

    图  3  不同预热情况下λ1=10μmλ2=40μm双模扰动所激发的谐波的密度幅值随时间的演化
    Figure  3.  Temporal evolution of the density amplitude of harmonics excited by the coupling of two modes (short-and long-wavelength λ1=10μm and λ2=40μm) in different preheating cases

    我们给出了波长为λ1=10μmλ2=40μm的扰动下t=7.5 ns时的密度等值线,如图4所示。图4(a)为SP情况,图4(b)为MP情况,图4(c)为WP情况。可以看到此时气泡展现出两种不同的波长的结构,其波长大致为40μm20μm,这是因为模耦合产生的和频与差频谐波(8μm13μm,前者增长不明显,文中没有给出其增长曲线)都为短波,耦合谐波的增长不大,演化中主要增长的模式为基模中的长波模(40μm)以及其二次谐波(20μm)。随着预热程度减小,这两种结构更加明显。这进一步说明在长波与短波的耦合情况下,长波结构占主导地位,且预热程度越高,气泡尖钉增长越慢。

    图  4  不同预热情况下λ1=10μmλ2=40μm双模扰动下在t=7.5 ns时的密度等值线
    Figure  4.  Density contours for λ1=10μm and λ2=40μm mode coupling under different preheating at t=7.5 ns

    短波与短波耦合所得到谐波幅值随时间的演化如图5所示。图5(a)为SP情况,图5(b)为MP情况,图5(c)为WP情况。在三种预热条件下都产生了波长为60302015μm的长波模,以及波长为12μm10μm的幅值涨幅较小的短波模。在SP情况下,5.5 ns左右波长为60μm30μm的长波模的幅值迅速增长,短波模的幅值出现震荡,长波模一直占据上风,随后在8 ns左右,波长为60μm的长波模达到峰值,波长为30μm的长波模继续增长。在MP与WP情况下,5 ns左右波长为60μm30μm的长波模迅速增长,长波模同样占据上风,而短波的幅值在整个过程中都处于小幅值震荡状态。MP和WP情况下,长波模幅值达到的最大幅值的时间比SP情况要早5 ns左右。这表明在短波与短波耦合中,会激发原本基模扰动中没有的长波模(该模态的波长满足两个短波波长的最小公倍数值),并且在弱非线性阶段长波模占据主导,在预热较多的SP情况下,弱非线性阶段谐波幅值达到峰值的时间要更长。

    图  5  不同预热情况下λ1=10μmλ2=12μm双模扰动所激发的谐波的密度幅值随时间的演化
    Figure  5.  Temporal evolution of the density amplitude of harmonics excited by the coupling of two short-wavelength modes (λ1=10μm and λ2=12μm) in different preheating cases

    波长为λ1=10μmλ2=12μm的扰动下,t =7.5 ns时的密度等值线如图6所示。图6(a)为SP情况,图6(b)为MP情况,图6(c)为WP情况。此时气泡展现出多种不同的波长的结构,其中最为明显的是大约60μm30μm的长波结构。在MP和WP情况下,60μm的伞状尖钉结构更加明显。这也进一步说明了在短波与短波的耦合情况下,依然是长波结构占主导地位,且预热程度越小,气泡尖钉增长越快。

    图  6  不同预热情况下λ1=10μmλ2=12μm双模扰动下在t=7.5 ns时的密度等值线
    Figure  6.  Density contours for λ1=10μm and λ2=12μm mode coupling under different preheating at t=7.5 ns

    本文采用数值求解流体方程的方法研究了不同预热条件下的双模扰动烧蚀-瑞利泰勒不稳定性的演化规律。通过一维稳态分析得到了强、中等和弱预热条件下的平衡流场,并在此基础上引入两个速度扰动,获得了双模耦合谐波幅值随时间的演化。在三种预热条件下,不论是长波与短波耦合还是短波与短波耦合,都会激发出许多新的高次谐波,并且谐波中的长波总是占主导。在我们设置的预热程度范围内,预热程度越大,长波结构增长越缓慢,这表明预热会抑制双模耦合高次谐波的增长,对烧蚀瑞利-泰勒不稳定性具有一定的致稳作用。

  • 图  1  理想馈源馈电的TCDA模型及等效电路

    Figure  1.  TCDA model and equivalent circuit for ideal feed feeding

    图  2  (a) TCDA结构(b)宽带巴伦结构(c)宽带巴伦S参数仿真结果

    Figure  2.  (a) TCDA structure, (b) broadband balun structure and (c) broadband balun S-parameter simulation results

    图  3  FSS等效电路

    Figure  3.  Equivalent circuit of FSS

    图  4  电介质层与FSS仿真结果

    Figure  4.  Dielectric layer and FSS simulation results

    图  5  TCDA加载FSS时有源驻波系数特性曲线

    Figure  5.  Active VSWR characteristic curves loaded with FSS

    图  6  TCDA轴向可实现增益

    Figure  6.  Broadside realized gain of the TCDA

    图  7  20×∞有限阵列

    Figure  7.  20×∞ finite array

    图  8  20×∞阵列归一化方向图,E面:0°/30°/45°/60°/70°/80°;H面:0°/30°/45°

    Figure  8.  20×∞ finite array normalized pattern E-plane: 0°/30°/45°/60°/70°/80°; H-plane: 0°/30°/45°

    表  1  优化参数值

    Table  1.   Optimized parameter value  (mm)

    dE dH w l W2 W3 W4 W5 l2 l3 l4 Wb
    14.3 14.3 3.4 7.7 0.55 1.25 3 3 5 18 3.2 8
    下载: 导出CSV
  • [1] Caminita F, Costanzo S, Massa G D, et al. Reduction of patch antenna coupling by using a compact EBG formed by shorted strips with interlocked branch-stubs[J]. IEEE Antennas and Wireless Propagation Letters, 2009, 8: 811-814. doi: 10.1109/LAWP.2009.2021589
    [2] 唐明春, 肖绍球, 高山山, 等. 新型电谐振人工异向介质抑制阵列天线单元间互耦[J]. 物理学报, 2010, 59(3): 1851-1856. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201003067.htm

    Tang Mingchun, Xiao Shaoqiu, Gao Shanshan, et al. New type of electrical resonance artificial anisotropic medium suppresses mutual coupling between array antenna elements. Acta Physica Sinica, 2010, 59(3): 1851-1856 https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201003067.htm
    [3] Munk B, Pryor J, Durham T, et al. A low-profile broadband phased array antenna[C]//Antennas & Propagation Society International Symposium. 2003: 448-451.
    [4] Cedric W L, Dimitris P, Asimina K, et al. Body-worn 67: 1 bandwidth antenna using 3 overlapping dipole elements[C]//11th European Conference on Antennas and Propagation. 2017: 1557-1558.
    [5] Zhong J, Alwan E A, Volakis J L. Ultra-wideband dual-linear polarized phased array with 60° scanning for simultaneous transmit and receive systems[C]//International Workshop on Antenna Technology: Small Antennas, Innovative Structures, and Applications. 2017: 10-141.
    [6] Xiao Shiwei, Yang Shiwen, Chen Yikai, et al. An ultra-wideband tightly coupled dipole array co-designed with low scattering characteristics[J]. IEEE Transactions on Antennas and Propagation, 2019, 67(1): 676-680. doi: 10.1109/TAP.2018.2880038
    [7] Doane J. Wideband low-profile antenna arrays: fundamental limits and practical implementations[J]. Dissertations & Theses-Gradworks, 2013, 76(5): 266.
    [8] Li Xue, Gu Changqing, Han Guodong, et al. A UWB wide-scan tightly coupled dipole array[C]//IEEE International Workshop on Electromagnetics: Applications and Student Innovation Competition. 2016: 1-3.
    [9] Yetisir E, Ghalichechian N, Volakis J L. Ultrawideband array with 70° scanning using FSS superstrate[J]. IEEE Transactions on Antennas and Propagation, 2016, 64(10): 4256-4265. doi: 10.1109/TAP.2016.2594817
    [10] Kasemodel J A, Chen C C, Volakis J L. Wideband planar array with integrated feed and matching network for wide-angle scanning[J]. IEEE Transactions on Antennas and Propagation, 2013, 61(9): 4528-4537. doi: 10.1109/TAP.2013.2266090
    [11] Zhou Wenyang, Chen Yikai, Yang Shiwen. Dual-polarized tightly coupled dipole array for UHF-X band satellite applications[J]. IEEE Antennas and Wireless Propagation Letters, 2016, 18(3): 467-471.
    [12] Novak M H, Volakis J L. Ultrawideband antennas for multiband satellite communications at UHF-Ku frequencies[J]. IEEE Transactions on Antennas and Propagation, 2015, 63(4): 1334-1341. doi: 10.1109/TAP.2015.2390616
    [13] 陈常青, 周志鹏, 张金平, 等. 一种基于紧耦合结构的超宽带天线阵列设计[J]. 现代雷达, 2018, 40(2): 73-76. https://www.cnki.com.cn/Article/CJFDTOTAL-XDLD201802015.htm

    Chen Changqing, Zhou Zhipeng, Zhang Jinping, et al. Design of ultra-wideband antenna array based on tightly coupled structure. Modern Radar, 2018, 40(2): 73-76 https://www.cnki.com.cn/Article/CJFDTOTAL-XDLD201802015.htm
    [14] Papantonis D, Volakis J. Dual polarized tightly coupled array with substrate loading[J]. IEEE Antennas & Wireless Propagation Letters, 2016, 15: 325-328.
    [15] Munk B A. Broadband wire arrays[M]. Finite Antenna Arrays and FSS, 2009.
    [16] Magill E, Wheeler H. Wide-angle impedance matching of a planar array antenna by a dielectric sheet[J]. IEEE Transactions on Antennas and Propagation, 1966, 14(1): 49-53. doi: 10.1109/TAP.1966.1138622
    [17] Mailloux R J. Phased array antenna handbook[M]. 2nd Ed. Beijing: Publishing House of Electronics Industry, 2008.
  • 期刊类型引用(2)

    1. 徐伟东,叶文怡,王炅,成文凭,赵伟康,徐蓉,严萍. 增强型电磁轨道发射技术现状及发展趋势. 高电压技术. 2023(02): 871-884 . 百度学术
    2. 杜翔宇,刘少伟,关娇,时建明. 基于电接触特性的电枢臂形状优化设计. 弹道学报. 2022(03): 103-110 . 百度学术

    其他类型引用(0)

  • 加载中
图(8) / 表(1)
计量
  • 文章访问数:  1967
  • HTML全文浏览量:  638
  • PDF下载量:  107
  • 被引次数: 2
出版历程
  • 收稿日期:  2019-05-13
  • 修回日期:  2019-08-12
  • 刊出日期:  2019-11-15

目录

/

返回文章
返回