喷灌均匀系数的三次样条两次插值计算方法

为提高喷灌水量分布均匀性评价的准确性,当雨量筒径向布置时,为考虑所有测点数据对插值点降水深的影响,采用径向和周向两次的三次样条插值计算出未知点的降水深,从而计算喷灌均匀系数。以美国雨鸟30PSH型喷头雨量筒间隔为1m和2m的喷洒试验数据,计算网格点取1m和0.25m,分别采用三次样条两次插值法和邻近四点距离线性插值法计算了克里斯琴森均匀系数。结果表明,均匀系数由高至低的顺序依次为采样间隔为2m的线性插值、采样间隔为2m的三次样条两次插值、采样间隔为1m的线性插值和采样间隔为1m的三次样条两次插值。采样间隔2m比1m计算出的均匀系数总体高3~4个百分点,三次样条两次插值法比邻近点距离线性插值法略低1个百分点,2种计算网格点间距下的均匀系数差值小于1个百分点。结果证明,采样间距、插值方法、计算网格间距对均匀系数的影响依次降低,三次样条两次插值法可以用来评价喷灌组合均匀系数。

AIN Monolithic Microchannel Cooled Heatsink for High Power Laser Diode Array

A novel AIN monolithic microchannel cooled heatsink for high power laser diode array is introduced.The high power stack laser diode array with an AIN monolithic microchannel heatsink is fabricated and tested.The thermal impedance of a 10 stack laser diode array is 0.121℃/W.The pitch between two adjacent bars is 1.17mm.The power level of 611W is achieved under the 20% duty factor condition at an emission wavelength around 808nm.

Design of spectrometer based on volume phase grating for near infrared range

In order to design and fabricate a spectrometer for the infrared range widely used in the different applications, Volume Phase Grating (VPG) with. low Polarization Dependence Loss (PDL) and high efficiency has been adopted as the dispersion element. VPG is constructed by coating an optical substrate with a thin film of dichromated. gelatin and exposing the film to two mutually coherent laser beams to form index modulation. The diffraction efficiency for a VPG is governed by Bragg effects. The depth (d) and index modulation contrast of the grating structure control the efficiency at which the light is diffracted when the Bragg condition is satisfied. Gradient index lens with high performance and low aberration are used as collimating system instead of standard lens. The spot diagrams and MTF curve of the collimating lens are shown in the paper. The receive system is InCaAs photodiode (PD) array including 512 pixels with 25 mum pitch. The spectrum resolution of the spectrometer reaches to 0.2nm and wavelength accuracy is 40pm.

Mode theory and analysis of planar array waveguides

Bloch modes can be excited in planar array due to its periodic lateral refractive index. The power coupled into each eigenmode of the array waveguides is calculated through the overlap integrals of the input field with the eigenmode fields of the coupled infinite array waveguides projected onto the x-axis. Low losses can be obtained if the transition from the array to the free propagation region is adiabatic. Due to the finite resolution of lithographic process the gap between the waveguides will stop abruptly, however, when the waveguides come into too close together. Calculation results show that losses will occur at this discontinuity, which are dependent on the ratio of the gap between the waveguides and grating pitch and on the confinement of field in the array waveguides. Tapered waveguides and low index contrast between the core and cladding layers can lower the transmitted losses.

First charge collection and position-precision data on the medium-resistivity silicon strip detectors before and after neutron irradiation up to 2x10(14) n/cm(2)

Test strip detectors of 125 mu m, 500 mu m, and 1 mm pitches with about 1 cm(2) areas have been made on medium-resistivity silicon wafers (1.3 and 2.7 k Ohm cm). Detectors of 500 mu m pitch have been tested for charge collection and position precision before and after neutron irradiation (up to 2 x 10(14) n/cm(2)) using 820 and 1030 nm laser lights with different beam-spot sizes. It has been found that for a bias of 250 V a strip detector made of 1.3 k Ohm cm (300 mu m thick) can be fully depleted before and after an irradiation of 2 x 10(14) n/cm(2). For a 500 mu m pitch strip detector made of 2.7 k Ohm cm tested with an 1030 nm laser light with 200 mu m spot size, the position reconstruction error is about 14 mu m before irradiation, and 17 mu m after about 1.7 x 10(13) n/cm(2) irradiation. We demonstrated in this work that medium resistivity silicon strip detectors can work just as well as the traditional high-resistivity ones, but with higher radiation tolerance. We also tested charge sharing and position reconstruction using a 1030 nm wavelength (300 mu m absorption length in Si at RT) laser, which provides a simulation of MIP particles in high-physics experiments in terms of charge collection and position reconstruction, (C) 1999 Elsevier Science B.V. All rights reserved.

地形可视性分析中观察点设置问题研究

观察点设置问题是地形可视性分析中的一类重要问题，对该问题的研究可以在空间信息辅助决策、通信、旅游、野生动物保护等领域发挥重大作用。本文在对地形可视性分析中观察点设置问题现有研究成果总结和分析基础上对该问题展开深入研究。 首先，针对现有解决方法只从智能算法或地形数据表示方法单一角度进行分析和研究的局限性，提出了一种问题相关的智能算法和数据表示方法相结合的解决问题新框架。该框架考虑了解决观察点设置问题时智能算法的优点和数据表示方式的特点相互配合问题，目的是充分发挥二者各自的优势以提高观察点设置问题解决的准确度与效率。 其次，在深入分析观察点设置问题本身特点的基础上，结合隶属云理论的基本理论和方法，对经典模拟退火算法从退温函数设计、温度产生过程、状态生成过程三方面进行了问题相关的改进，提出了一种适于观察点设置问题的改进模拟退火算法(Improved Simulated Annealing algorithm, ISA)。该算法一方面保持了经典模拟退火算法的稳定倾向特性，保证了算法满足伴随退火温度的不断下降，对恶化的新状态越来越难于接受这一模拟退火算法的最基本特征；另一方面其退火温度的连续性随机变化特性和隐含的“回火升温”过程，则有利于算法有效拒绝恶化解，加速算法收敛，能够更好地满足观察点设置问题对于算法收敛速度的要求。 再次，在分析地形数据的精度、误差等因素对观察点设置问题的解决准确性和解决效率影响程度的基础上，提出了一种基于离散余弦变换的地形数据内插方法(Discrete Cosine Transformation Interpolation method, DCTI)。新方法将传统空域上的地形内插转换到变换域上进行，同时充分利用了离散余弦变换的熵保持特性和能量压缩特性，简化了变换域上的内插过程，提高了地形数据内插的效率和精度。DCTI方法与其他现有典型地形数据内插方法相比，对地形可视性信息获取的准确性和效率影响最小，为平衡观察点设置问题解决过程中时间效率和准确度之间的关系，最终有效地解决观察点设置问题提供了数据基础。 最后，从智能算法和地形数据相结合的角度出发，提出了一种基于ISA和DCTI相结合的观察点设置问题多分辨率处理方法(Multi-Resolution Processing method, MRP)。新方法将模拟退火算法的逐次退火特点和地形数据的多分辨率表示充分结合，达到了发挥算法数据相结合的综合优势的目的。与现有单纯基于模拟退火算法的解决方法相比，在问题解决准确度保持不变的前提下，基于MRP方法的观察点设置问题解决的平均耗时减少85%~95%，为实际工程应用问题的解决提供了一条重要途径。

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.