高斯涡旋光束的光束传输因子和峭度参数

基于强度二阶矩定义, 导出了高斯涡旋光束光束传输因子即M² 因子的解析表达式, 高斯涡旋光束的M² 因子唯一取决于拓扑电荷数n. 数值计算表明, 高斯涡旋光束的M² 因子随着拓扑电荷数n的增大而增大. 基于强度高阶矩, 还导出了高斯涡旋光束经傍轴ABCD光学系统传输时峭度参数的解析表达式, 高斯涡旋光束的峭度参数取决于拓扑电荷数n、参数δ、矩阵元A和矩阵元D. 在自由空间传输时, 高斯涡旋光束的峭度参数仅取决于拓扑电荷数n和参数δ. 自由空间传输时, 高斯涡旋光束峭度参数的变化规律为: 峭度参数随参数δ的增大先减小而后趋向于一最小值, 随拓扑电荷数n的增大而减小. 这一研究有助于高斯涡旋光束的实际应用.     

格子Boltzmann方法模拟填充纳米流体三角形腔内的自然对流

采用格子Boltzmann方法研究填充水-氧化铝纳米流体的等腰直角三角形腔体中的自然对流。讨论瑞利数、颗粒体积分数、热源位置等因素对对流换热的影响,以及不同纳米流体模型对模拟结果的影响。结果表明：在低瑞利数下,随着热源在左壁面向上移动,换热效率逐渐增加。而在高瑞利数（Ra=10⁶）时,观察到相反的现象;采用单相纳米流体模型,模拟表明热壁面平均努塞尔数比率随着体积分数的增加近似线性增加。采用改进的纳米流体模型,结果显示平均努塞尔数比率随着体积分数的增加而增大,但是平均努塞尔数比率的变化斜率逐渐减小。改进模型模拟的换热效率比单相模型高,这是因为改进模型考虑了粒子间作用力及换热,更符合实际情况。     

Er3+：GdVO4中Er3+离子的光谱参数计算和晶场中能级分裂的讨论

对Er³⁺：GdVO₄样品的光谱参数以及Er³⁺在晶场中能级的分裂情况进行了研究. 首先对样品进行了吸收光谱的测量，接着用Judd-Ofelt理论拟合出了Er³⁺在GdVO₄晶体中的强度参量Ω_t，并由此计算了跃迁的振子强度、自发辐射跃迁速率、荧光分支比和积分发射截面. 通过计算结果可以发现有较多能级之间的跃迁都有大于10^-6的振子强度和大于10^-18cm的积分发射截面，并且具有较高的荧光分支比，特别是²H_11/2→⁴I_15/2，⁴S_3/2→⁴I_15/2，⁴F_9/2→⁴I_15/2和⁴I_13/2→⁴I_15/2等几个强发光能级除了具有较大的振子强度和积分发射截面外还有很好的应用前景，因此也更加值得关注. 最后还利用群论讨论了Er³⁺离子在GdVO₄晶场中各能级的分裂情况并对各Stark子能级的J_z混杂情况进行了分析.     

飞秒激光激发空气电离的阈值研究

报道在脉宽50fs—22ps，波长800nm脉冲激光作用下的空气电离阈值的研究结果.利用探测等离子体发光信号的方法，实验测量了激发空气电离所需的阈值激光强度.结果表明，当激光脉冲宽度从50fs增加到22ps时，阈值光强I_th从8.7×10¹⁴W/cm²下降到2.7×10¹³W/cm²；I_th经历了由迅速降低逐渐发展为缓慢降低的过程.在50fs—1p     

基于交通流量的病毒扩散动力学研究

不同于经典扩散模型中节点传染力等同于节点度k的假定, 基于交通流量的病毒扩散模型中, 各个节点的传染力可以等同于节点实际介数b_k. 利用平均场近似方法, 提出基于交通流量SIS病毒修正扩散模型. 根据修正SIS模型, 以最小搜索信息路由为例, 重新研究病毒传播率β, 平均发包率λ同传播阈值β_c, 平稳状态病毒密度ρ之间的关系. 理论分析与实验结果均表明, 当网络拓扑和路由策略一定时, 传播阈值β_c为实际介数b_k的均值<b_k>与其平方的均值<b_k²>的比值. 而稳定状态时感染密度ρ同感染同病毒传播率β, 平均发包率λ 以及λ =1时节点实际介数的均值<b_λ=1> 的乘积倒数存在幂率关系.     

第Ⅵ类两态叠加多模叠加态光场的等阶N次方Y压缩特性研究

本文构造了由多模相干态|{Z_j}〉q与多模虚相干态|{iZ_j}〉q这两者的线性叠加所组成的第Ⅵ类两态叠加多模叠加态光场|Ψ₆⁽²⁾〉q.利用多模压缩态理论,研究了态|Ψ₆⁽²⁾〉q的广义非线性等阶N次方Y压缩特性.结果发现:1)在压缩阶数N＝2p且p＝2l(l＝1,2,3,…,…)的条件下,无论各模的初始相位φ_j(j＝1,2,3,…,q)、态间的初始相位差(θ_pq^(R)-θ_pq^(I))以及各单相干态光场的平均光子数之和 R_j²等如何变化,态|ψ₆⁽²⁾〉q总是恒处于偶数阶等阶N-Y最小测不准态.2)在压缩阶数N＝2p且p＝2l+1(l＝0,1,2,3,…,…)的条件下,当各模的初始相位φj、态间的初始相位差(θ_pq^(R)-θ_pq^(I))以及各单模相干态光场的平均光子数之和 R_j²等分别满足一定的量子化条件(或者在特定的区间内连续取值)时,态|ψ₆⁽²⁾〉q的第一及第二这两个正交分量总可分别呈现出周期性变化的偶数阶的等阶N次方Y压缩效应.     

薄膜型超导红外(光)探测器

在高Tc GdBaCuO超导薄膜上,采用光刻技术分别制成两种不同结构的辐射、热测量器件及2×4集成阵列式微桥器件红外(光)探测器.探测器芯片安装在STD-3型红外探测器杜瓦冷指上.用黑体及波长为0.6328μm的He-Ne激光器辐照器件,系统观测各种器件的特性,其中最好的结果:在10Hz时的噪音等效功率NEP(500,10,1)=3.6×10^-12 WH_z^1/2;探测率D^*(500,10,1)=1.6×10¹⁰ cmH_z^1/2W^-1;响应率R_v=8.2×103VW^-1.另外,多元串接微桥器件出现的多台阶式的特性,可望在红外探测计量及高频方面获得重要应用.     

激光光控超导YBa2Cu3Ox薄膜开关的研制

在对光控热电效应开关进行理论分析的基础上,木文提出用YBa₂Cu₃O_x薄膜制作光控开关,并测试了在液氮温度下薄膜开关在不同激光波长下的特征参数,测试的最好结果是响应度R_v(632.8nm,10kHz,1Hz)为217V/W,归一化探测率D_*(632.8nm,10kHz,1Hz)为2.3×10¹¹cm.Hz^1/2/W,响应时间τ为0.21ms.     

光学微腔中少光子数叠加态的耗散动力学

通过考察耗散光学腔中少光子数叠加态的Wigner函数随时间 的变化行为, 揭示其非经典特性的动力学演化. 结果表明, 初始时Wigner函数为负的少光子数叠加态, 在耗散过程中其负性逐渐减小 直至消失, 并最后达到一个稳定的正值. 但这并不意味着耗散量子态非经典特性的完全消失. 实际上, 作为非经典特性的另一个重要参量, 光子的二阶关联函数g⁽²⁾(0) (g⁽²⁾(0)<1意味着光子是非经典地反聚束) 是一个耗散动力学不变的物理量. 我们证明, 光子的二阶反关联函数g^(2A)(0)却是一个随着耗散而改变的物理参量, 从而可以用于描述光学微腔中光量子态的耗散动力学行为. 最后, 我们给出一个在实验上如何制备少光子数叠加态并对其Wigner函数进行探测的方案.     

超短脉冲强激光与固体靶相互作用中Kα射线的实验研究

在相对论激光强度下，对p偏振30 fs激光与固体Cu靶相互作用中产生的Kα射线进行了实验研究.采用刀边成像技术和单光子计数X射线CCD相结合的探测装置，在单发激光脉冲打靶时同时得到X射线源的尺寸、能谱以及Kα光子的转换效率等多种信息.实验结果与Reich等人的理论计算结果有明显的差异，Kα光子的能量转换效率在激光功率密度为1.6×10¹⁸W/cm²的条件下达到最大值7.08×10^-6/sr.根据这一结果并结合蒙特卡罗程序，推断出在这一聚焦光强下激光能量转换为前向超热电子的效率约为10%.     

Scaling in convective evaporation and sidewall boundary layer

Turbulent convection at aspect ratios from 0.06 to 2 is investigated in the laboratory with evaporation experiments from vertical cylinders having different diameters and liquid levels. With alcohol, only diffusive evaporation takes place. With water, for small diameters, evaporation proceeds by diffusion whereas convective evaporation develops when the diameter is increased. This onset can be effectively interpreted in terms of a viscous sidewall boundary layer, whose thickness δ varies with respect to the available height h according to δ/h = 3.4 Ra^-0.28±0.01 versus Rayleigh number Ra. The Sherwood number Sh, analog of the Nusselt number, exhibits a power law variation Sh = 0.6 Ra^0.27±0.02 for Ra varying from 10⁴ to 3 ×10⁸. The scaling observed in this case of an open boundary is thus similar to the scaling measured in confined Rayleigh-Bénard convection.     

Application of the multi distribution function lattice Boltzmann approach to thermal flows

Numerical methods able to model high Rayleigh (Ra) and high Prandtl (Pr) number thermal convection are important to study large-scale geophysical phenomena occuring in very viscous fluids such as magma chamber dynamics (10⁴ < Pr < 10⁷ and 10⁷ < Ra < 10¹¹). The important variable to quantify the thermal state of a convective fluid is a generalized dimensionless heat transfer coefficient (the Nusselt number) whose measure indicates the relative efficiency of the thermal convection. In this paper we test the ability of Multi-distribution Function approach (MDF) Thermal Lattice Boltzmann method to study the well-established scaling result for the Nusselt number (Nu ∝ Ra ^1/3) in Rayleigh Bénard convection for 10⁴ ≤ Ra ≤ 10⁹ and 10¹ ≤ Pr ≤ 10⁴. We explore its main drawbacks in the range of Pr and Ra number under investigation: (1) high computational time N _c required for the algorithm to converge and (2) high spatial accuracy needed to resolve the thickness of thermal plumes and both thermal and velocity boundary layer. We try to decrease the computational demands of the method using a multiscale approach based on the implicit dependence of the Pr number on the relaxation time, the spatial and temporal resolution characteristic of the MDF thermal model.     

Long relaxation times and tilt sensitivity in Rayleigh Bénard turbulence

We study the influence of a small tilt angle ( rd) on the Nusselt number in a 1/2 aspect ratio Rayleigh-Bénard cell, at high Rayleigh number (5 x 10¹¹ < Ra < 4 x 10¹²). The small decrease observed is interpreted as revealing a two rolls structure of the flow. Transitions between different global flows are also observed, on very long times, comparable to the diffusion time on the whole cell. The consequence is that the Nusselt number observed in most high Ra experiments should significantly depend on initial conditions.Received: 19 May 2004, Published online: 12 August 2004PACS: 92.60.Ek Convection, turbulence, and diffusion - 47.27.Te Convection and heat transfer - 44.25. + f Natural convection     

Temperature gradient spectra and temperature dissipation rate in a turbulent convective flow

Determining mixing coefficients in oceanographic flows relies on the form of temperature gradient spectra in turbulent water flows at large wavenumbers. Several recent investigations concluded that these spectra are best described by the functional form proposed by Kraichnan rather than by the Batchelor form, more commonly used in oceanography. In this study, we provide additional support for this conclusion using laboratory measurements of the temperature gradient spectra in a Rayleigh–Bénard convective cell, in order to avoid difficulties inherent in oceanographic field measurements. The range of Rayleigh numbers in experiments is between Ra = 3×10⁷ and Ra = 5×10⁹. In addition to a traditional method of traversing thermistors, a novel optical technique recently introduced for oceanic measurements was used to obtain the spectra; comparison between these two methods serves as a validation test for the new optical technique. The temperature measurements were also augmented by 2D particle image velocimetry (PIV) observations. The measured dependence of the Nusselt number on the Rayleigh number followed Nu ∝ Ra^0.29 at Pr = 6 and was consistent with the literature data. We observed the temperature dissipation rate to vary by an order of magnitude over a horizontal transect at Ra > 10⁹. The temperature dissipation spectra obtained by both methods were in agreement over the Ra interval considered. The location of the temperature dissipation peaks was also consistent with PIV measured energy dissipation rates. Our data suggest increasing importance of top/bottom boundaries for the momentum and the temperature dissipation with increasing Ra number. Applied to oceanic upper ocean convection, our results imply that most of the dissipation occurs close to the air–sea boundary. Thus, attempts to parameterise or measure air–sea turbulent convective fluxes have to reflect the dominant role of near boundary dissipation at large Ra.     

Onset of Unsteady Horizontal Convection in Rectangular Tank at Pr=1

The horizontal convection within a rectangular tank is numerically simulated. The flow is found to be unsteady at high Rayleigh numbers. There is a Hopf bifurcation of Ra from steady solutions to periodic solutions, and the critical Rayleigh number Rac is obtained to be Rac = 5.5377×10^8 for the middle plume forcing at Pr = 1, which is much larger than the value previously obtained. In addition, the unstable perturbations are always generated from the central jet, which implies that the onset of instability is due to velocity shear （shear instability） other than thermally dynamics （thermal instability）. Finally, Paparella and Young＇s first hypotheses [J. Fluid Mech. 466 （2002） 205] about the destabilization of the flow is numerically proven, i.e. the middle plume forcing can lead to a destabilization of the flow.     

Laminar natural convection between vertical isothermal heated plates with different temperatures

Numerical investigation of laminar free convection heat transfer in the vertical parallel plate channel with asymmetric heating is presented. Both inlet and exit effects are included into the analysis. A numerical solution is obtained for a Prandtl number of 0.71 and for modified Rayleigh number [`(Ra)]\overline {Ra} = 10⁻¹ ÷ 10⁵, and varying heating ratio T_R = 0 ÷ 1 and aspect ratio A = 10. Fully elliptic Navier-Stokes and energy equations are solved using the finite volume techniques with staggered grid arrangements. The obtained results show a strong influence of the temperature ratio on local and average heat transfer coefficient on the hot and cold plates. With reduction of T_R the heat transfer parameter on the hot wall grows, and on the cold one, on the contrary, it decreases. As a result, the total heat exchange from two plates depends poorly on the parameter T_R.     

Determination of 226Ra, 232Th and 40K in teeth by use of gamma spectroscopy

This study entails the measurement of the specific activity of natural radionuclides (²²⁶Ra, ⁴⁰K and ²³²Th) in 18 tooth samples obtained from the clinic of the Universiti Sains Malaysia (USM), Penang, by using an HPGe detector. The specific activity of ²²⁶Ra, ⁴⁰K and ²³²Th was measured to estimate the hazard index of the radionuclides, radium equivalent activities (Ra_eq), external, internal hazard indices (H_ex, H_in), and absorbed dose (D_out, D_in). The maximum values of concentration of ²²⁶Ra, ²³²Th and ⁴⁰K in the tooth samples were found to be 60.82, 60.29 and 594.22?Bq kg^?1, respectively. Maximum values of Ra_eq, H_ex, H_in, D_out and D_in were found to be 192.78?Bq kg^?1, 0.520, 0.685, 89.29 and 169.81?nGy h^–1, Iγ and Iα as 0.702 and 0.304, respectively. The results were lower than the average world value (UNSCEAR). In addition, a strong correlation was found between the concentrations of ²²⁶Ra and Ra_eq, between energy and net area, as well as between radionuclides (²²⁶Ra, ⁴⁰K and ²³²Th) in tooth samples and age of volunteers. This study showed that the concentrations and hazard indices of tooth samples are below the recommended safe levels; therefore, the study area is considered safe in terms of radiological health hazards.     

Large-scale patterns in Rayleigh-Bénard convection

Rayleigh-Bénard convection at large Rayleigh number is characterized by the presence of intense, vertically moving plumes. Both laboratory and numerical experiments reveal that the rising and descending plumes aggregate into separate clusters so as to produce large-scale updrafts and downdrafts. The horizontal scales of the aggregates reported so far have been comparable to the horizontal extent of the containers, but it has not been clear whether that represents a limitation imposed by domain size. In this work, we present numerical simulations of convection at sufficiently large aspect ratio to ascertain whether there is an intrinsic saturation scale for the clustering process when that ratio is large enough. From a series of simulations of Rayleigh-Bénard convection with Rayleigh numbers between 10⁵ and 10⁸ and with aspect ratios up to 12π, we conclude that the clustering process has a finite horizontal saturation scale with at most a weak dependence on Rayleigh number in the range studied.     

Buoyancy induced flow in a nanofluid filled enclosure partially exposed to forced convection

A numerical study was performed on natural convection for water–CuO nanofluid filled enclosure where the top surface was partially exposed to convection. The cavity has a square cross-section and differentially heated. Except exposed convection part on the top, all sides are adiabatic on horizontal walls. Effects of Rayleigh number (10³ ? Ra ? 10⁵), Biot number (0 ? Bi ? ∞), length of partial convection (0.0 ? L ? 1.0) and volume fraction of nanoparticles (0.0 ? φ ? 0.1) on heat and fluid flow were investigated. The results showed that for the case of high Biot number that heat transfer along the heated was enhanced by increasing the Rayleigh number mainly at the upper portion of the heated wall. When the top wall was totally exposed to convection, the results prevail that the heat transfer was more effective at high Biot number especially at the upper portion of the heated wall. For the case of high Biot number, the results prevailed that the heat transfer at the upper portion of the heated wall increases considerably at high exposed length to convection (L); however, for L ? 0.75 the effect of L was less pronounced. Contour maps for percentage of heat transfer enhancement were presented and it was shown that the location of maximum enhancement in heat transfer was sensitive to Ra, φ and L.     

Computational approaches to aspect-ratio-dependent upper bounds and heat flux in porous medium convection

Direct numerical simulation (DNS) has shown that Rayleigh–Bénard convection in a fluid-saturated porous medium self-organizes into narrowly spaced plumes at (ostensibly) asymptotically high values of the Rayleigh number Ra. In this Letter a combination of DNS and upper bound theory is used to investigate the dependence of the Nusselt number Nu on the domain aspect ratio L at large Ra  . A novel algorithm is introduced to solve the optimization problems arising from the upper bound analysis, allowing for the best available bounds to be extended up to Ra≈2.65×10⁴$\mathit{Ra}\approx 2.65\times {10}^4$. The dependence of the bounds on L(Ra)$L(\mathit{Ra})$ is explored and a “minimal flow unit” is identified.