利用GPU实现基于物理模型的流体运动仿真

纳维-斯托克斯方程是描述流体运动的基础物理模型。首先简要介绍了在计算机上实现该模型的数值解算方法。在流体运动的计算机视觉仿真中,密度渲染才是最终结果,因此,又阐述了流体的密度场与速度场的关系,并且给出了流体仿真的基本流程。由于可编程图形处理器(GPU)具有并行性和速度快等特点,流体模拟的解算和渲染在图形处理器上运行更加高效,但是利用GPU编程实现流体仿真必须要实现3个核心转换,最后以支持3维纹理的Direct3D10为例,给出了利用GPU实现流体运动仿真两个简单实例。     

山区遥感图像自动分类研究——以西南地区为例(英文)

In mountainous areas, it is the undulant terrain, various types of geomorphic and land use that make the remote sensing images great metamorphism. Moreover, due to the elevation, there are many areas covered with shadow, clouds and snow that make the images more inaccurate. As a result, it would be very difficult to carry out auto-classification of RS images in these areas. The study took Southwest China as the case study area and the TM images, SPOT images as the basic information sources assisted by the a...     

A constrained LAMBDA method for GPS attitude determination

An improved method to obtain fixed integer ambiguity in GPS attitude determination is presented. Known conditions are utilized as constraints to acquire attitude information when the float solution and its variance–covariance matrix are not accurate enough. The searching ellipsoidal region is first expanded to compensate for errors caused by the inaccurate float solution. Then the constraints are used to shrink the region to a proper size, which maintains the true integer ambiguity. Experimental results demonstrate that this scheme gives a fast search time and a higher success rate in determining the fixed integer ambiguity than the unconstrained method. The accuracy of attitude angles is also improved.

A finite differences formulation for the linear and nonlinear dynamics of 2D catenary risers

A finite differences (FD) solution method is proposed for the numerical treatment of the dynamic equilibrium problem of 2D catenary risers. The method is based on the so-called Box approximation, which in the scope of the present contribution is applied to the complete nonlinear model as well as to the reduced linearized formulation. The application of the Box method transforms the original governing systems into convenient sets of algebraic equations, which in turn are solved efficiently by the relaxation method. Extensive numerical calculations are presented that describe the dynamic behaviour of the structure and evaluate the amplification in loading due to the dynamic components. The effect of the geometric nonlinearities is assessed through comparative calculations that concern both mathematical formulations examined in the present, i.e. the complete nonlinear, and the reduced linearized model. Special attention is paid to the heave excitations as they amplify significantly the magnitudes of the loading components.     

Assessment of boundary-element method for modelling a free-floating sloped wave energy device. Part 1: Numerical modelling

The boundary-element method has been widely used as a design tool in the offshore and ship building industry for more than 30 years. Its application to wave energy conversion is, however, more recent. This paper deals with the numerical modelling of a free-floating sloped wave energy device. The power take-off mechanism of the device consists of an immersed tube with a piston sliding inside. The modelling is done using the boundary-element method package WAMIT. The model is first worked out for the case where the axis of the tube is vertical. It is then derived for the tube inclined and successfully verified against numerical benchmark data. A companion paper presents results of a detailed comparison with a physical model study.     

Numerical calculation of free-surface potential flow around a ship using the modified Rankine source panel method

A modified Rankine source panel method is presented for solving a linearized free-surface flow problem with respect to the double-body potential. The method of solution is based on the distribution of Rankine sources on the hull as well as its image and on the free surface. An iterative algorithm is used for determining the free surface and wave resistance using upstream finite difference operator. A verification of numerical modeling is made using the Wigley hull and the validity of the computer program is examined by comparing the details of wave profiles and wave making resistance with Series 60 model.     

Verification of significant wave representation method

The significant wave representation method is the simplest method for computing the transformation of significant wave height across-shore. However, many engineers are reluctant to use this method because many researchers have pointed out that the method possibly contains a large estimation error. Nevertheless, Rattanapitikon et al. [Rattanapitikon, W., Karunchintadit, R., Shibayama, T., 2003. Irregular wave height transformation using representative wave approach. Coastal Engineering Journal, JSCE 45(3), 489–510.] showed that the wave representation method could be used to compute the transformation of root mean square wave heights. It may also be possible to use it for computing the significant wave height transformation. Therefore, this study was carried out to examine the possibility of simulating significant wave height transformation across-shore by using the significant wave representation method. Laboratory data from small- and large-scale wave flumes were used to calibrate and examine the models. Six regular wave models were applied directly to irregular waves by using the significant wave height and spectral peak period. The examination showed that three regular wave models (with new coefficients) could be used to compute the significant wave height transformation with very good accuracy. On the strength of both accuracy and simplicity of the three models, a suitable model is recommended for computing the significant wave height transformation. The suitable model was also modified for better predictions. The modified model (with different coefficients) can be used to compute either regular wave height or significant wave height transformation across-shore.     

Solution of the propagation of the waves in open channels by the transfer matrix method

Many problems in mechanics can be solved by the use of the transfer matrix method. The use of this method in hydraulics engineering is not widespread and only limited studies are available. In this study, linearized St. Venant equations were used and the use of transfer matrix in ocean engineering was investigated for long waves in open channels, and numerical application was carried out. The results obtained through the transfer matrix method, which is quite easy to use, program and comprehend, showed similar results obtained from the characteristics method and finite differences method.     

极限波浪运动特性的非线性数值模拟

利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究.     

Microbial photosynthesis in coral reef sediments (Heron Reef, Australia)

We investigated microphytobenthic photosynthesis at four stations in the coral reef sediments at Heron Reef, Australia. The microphytobenthos was dominated by diatoms, dinoflagellates and cyanobacteria, as indicated by biomarker pigment analysis. Conspicuous algae firmly attached to the sand grains (ca. 100 μm in diameter, surrounded by a hard transparent wall) were rich in peridinin, a marker pigment for dinoflagellates, but also showed a high diversity based on cyanobacterial 16S rDNA gene sequence analysis. Specimens of these algae that were buried below the photic zone exhibited an unexpected stimulation of respiration by light, resulting in an increase of local oxygen concentrations upon darkening. Net photosynthesis of the sediments varied between 1.9 and 8.5 mmol O₂ m⁻² h⁻¹ and was strongly correlated with Chl a content, which lay between 31 and 84 mg m⁻². An estimate based on our spatially limited dataset indicates that the microphytobenthic production for the entire reef is in the order of magnitude of the production estimated for corals. Photosynthesis stimulated calcification at all investigated sites (0.2–1.0 mmol Ca²⁺ m⁻² h⁻¹). The sediments of at least three stations were net calcifying. Sedimentary N₂-fixation rates (measured by acetylene reduction assays at two sites) ranged between 0.9 to 3.9 mmol N₂ m⁻² h⁻¹ and were highest in the light, indicating the importance of heterocystous cyanobacteria. In coral fingers no N₂-fixation was measurable, which stresses the importance of the sediment compartment for reef nitrogen cycling.