CN108871595B - 超时间分辨冲击波速度计算方法 - Google Patents
超时间分辨冲击波速度计算方法 Download PDFInfo
- Publication number
- CN108871595B CN108871595B CN201810839877.8A CN201810839877A CN108871595B CN 108871595 B CN108871595 B CN 108871595B CN 201810839877 A CN201810839877 A CN 201810839877A CN 108871595 B CN108871595 B CN 108871595B
- Authority
- CN
- China
- Prior art keywords
- light
- phase
- formula
- moment
- streak camera
- 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.)
- Active
Links
- 230000035939 shock Effects 0.000 title claims abstract description 60
- 238000004364 calculation method Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000003384 imaging method Methods 0.000 claims abstract description 10
- 239000000523 sample Substances 0.000 claims description 23
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 claims description 20
- 230000003068 static effect Effects 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 9
- 230000003111 delayed effect Effects 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 3
- 230000005855 radiation Effects 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000003416 augmentation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000009514 concussion Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J11/00—Measuring the characteristics of individual optical pulses or of optical pulse trains
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种超时间分辨冲击波速度计算方法,包括以下步骤:S1:推导出冲击波速度与干涉条纹相位之间的解析关系;S2:推导出冲击波速度与干涉条纹相位之间的数值关系;S3:滤除数值方程解的噪声。采用以上方法,从波动方程解出发,精密证明条纹移动与速度变化之间的关系,为冲击波产生瞬间或追赶瞬间的速度计算提供了数值计算方法,将为研究极端条件下辐射冲击问题、流体力学不稳定性行为瞬时观测提供了与VISAR成像系统配套的数值计算方法。
Description
技术领域
本发明属于激光测量技术领域,具体涉及一种超时间分辨冲击波速度计算方法。
背景技术
冲击波的传播经常伴随着介质的压强、温度、密度等物理性质的跳跃式改变。测量得到冲击波速度将为物质状态方程求解提供边界条件,是广泛用于物态研究的物理参量。尤其在惯性约束聚变领域,烧蚀内爆过程产生的冲击波是一种高速流逝过程,其产生传播、反射、衰减、消失全过程只在纳秒、皮秒量级,且其冲击传播难免受到激光等离子体相互作用、流体力学不稳定性的影响。因此,具备更高冲击波速度求解方法将为此类物理过程观察提供途径,可以为研究者进一步完善基础物理理论模型,以及从事激光脉冲整形、调束提供可靠实验参数。
任意反射面速度干涉仪(VISAR)是用于冲击波速度信息采集的重要工具。其原理可以概括为:探针激光在冲击波界面发生反射,由于光学多普勒效应,反射光频率发生变化,因而承载冲击波速度信息;VISAR成像系统收集信息光,然后信息光在干涉仪中被分为两束,其中一束被延时,一束不被延时;最后两束光在相机阴极附近发生干涉,通过观察干涉条纹的移动便可以反推出冲击波速度。为满足不同条件下的诊断需要,VISAR装置在不断发展的同时,关于VISAR的数值计算分析也一直在完善。数值计算分析对象包含成像系统收光立体角、窗口材料、标准具的光学延时、噪声涨落等等。对于现有的VISAR技术,其时间分辨取决于系统本身的时间分辨(即标准具延迟时间,大于百皮秒)、像面记录装置的时间分辨、信号传输系统响应时间、记录装置的响应时间。随着工艺的不断改进,后面三项对VISAR装置时间分辨影响只有几皮秒,于是系统本身的时间分辨成为制约VISAR探测技术时间分辨提升的主要因素,经研究发现,系统本身的时间分辨其实是由现有算法的局限性导致,解决以上问题成为当务之急。
发明内容
为解决现有算法的局限性限制系统本身时间分辨的技术问题,本发明提供一种超时间分辨冲击波速度计算方法,大幅提高了冲击波速度时间分辨的计算精度。
为实现上述目的,本发明技术方案如下:
一种超时间分辨冲击波速度计算方法,包括以下步骤:
S1:推导出冲击波速度与干涉条纹相位之间的解析关系,具体步骤包括:
S11:表征条纹相机狭缝处的波动方程的解:
条纹相机狭缝处波动方程的解为
为干涉条纹振幅,Aj(x,y,t)为第j支路在条纹相机狭缝上的最大光振幅,i2=-1,t0为激光器发出光的时刻,为光源的初始相位,为物面上入射光与反射光之间的相位差,υ0为探针光的初始频率,Tj为支路j产生信息光的对应时刻,sj为Tj时刻光源与冲击波面之间的光程,λ0为探针光的初始波长,支路1为延时支路,支路2为非延时支路;
式(1)中
u为冲击波速度,c为光速,T条纹相机记录时刻,T0为物处于静态的任意时刻,sj为Tj时刻光源与冲击波面之间的光程,s0为T0时刻光源与物面之间的光程,Lj为Tj时刻条纹相机表面与物面之间的光程,L0为T0时刻条纹相机表面与物面之间的光程,τ为动态探测时探针光在标准具中的延迟时间,k为VISAR成像系统的纵向放大倍率,(xj,yj)为支路j在条纹相机狭缝处动态像的任意点坐标,(x0j,y0j)为支路j在条纹相机狭缝处静态像的坐标;
S12:表征相位方程解的形式:
光场强度I(x,y,t)满足以下关系:
A1为支路1在条纹相机狭缝上的最大光振幅,A2为支路2在条纹相机狭缝上的最大光振幅,由于实验条纹图样到计算使用的条纹图样,会经历傅里叶滤波过程,此时计算使用的条纹图样分布主要由相位变化决定,故A1=A2;
令x轴与条纹相机狭缝重合,将联立旋转变换公式以及式(1)代入式(2)得到的相位方程解的形式:
其中,(x0,y0)为像面传播方向与条纹相机狭缝完全垂直时的静态像面坐标,(0,m)为支路1与支路2的静态像面交点;
S13:将式(12)和式(13)进行化简得:
其中
令标准具修正量常数
则
h为标准具厚度,n为动态探测时探针光在标准具中的折射率,n0为静态探测时探针光在标准具中的折射率,τ0为静态探测时探针光在标准具中的延迟时间,λ为动态探测时探针光的波长;
对式(4)求微分得则
其中,为T1时刻的冲击波速度,为T2时刻的冲击波速度;
S14:用解析式表征冲击波速度与干涉条纹相位的关系:
对式(3)
上式中,由于中括号内的省略引起的相位变化远小于条纹相机记录的分辨率,中括号外的系数省略使得相位变化小于千分之一,因此在数值计算中直接使用约等号之后的形式,并在时间标度上对齐:
t为激光器发出光的时间变量,为t+τ0时刻的冲击波速度,ut为t时刻的冲击波速度;
S2:推导出冲击波速度与干涉条纹相位之间的数值关系,具体步骤包括:
S21:令τ0=4.378ι,ι为条纹相机的时间分辨率,先对式(5)求时间分辨上的算术平均,然后在整个时间域上求和并化简:
S22:由于在相位记录过程中,存在着数据涨落情况,故将涨落计入到相位值中,真实相位等于记录相位减去涨落,方程(6)的差分形式重新写为:
ω为随机涨落,将速度、相位值分别看成以时间变化的向量,上式写为:
其中d和p均为正整数;
故冲击波速度与干涉条纹相位的数值关系为
S3:滤除式(7)中解的噪声。
作为优选:步骤S3中,对式(7)即进行离散傅里叶变换:
其中Wm=e(-2πi)/m,m、f和ξ均为正整数,表示第ξ列的振幅,利用形如的窗口进行滤波。
采用以上方法,中,相当于对函数进行了加权,因此某一频谱的噪声在计算中被加强,利用形如的窗口进行滤波,故提出了一种新的速度谱的滤波方法,通过理论计算将噪声对速度谱影响的实质展现出来,并设计滤波窗口对速度谱进行去噪。
作为优选:光源相位在诊断时间内连续变化。因为需要保证在静态探测时具有稳定的干涉条纹,这样动态条纹相位的变化才仅与冲击波速度相关。既满足连续探测的要求,又在建立相位变化与速度的关系时有了实际意义。
与现有技术相比,本发明的有益效果是:
采用本发明提供的超时间分辨冲击波速度计算方法,从波动方程解出发,精密证明条纹移动与速度变化之间的关系,为冲击波产生瞬间或追赶瞬间的速度计算提供了数值计算方法,将为研究极端条件下辐射冲击问题、流体力学不稳定性行为瞬时观测提供了与VISAR成像系统配套的数值计算方法。
附图说明
图1为VISAR成像系统的原理图;
图2为实验数据解相位过程示意图;
图3为窗口滤波在速度处理中的演示实例;
图4为本发明的边值问题对数值计算的影响示意图;
图5为本发明与传统算法的速度分布示意图。
具体实施方式
以下结合实施例和附图对本发明作进一步说明。
如图1所示,VISAR成像系统的基本原理如下:
探针光通过分束镜照射在冲击波面上,由于光学多普勒效应,反射光携带冲击波速度信息,反射光经成像系统进入干涉仪,在干涉仪内被分为两束,其中一束被延时,一束不被延时,最后两束光在条纹相机狭缝附近发生干涉,干涉条纹移动反映速度信息,因此可以求解速度与条纹移动(相位变化)的关系。
一种超时间分辨冲击波速度计算方法,包括以下步骤:
S1:推导出冲击波速度与干涉条纹相位之间的解析关系,具体步骤包括:
S11:表征条纹相机狭缝处的波动方程的解:
条纹相机狭缝处波动方程的解为
为干涉条纹振幅,Aj(x,y,t)为第j支路在条纹相机狭缝上的最大光振幅,i2=-1,t0为激光器发出光的时刻,为光源的初始相位,为物面上入射光与反射光之间的相位差,υ0为探针光的初始频率,Tj为支路j产生信息光的对应时刻,sj为Tj时刻光源与冲击波面之间的光程,λ0为探针光的初始波长,支路1为延时支路,支路2为非延时支路,光源相位在诊断时间内变化连续;
式(1)中
u为冲击波速度,c为光速,T条纹相机记录时刻,T0为物处于静态的任意时刻,sj为Tj时刻光源与冲击波面之间的光程,s0为T0时刻光源与物面之间的光程,Lj为Tj时刻条纹相机表面与物面之间的光程,L0为T0时刻条纹相机表面与物面之间的光程,τ为动态探测时探针光在标准具中的延迟时间,k为VISAR成像系统的纵向放大倍率,(xj,yj)为支路j在条纹相机狭缝处动态像的任意点坐标,(x0j,y0j)为支路j在条纹相机狭缝处静态像的坐标;
S12:表征相位方程解的形式:
光场强度I(x,y,t)满足以下关系:
A1为支路1在条纹相机狭缝上的最大光振幅,A2为支路2在条纹相机狭缝上的最大光振幅,由于实验条纹图样到计算使用的条纹图样,会经历傅里叶滤波过程,此时计算使用的条纹图样分布主要由相位变化决定,故A1=A2。
请参见图2,图2为一次实验数据解相位过程,其中图2(a)为获取图样,图2(b)为取样,图2(c)为滤波,图2(d)为解相位。
令x轴与条纹相机狭缝重合,将联立旋转变换公式以及式(1)代入式(2)得到的相位方程解的形式:
其中,(x0,y0)为像面传播方向与条纹相机狭缝完全垂直时的静态像面坐标,(0,m)为支路1与支路2的静态像面交点;
以上存在同一时间标度下的三个时刻,条纹相机记录时刻T,波面1的作用时刻T1,波面2的作用时刻T2,需要将这三个时刻统一起来,以下是三者的关系:
S13:将式(12)和式(13)进行化简得:
其中
令标准具修正量常数
则
h为标准具厚度,n为动态探测时探针光在标准具中的折射率,n0为静态探测时探针光在标准具中的折射率,τ0为静态探测时探针光在标准具中的延迟时间,λ为动态探测时探针光的波长;
因此,当T1为均匀变化时,对式(4)求微分得:
则
其中,为T1时刻的冲击波速度,为T2时刻的冲击波速度,能够测到的均取VISAR成像系统能够测得的最大速度200km/s,△T和△τ均为通过VISAR成像系统测得的实验数据;
S14:用解析式表征冲击波速度与干涉条纹相位的关系:
对式(3)
可以计算
由于中括号内的省略引起的相位变化远小于条纹相机记录的分辨率,中括号外的系数省略使得相位变化小于千分之一,因此在数值计算中直接使用约等号之后的形式,只有当记录系统达到0.001的相位分辨时,该约等将会被重新考虑。
进一步在时间标度上对齐:
t为激光器发出光的时间变量,为t+τ0时刻的冲击波速度,ut为t时刻的冲击波速度;
S2:推导出冲击波速度与干涉条纹相位之间的数值关系,具体步骤包括:
S21:解的时间精度原则上不能小于条纹相机时间分辨,因此对式(5)求时间分辨上的算术平均,然后在整个时间域上求和并化简,其中令τ0=4.378ι,ι为条纹相机的时间分辨率。
S22:在真实的相位记录过程中,总是存在着数据涨落情况,这些可能来自于条纹相机记录响应的涨落,也有可能来自于噪声涨落,还有可能来自与相位解缠过程中引入的算法涨落等等。这些涨落相对于相位值来说很小,使用滤波的方法也难以消除,但是,它们会对最终的数值运算解产生影响,涨落中的特定频率会被迭代求解放大,导致最终速度解失真。将涨落计入到相位值中,真实相位等于记录相位减去涨落,式(6)的差分形式重新写为:
ω为随机涨落,将速度、相位值分别看成以时间变化的向量,上式写为:
其中d和p均为正整数;
故冲击波速度与干涉条纹相位的数值关系为
S3:滤除式(7)中解的噪声,具体如下:
对式(7)即进行离散傅里叶变换:
其中Wm=e(-2πi)/m,m、f和ξ均为正整数,的绝大多数列的离散傅里叶变换振幅相等。
请参见图3,图3(a)从上到下分别是图二采样的第200、400、600行的值,对其做快速傅里叶变换,傅里叶变换频谱如图3(b)。式(6)求出来的速度以及对速度的傅里叶变换频谱如图3(c)。图3(d)是两种频谱的对数坐标下的叠加。
因此,用第ξ列对振幅项进行表示,即ei2πζ(s)是叠加后的相位分布,ζ(s)≈0。当涨落能通过标定测量得到,式(8)可以精确求解。涨落引起的相位变化相对于真实φ的值,在动态过程可以忽略,而静态过程(静态条纹)对真实φ的获取影响大。但是静态过程的速度大小理想条件下是0,无需求解运算。
其中,相当于对函数进行了加权,故某一频谱的噪声在计算中被加强。因此,使用形如的窗口对第一项进行滤波。
请参见图4,在滤波过程中有两个问题需要注意:
一是最初的几行傅里叶变换频谱与“标准窗口”差异大,比如图4(a)是第10行的值与傅里叶变换结果。导致用标准谱滤波后最初时间和最后时间附近的速度值出现较大震荡,如图4(b),即该方法其实牺牲了边界值附近的精确速度解。
二是窗口的系数选择。系数值与随机变量Σωp相关。观察静态条纹中Σωp取值随时间的额变化,在前100个像素点内变化了约0.1,在前200个像素点内变化了约0.5,于是Σωp总体取值可在0.6-2.0之间。图4(b)列出了第60行几种不同的滤波情况下的最终速度曲线。在不能完全确定条纹相机的涨落情况下,不妨对Σωp取为1。
请参见图5,比较本发明提供的算法与传统算法的速度分布示意图,在不考虑新速度图的边值问题,图5(a)和图5(b)求得的速度基本一致,符合理论预期,但是区别在于新的数值计算方法将冲击加载一瞬间的速度分布求了出来(体现了时间分辨的提升),这为以后空间内辐射冲击转换为机械冲击、冲击波追赶或加载瞬间的速度分布求解提供了途径,甚至为精密观测流体力学不稳定行为提供了观测基础;为精密诊断冲击加载过程尤其是速度变化较大情况时的速度值精密求解提供了可能;可以精确观测得到具体哪个时段产生了不稳定性的增长,以及不稳定性的演化。
需要指出的是,本超时间分辨冲击波速度计算方法基于以下假设:
1.观测点反射回的中心光束频率替代返回的所有光的频率。实际收集到的每一点的返回探针光立体角在8度以内,计算该点返回频率无需进一步校正。
2.波动方程解建立在条纹相机的狭缝附近。因为干涉条纹形成区域是所关心的。
3.狭缝附近的光场传播方向即为像面传播方向。冲击压缩中,像面始终处于景深范围内,因而像面相对位置不产生变化,但相位分布随时间和位置变化。
最后需要说明的是,上述描述仅仅为本发明的优选实施例,本领域的普通技术人员在本发明的启示下,在不违背本发明宗旨及权利要求的前提下,可以做出多种类似的表示,这样的变换均落入本发明的保护范围之内。
Claims (2)
1.一种超时间分辨冲击波速度计算方法,其特征在于,包括以下步骤:
S1:推导出冲击波速度与干涉条纹相位之间的解析关系,具体步骤包括:
S11:表征条纹相机狭缝处的波动方程的解:
条纹相机狭缝处波动方程的解为
为干涉条纹振幅,Aj(x,y,t)为第j支路在条纹相机狭缝上的最大光振幅,i2=-1,t0为激光器发出光的时刻,为光源的初始相位,为物面上入射光与反射光之间的相位差,υ0为探针光的初始频率,Tj为支路j产生信息光的对应时刻,sj为Tj时刻光源与冲击波面之间的光程,λ0为探针光的初始波长,支路1为延时支路,支路2为非延时支路;
式(1)中
u为冲击波速度,c为光速,T条纹相机记录时刻,T0为物处于静态的任意时刻,sj为Tj时刻光源与冲击波面之间的光程,s0为T0时刻光源与物面之间的光程,Lj为Tj时刻条纹相机表面与物面之间的光程,L0为T0时刻条纹相机表面与物面之间的光程,τ为动态探测时探针光在标准具中的延迟时间,k为VISAR 成像系统的纵向放大倍率,(xj,yj)为支路j在条纹相机狭缝处动态像的任意点坐标,(x0j,y0j)为支路j在条纹相机狭缝处静态像的坐标;
S12:表征相位方程解的形式:
光场强度I(x,y,t)满足以下关系:
A1为支路1在条纹相机狭缝上的最大光振幅,A2为支路2在条纹相机狭缝上的最大光振幅,由于实验条纹图样到计算使用的条纹图样,会经历傅里叶滤波过程,此时计算使用的条纹图样分布主要由相位变化决定,故A1=A2;
令x轴与条纹相机狭缝重合,将联立旋转变换公式以及式(1)代入式(2)得到的相位方程解的形式:
其中,(x0,y0)为像面传播方向与条纹相机狭缝完全垂直时的静态像面坐标,(0,m)为支路1与支路2的静态像面交点;
S13:将式(12)和式(13)进行化简得:
其中
令标准具修正量常数
则
h为标准具厚度,n为动态探测时探针光在标准具中的折射率,n0为静态探测时探针光在标准具中的折射率,τ0为静态探测时探针光在标准具中的延迟时间,λ为动态探测时探针光的波长;
对式(4)求微分得则
其中,为T1时刻的冲击波速度,为T2时刻的冲击波速度;
S14:用解析式表征冲击波速度与干涉条纹相位的关系:
对式(3)求微分
上式中,由于中括号内的省略引起的相位变化远小于条纹相机记录的分辨率,中括号外的系数省略使得相位变化小于千分之一,因此在数值计算中直接使用约等号之后的形式,并在时间标度上对齐:
t为激光器发出光的时间变量,为t+τ0时刻的冲击波速度,ut为t时刻的冲击波速度;
S2:推导出冲击波速度与干涉条纹相位之间的数值关系,具体步骤包括:
S21:令τ0=4.378ι,ι为条纹相机的时间分辨率,先对式(5)求时间分辨上的算术平均,然后在整个时间域上求和并化简:
S22:由于在相位记录过程中,存在着数据涨落情况,故将涨落计入到相位值中,真实相位等于记录相位减去涨落,方程(6)的差分形式重新写为:
ω为随机涨落,将速度、相位值分别看成以时间变化的向量,上式写为:
其中d和p均为正整数;
故冲击波速度与干涉条纹相位的数值关系为
S3:滤除式(7)中解的噪声:
对式(7)即进行离散傅里叶变换:
其中m、f和ξ均为正整数,表示第ξ列的振幅,利用形如的窗口进行滤波。
2.根据权利要求 1所述的超时间分辨冲击波速度计算方法,其特征在于:光源相位在诊断时间内连续变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810839877.8A CN108871595B (zh) | 2018-07-27 | 2018-07-27 | 超时间分辨冲击波速度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810839877.8A CN108871595B (zh) | 2018-07-27 | 2018-07-27 | 超时间分辨冲击波速度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108871595A CN108871595A (zh) | 2018-11-23 |
CN108871595B true CN108871595B (zh) | 2019-11-15 |
Family
ID=64306078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810839877.8A Active CN108871595B (zh) | 2018-07-27 | 2018-07-27 | 超时间分辨冲击波速度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108871595B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112215770B (zh) * | 2020-10-10 | 2022-08-02 | 成都数之联科技股份有限公司 | 一种图像处理方法及系统及装置及介质 |
CN112927824B (zh) * | 2021-01-27 | 2024-01-19 | 中国人民解放军火箭军工程大学 | 基于勒让德模型的广角visar条纹分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093962A (zh) * | 2016-08-03 | 2016-11-09 | 中国工程物理研究院流体物理研究所 | 一种干涉测速系统及方法 |
CN107179132A (zh) * | 2017-05-09 | 2017-09-19 | 中国工程物理研究院激光聚变研究中心 | 光纤传像束速度干涉仪及冲击波速度计算方法 |
CN108036863A (zh) * | 2017-12-19 | 2018-05-15 | 中国工程物理研究院激光聚变研究中心 | 宽量程冲击波速度诊断装置及测量方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0271111A (ja) * | 1988-09-07 | 1990-03-09 | Kokusai Electric Co Ltd | 超音速飛行体の通過位置検知方法 |
US10018647B2 (en) * | 2014-05-16 | 2018-07-10 | Defence Research Development Organisation | Velocity interferometer for any reflector with variable sensitivity range and time resolution |
-
2018
- 2018-07-27 CN CN201810839877.8A patent/CN108871595B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093962A (zh) * | 2016-08-03 | 2016-11-09 | 中国工程物理研究院流体物理研究所 | 一种干涉测速系统及方法 |
CN107179132A (zh) * | 2017-05-09 | 2017-09-19 | 中国工程物理研究院激光聚变研究中心 | 光纤传像束速度干涉仪及冲击波速度计算方法 |
CN108036863A (zh) * | 2017-12-19 | 2018-05-15 | 中国工程物理研究院激光聚变研究中心 | 宽量程冲击波速度诊断装置及测量方法 |
Non-Patent Citations (1)
Title |
---|
预热效应影响下冲击波速度推算方法;王峰等;《光学学报》;20100531;第30卷(第5期);第1327-1331页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108871595A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109298457B (zh) | 适用于原子干涉重力仪的振动噪声修正补偿方法 | |
CN108871595B (zh) | 超时间分辨冲击波速度计算方法 | |
CN104006765A (zh) | 单幅载频干涉条纹相位提取方法及检测装置 | |
CN102175580A (zh) | 动态散斑法测量浑浊介质微粒运动的装置和方法 | |
CN103592108A (zh) | Ccd芯片调制传递函数测试装置及方法 | |
Limacher et al. | On the calculation of force from PIV data using the generalized added-mass and circulatory force decomposition | |
Kaufmann et al. | Reconstruction of surface-pressure fluctuations using deflectometry and the virtual fields method | |
Zaletelj et al. | Full-field DIC-based model updating for localized parameter identification | |
Zhang et al. | A high precision signal processing method for laser Doppler velocimeter | |
Gordeyev et al. | Aero-Optical Effects: Physics, Analysis and Mitigation | |
CN106383367A (zh) | 绝对重力测量方法和装置 | |
Kaufmann et al. | Reconstruction of surface pressures on flat plates impacted by blast waves using the Virtual Fields Method | |
CN111610572B (zh) | 一种原子干涉重力仪波前相移快速评估方法及系统 | |
Cadel et al. | Investigation of measurement sensitivities in cross-correlation Doppler global velocimetry | |
Haonian et al. | Method improving low Signal-to-noise ratio of velocity test signals for Laser-induced shock waves | |
CN103698022A (zh) | 横向剪切干涉仪波前测量方法 | |
CN103267485A (zh) | 一种点衍射三维绝对位移测量方法 | |
CN101285711A (zh) | 基于面阵ccd的线性相位反演波前传感器 | |
Rasouli et al. | Applications of 2-D moiré deflectometry to atmospheric turbulence | |
Pfeifer et al. | Application of the Tomographic Laser Doppler Anemometry (TDLA) to a fuel spray | |
Dai et al. | Non-destructive strain determination based on phase measurement and radial basis function | |
Zhang et al. | Signal processing of double-beam and double-scattering laser Doppler velocimeter | |
Wang et al. | Sequential observation of rebound shock wave generated by collapse of vapor bubble in BOS system | |
CN202024946U (zh) | 动态散斑法测量浑浊介质微粒运动的装置 | |
CRAIG et al. | Propagation diagnostic technique for turbulent transonic flow |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |