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

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

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.     

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

礁坪上波浪传播的数值模拟

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。     

Nonlinear dynamo effects for an inhomogeneously turbulent rotating fluid

Im Rahmen der Magnetohydrodynamik der mittleren Felder werden die turbulenzbedingte elektromotorische Kraft und der turbulenzbedingte Spannungstensor für eine leitende Flüssigkeit mit inhomogener Turbulenz auf einem rotierenden Körper berechnet. Dabei wird der Einflu$sZ des Magnetfeldes auf der Turbulenz berücksichtigt. So erscheinen nichtlineare Effekte bei den Vorgängen, die für das Vorhandensein von Magnetfeld und differentieller Rotation verantwortlich sind; die dafür ma$sZgebenden Koeffizienten sind Funktionen des Magnetfeldes. Die Nichtlinearitäten werden bis zur zweiten Ordnung im mittleren Magnetfeld (d. h. als schwache Nichtlinearitäten) erfa$sZt. Der $aL-Effekt beim $aL$OM-Dynamo erleidet sowohl bei langsamer als auch bei rascher Rotation eine Schwächung durch das Magnetfeld. Es werden einige astrophysikalische Anwendungen der Ergebnisse erörtert. Insbesondere wird auf die Möglichkeit hingewiesen, da$sZ die bei der Sonne beobachteten Torsionsschwingungen als Folge des im Aktivitätszyklus schwankenden magnetischen Einflusses auf den für die differentielle Rotation ma$sZgebenden A-Effekt erscheinen.     

Diel and seasonal variability in the fish community structure of a mud‐bottom estuarine habitat in the Gulf of Mexico

We hypothesized that temporal variation in fish species composition and community structure in a low complexity habitat in the Pueblo Viejo Lagoon, Mexico, is influenced by diel light/dark cycles and tidal stage, and by seasonal changes in salinity and temperature. We collected a total of 17,661 individuals during 2‐h interval sampling over six bi‐monthly 24‐h sampling cycles representing 53 species, of which 11 (～20%) were previously unknown in the system. Diel variation indicated that significantly higher numbers of individuals and species were caught at night, whereas diversity and evenness were higher during the day. Species richness was significantly higher in July and January, whereas diversity and evenness peaked around May; both were correlated with temperature. Diel variation in species composition was influenced primarily by the light/dark cycle. Cluster analyses of each diel cycle separated fish assemblages from midday samples from those of nocturnal samples, separated by an extended wide transition period as fish moved at dawn and during the late afternoon/dusk. Significant shifts (as determined by MANOVA) in assemblage structure occurred between months. Canonical correspondence analysis showed that temperature and day/night effects were the most important environmental variables structuring the fish community. This constrained ordination also defined species with specific habitat preferences as follows: (i) diurnal, warm temperature species (mainly planktivores) (Brevoortia gunteri, Cetengraulis edentulus, Diapterus auratus, and Membras martinica); (ii) nocturnal, warm temperature species (mainly predators) (Citharichthys spilopterus, Cathorops melanopus, and Bairdiella spp.); and (iii) low temperature, diurnal species (Brevoortia patronus and Mugil curema) or those with twilight and nocturnal distributions (Anchoa mitchilli, the most numerically abundant species). Our results indicate that diel and seasonal changes in fish community structure were mainly related to day/night cycles and temperature regimes.     

微扰Kepler系统轨道微分方程的近似Lie对称性与近似不变量

把极角θ视为独立变量,得到Kepler系统的轨道微分方程. 首先讨论Kepler系统轨道微分方程的Lie对称性和不变量,微扰Kepler系统轨道微分方程的精确Lie对称性和精确不变量,其次讨论微扰Kepler系统轨道微分方程的近似Lie对称性和近似不变量,并得到了微扰Kepler系统的9个一阶近似Lie对称性和6个一阶近似不变量,其中1个实为精确不变量,而其余5个分别等于微扰系数ε乘以Kepler系统相应的5个不变量.     

New finite-gap solutions for the coupled Burgers equations engendered by the Neumann systems

On the tangent bundle TS^N-1 of the unit sphere S^N-1, this paper reduces the coupled Burgers equations to two Neumann systems by using the nonlinearization of the Lax pair, whose Liouville integrability is displayed in the scheme of the r-matrix technique. Based on the Lax matrix of the Neumann systems, the Abel--Jacobi coordinates are appropriately chosen to straighten out the restricted Neumann flows on the complex torus, from which the new finite-gap solutions expressed by Riemann theta functions for the coupled Burgers equations are given in view of the Jacobi inversion.     

雨滴谱数据在S波段双偏振天气雷达数据质量评估中的应用

结合T-Matrix方法和雨滴谱仪数据反演双偏振参量,设计了一种雨滴谱反演评估方法(雨滴谱法),利用其反演结果可对双偏振天气雷达数据进行评估。将此方法应用于不同类型降水过程的雷达偏振参量数据质量评估,结果表明:满足微雨条件时,雨滴谱法与传统的微雨滴法均能有效监测由雷达系统误差导致的观测偏差,且评估结果十分接近;不满足微雨条件时,雨滴谱法评估的不同类型降水过程的雷达数据质量与满足微雨条件下的评估结果十分接近,客观证明了系统调整后雷达运行良好、数据质量稳定;不同雨强(R)情况下,雨滴谱法的评估结果具有一致性。10 mm·h^-1≤R<20 mm·h^-1时,雷达观测的数据质量最为稳定。该雨滴谱法适用范围广,可动态评估双偏振雷达偏振参量质量,有助于雷达调整和数据校正。