海洋-声学耦合模式捕捉水声环境不确定性

考虑到海洋环境的时空变化对水声环境不确定性的影响,建立了海洋-声学耦合数值模式,实现了并行计算。该模式将声学计算纳入到运动的海洋中,从而实现了对水声环境的动态预报和估计。同时,采用集合预报方法对典型断面的温度垂直结构、实验海区声速剖面和传播损失进行预报,并给出了声速剖面的预报误差、不同深度与频率下,传播损失90%的概率区间以及声速、传播损失、声呐作用距离的不确定性直方图。结果反映了海洋时空变化对水声环境不确定性的影响,量化了水声环境中不确定性的大小。实验结果表明该方法可以刻画海洋动态变化引起的水声环境不确定性,并对其进行了量化和描述。      

三维复杂温度场高精度声学测量方法

声学温度场检测技术通过多路径声波传播时间数据,反演被测区域的温度分布。提供了一种高精度的三维复杂温度场的声学测量方法。首先从射线声学角度给出了三维非均匀温度场中声波传播路径的数学模型。在此基础上,将三维温度场的重建问题转化为声波传播路径的求解和温度场的反演问题,建立了基于多项式修正径向基函数(RBF-PR)和改进的Tikhonov正则化三维温度场重建算法。采用两种典型的炉膛三维温度场模型,在信噪比SNR=35 dB下进行了数值模拟,分析了声波传播路径在非均匀温度场中的弯曲特性、算法的重建质量和抗噪性,同时进行了实际炉膛内二维温度场的重建。结果表明了提出的考虑声线弯曲的温度场重建算法具有精度高,抗噪性强、适用性好的特点,为实现高精度的复杂温度场的声学测量提供了有效方法。      

海洋-声学耦合模式捕捉水声环境不确定性

考虑到海洋环境的时空变化对水声环境不确定性的影响,建立了海洋-声学耦合数值模式,实现了并行计算。该模式将声学计算纳入到运动的海洋中,从而实现了对水声环境的动态预报和估计。同时,采用集合预报方法对典型断面的温度垂直结构、实验海区声速剖面和传播损失进行预报,并给出了声速剖面的预报误差、不同深度与频率下,传播损失90%的概率区间以及声速、传播损失、声呐作用距离的不确定性直方图。结果反映了海洋时空变化对水声环境不确定性的影响,量化了水声环境中不确定性的大小。实验结果表明该方法可以刻画海洋动态变化引起的水声环境不确定性,并对其进行了量化和描述。     

界面参数估计及传播损失的射线声学预报

利用射线声学建立的海洋信道模型,根据给定海区声传播实验的实测数据,反演得到界面声学参数,并进行声传播损失预报。预报结果与实验数据比较表明,信道模型正确,界面参数选取与估计方法合理,传播损失预报性能良好。     

声学CT复杂温度场重建研究

为提高声学CT复杂温度场重建能力,提出一种利用Markov径向基函数逼近和Tikhonov正则化的温度场重建算法,简称MTR算法。该算法首先用Markov径向基函数的线性组合,逼近介质中的复杂声速场分布,然后利用介质中多路径声波传播时间和Tikhonov正则化法,求解声速场分布,进而利用声速与温度的关系获得温度分布。对单热点、三热点和五热点温度场模型进行了仿真重建,结果表明MTR算法热点定位精度高,重建误差小。开发了声学CT温度场重建实验系统,用电加热器在内装1200 kg大豆的实验粮仓中形成热点,MTR重建结果能正确反映热点位置,热点温度重建误差1.3%。可见,MTR算法复杂温度场重建能力强,可望用于实际储粮温度分布监测。      

动态水声环境不确定性的预报与估计

水声环境具有强烈的时空易变性,为了解其变化中的规律和产生的不确定性,将POM海洋数值模型和水声传播模型FOR3D进行耦合,建立海洋—声学耦合数值模式。对海区的温盐等环境参数和水下声场进行预报,给出了典型断面的温度垂直结构、实验海区声速剖面和传播损失。同时,采用集合卡尔曼滤波(EnKF)方法结合海洋—声学耦合数值模式分析了水声环境的不确定性,对不同深度与频率下,传播损失的均值与满足90%的可能区间进行了计算,同时给出了不同点声速、传播损失、声纳作用距离的不确定性直方图。和实验数据的对比结果表明该方法具有较高的有效性,能对水声环境的不确定性进行合理的预报与估计。     

海洋声学传播新模型

海洋声学模型将我们对海洋中传播的声音在物理上的理解转换为数学表示,进而利用计算机进行模拟计算。海洋声学模型主要包括环境模型、传播模型、噪声模型、混响模型和声纳性能模型。水平变化传播模型研究,包括二维(距离和深度)和三维(距离、深度和方位角)模型研究,是目前海洋声学的重要研究方向之一。     

海底反射损失的高分辨率估计方法

在海洋声学中,海底反射损失是一个重要的物理量,不仅可以用于传播损失的计算,也可以用来反演海底参数(如声速、密度、吸收系数等)。海底反射损失与声波频率和掠射角有关,近年来出现了多种估计方法,其中一种简单的被动方法是对环境噪声做垂直波束形成,来自海底与来自海面的噪声比例能够反映出海底反射损失的大小[1]。该方法的一个主要优点是能够直接得到反射损失,不需要任何反演机制。但有一个明显的缺点,当垂直接收     

海洋声学环境下水中声源的时延估计法定位精度分析

水中声源的定位精度受到海洋声学环境的重要影响。结合海上试验的实际应用,分析了水下观测平台采用时延估计法对声源的定位精度问题。根据理论分析,计算了时延估计误差、海洋中声速不均匀、平台非稳性、及声传播起伏等因素引起的俯仰角和方位角误差。利用误差传递公式,获得了上述因素引起的不同平台深度下,不同距离声源的定位误差。比较了采用平面阵与立体阵、是否补偿声线弯曲效应等条件下定位误差的变化,并通过海上试验结果进行了部分验证。研究结果表明,海洋声速不均匀对定位误差的贡献最大。采用立体阵代替平面阵、测量海洋声速剖面并补偿声线弯曲引起的定位误差,在1000m距离上可使定位相对误差从最大30%降低到约10%,有效提高了较远距离上的定位精度。研究结果对于采取措施提高水中声源的定位精度有指导意义。      

重塑砂质含气沉积物的声响应特性与模型预测分析

海洋沉积物中气泡的存在对沉积物的声学特性有显著影响.为实现在实验室内对不同初始应力条件下含气沉积物的声学特性测量,研制了一套可用于土工三轴仪的双探针声学测量系统,基于CT扫描试验获得重塑含气砂样品中气泡群的尺寸分布,确定其共振频率范围.试验结果表明:细砂沉积物声速随气体含量的增加呈指数型降低,衰减系数随气体含量的增加呈...     

Propagation modeling of ocean-scattered low-elevation GPS signals for maritime tropospheric duct inversion

The maritime tropospheric duct is a low-altitude anomalous refractivity structure over the ocean surface,and it can significantly affect the performance of many shore-based/shipboard radar and communication systems.We propose the idea that maritime tropospheric ducts can be retrieved from ocean forward-scattered low-elevation global positioning system(GPS) signals.Retrieval is accomplished by matching the measured power patterns of the signals to those predicted by the forward propagation model as a function of the modified refractivity profile.On the basis of a parabolic equation method and bistatic radar equation,we develop such a forward model for computing the trapped propagation characteristics of an ocean forward-scattered GPS signal within a tropospheric duct.A new GPS scattering initial field is defined for this model to start the propagation modeling.A preliminary test on the performance of this model is conducted using measured data obtained from a 2009-experiment in the South China Sea.Results demonstrate that this model can predict GPS propagation characteristics within maritime tropospheric ducts and serve as a forward model for duct inversion.     

光纤中扰动的小信号增益

从非线性薛定谔方程出发,在小信号近似下,推导并求解了光纤中扰动相位和幅度的演化方程,利用得到的扰动相位及功率增益的表达式,研究了初相位和频率对传输过程中扰动增益的影响。研究表明：扰动的初相位对扰动增益的初值和初始阶段的演化规律有重要影响;取决于扰动初相位,任何一个频率的扰动增益都有可能达到一个共同的最大值;在被认为无调制不稳定的正色散区和扰动频率大于截止频率的负色散区,扰动增益随距离是振荡的;在被认为有调制不稳定的扰动频率小于截止频率的负色散区,频率相同而初相位不同的扰动增益将经历不同形式的演化后趋于同一正值。     

Shock wave formation in hot nuclear matter

Assuming that nuclear matter can be treated as a perfect fluid, we study the propagation of perturbations in the baryon density at high temperature. The equation of state is derived from the non-linear Walecka model. The expansion of the Euler and continuity equations of relativistic hydrodynamics around equilibrium configurations lead to the breaking wave equation for the density perturbation. We solve it numerically for this perturbation and follow the propagation of the initial pulses.     

Kinetics of an irreversible transition to a normal state by a multiconductor superconducting composite

A model that takes into account the discrete character of a multiconductor superconducting medium is used to study the characteristic features of transition processes which occur when a region of normal conductivity appears. This analysis is based on numerical simulations of the propagation of a normal zone inside a superconducting region consisting of rectilinear superconductors, with a model for the kinetics of the irreversible loss of the superconducting properties of a single-ply, large radius winding. The thermophysical aspects of the irreversible propagation of a thermal perturbation within the transverse cross section of a multiconductor superconducting medium are analyzed. The major features of the phenomena taking place during the loss of superconductivity by a discrete region are formulated as functions of the conditions for heat transfer between its elements, the character of the perturbation, the magnitude of the flowing current, and the conditions for stabilization. The influence of longitudinal heat conduction on the magnitude of the transverse propagation velocity of a normal zone is examined. Zh. Tekh. Fiz. 69, 40–51 (May 1999)     

Simulation of the destruction of a woodland massif by the shock wave generated by the impact of a large cosmic body

The interaction of shock waves with woodland massifs with consideration given to vegetation destruction is examined. The results of numerical simulations of the propagation of a perturbation through a vegetation massif based on solving the Euler equations for two-dimensional and axisymmetric cases are presented. This approach was used by the present author in simulating the propagation of wave fronts through woodland massifs without regard for forest stand destruction. The effect of vegetation is reproduced by introducing bulk forces in the volume occupied by forest. An initial perturbation of explosion-in-flight type was specified in the form of an increased-pressure region located at various heights from the ground surface. Simple threshold models are used to examine the effect of the dynamic pressure on vegetation destruction. A refined model was developed to investigate the change in the configuration of woodland massifs with consideration given to vegetation fall and the effect of this process on the rate of perturbation dissipation. The system of equations was solved by the Godunov method implemented using parallel programming. Numerical experiments were conducted with the use of two types of models of forest stand destruction, with various threshold values of the destruction factor, various initial perturbation intensities, and various heights above the vegetation massif. The effects arising during the interaction of a high-intensity perturbation with a woodland massif were identified. An analysis of the results showed that the destruction of vegetation occurs during its interaction with the fronts of the incident and reflected shock wave and with the jet core of the perturbation. In some cases, the model taking into account the accumulation of fallen vegetation predicts situations where destruction will occur only in the upper part of the forest stand, without touching the canopy. This effect can explain the existence of such regions in the area of the Tunguska event. It was also demonstrated that estimates of the sizes of the destruction zones obtained with the use of the models taking into account the influence of vegetation on the perturbation and destructibility differ substantially from those obtained within the framework of the models without regard for vegetation resistance. The highest range of perturbation dissipation was predicted by the models simulating the formation of obstructions.     

发散光束小尺度自聚焦特性的研究

研究了发散光束的小尺度自聚焦效应。从非线性傍轴波动方程出发,利用坐标变换,推导出发散光束小尺度扰动的传输方程,进而得到小尺度扰动增长的临界频率、最大增长频率和相应B积分值的变化规律。研究了发散光束初始曲率半径对小尺度自聚焦效应的影响。结果表明,对于一定的传输距离,随着发散光束初始曲率半径的减小,小尺度扰动的最大增长频率减小,相应的最大增益减小,即B积分值也减小。对于一定的初始曲率半径,随着传输距离的增大,B积分值增长变缓,并最终停止。利用局部能量守恒定律研究了发散光束的成丝距离,发现小的初始曲率半径可以延长成丝距离。     

Blowup of Smooth Solutions for Relativistic Euler Equations

We study the singularity formation of smooth solutions of the relativistic Euler equations in (3 + 1)-dimensional spacetime for both finite initial energy and infinite initial energy. For the finite initial energy case, we prove that any smooth solution, with compactly supported non-trivial initial data, blows up in finite time. For the case of infinite initial energy, we first prove the existence, uniqueness and stability of a smooth solution if the initial data is in the subluminal region away from the vacuum. By further assuming the initial data is a smooth compactly supported perturbation around a non-vacuum constant background, we prove the property of finite propagation speed of such a perturbation. The smooth solution is shown to blow up in finite time provided that the radial component of the initial ``generalized' momentum is sufficiently large.     

Simulation of Internal Gravity Wave Propagation Due to Sudden Stratospheric Warming

The results of modeling the effect of stratospheric internal gravity waves (IGWs) excited in the region of sudden stratospheric warming (SSP) on the state of the upper atmosphere are presented. The numerical experiment used a two-dimensional model of atmospheric wave propagation including the dissipative and nonlinear processes accompanying the wave propagation. As a perturbation source, perturbations of temperature and density in a localized region at stratospheric altitudes during the SSW were considered. The amplitude and frequency of the perturbation source were evaluated from the results of observations and IGW theory. The results of numerical calculations showed that the heat source localized in the stratosphere excites IGWs with periods of up to several hours, which reach the thermospheric altitudes within a few hours. The maximum relative perturbations created by these waves relative to the unperturbed conditions were noted at altitudes of 100–200 km and at distances of up to ~1000 km from the source center.     

Multifilamentation of femtosecond laser pulses propagating in?turbulent air near the ground

The propagation of femtosecond laser pulses in turbulent air near the ground is analyzed. Confining to a power regime distinctly above the critical power for self-focusing, i.e. P≈100P _cr, and concentrating on initial peak intensities around 2.5×10¹¹W/cm², the onset and early evolution of multiple filaments are addressed. Making use of the turbulence phase-screen method, numerical simulations of the pulse propagation indicate that turbulence fields with spatial scales below 6 mm are able to induce the onset of multifilamentation. An analytical linear plane wave perturbation model of the underlying modulation instability of the pulse front is introduced in support of the computational results. By this means, insight into the amplification of an initial perturbation of the pulse front from the point of view of the spatial frequency domain is given.