一种线性走时插值射线追踪改进算法

针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍.      

Sensitivity of simulated flow fields and bathymetries in meandering channels to the choice of a morphodynamic model

Morphodynamic models are used by river practitioners and scientists to simulate geomorphic change in natural and artificial river channels. It has long been recognized that these models are sensitive to the choice of parameter values, and proper calibration is now common practice. This paper investigates the less recognized impact of the choice of the model itself. All morphodynamic models purport to simulate the same flow and sediment dynamics, often relying on the same governing equations. Yet in solving these equations, the models have different underlying assumptions, for example regarding spatial discretization, turbulence, sediment inflow, lateral friction, and bed load transport. These differences are not always considered by the average model user, who might expect similar predictions from calibrated models. Here, a series of numerical simulations in meandering channels was undertaken to test whether six morphodynamic codes (BASEMENT, CCHE‐2D, NAYS, SSIIM‐1, TELEMAC‐2D and TELEMAC‐3D) would yield significantly different equilibrium bathymetries if subjected to identical, initial flow conditions. We found that, despite producing moderately similar velocity patterns on a fixed‐flat bed (regression coefficient r of 0.77 ± 0.20), the codes disagree substantially with respect to simulated bathymetries (r = 0.49 ± 0.31). We relate these discrepancies to differences in the codes' assumptions. Results were configuration specific, i.e. codes that perform well for a given channel configuration do not necessarily perform well with higher or lower sinuosity configurations. Finally, limited correlation is found between accuracy and code complexity; the inclusion of algorithms that explicitly account for the effects of local bed slope and channel curvature effects on transport magnitude and direction does not guarantee accuracy. The range of solutions obtained from the evaluated codes emphasizes the need for carefully considering the choice of code. We recommend the creation of a central repository providing universal validation cases and documentation of recognized fluvial codes in commonly studied fluvial settings. Copyright © 2016 John Wiley & Sons, Ltd.     

Jet flow over foredunes

Jet flows, which are localized flows exhibiting a high speed maxima, are relatively common in nature, and in many devices. They have only been occasionally observed on dunes, and their dynamics are poorly known. This paper examines computational fluid dynamic (CFD) two‐dimensional (2D) modelling of jet flow over a foredune topography. Flow was simulated in 10° increments from onshore (0°) to highly oblique alongshore (70°) incident wind approach angles. CFD modelling reveals that the formation of a jet is not dependent on a critical wind speed, and an increase in incident wind velocity does not affect the magnitude of jet flow. A jet is first formed at ~1.0 m seawards of the foredune crest on the Prince Edward Island foredune morphology example examined here. A jet is not developed when the incident wind is from an oblique approach angle greater than ~50° because there is significantly less flow acceleration across a much lower slope at this incident angle. The presence of a scarp does influence the structure of the crest jet, in that the jet is more pronounced where a scarp is present. Surface roughness affects the magnitude of jet expansion and jets are better developed on bare surfaces compared to vegetated ones. Copyright © 2016 John Wiley & Sons, Ltd.     

Modal model identification of structures under unmeasured seismic excitations

Techniques developed for structural identification of a structural model are typically based on information regarding the response and the forcing actions. However, in some situations it can be necessary, or simply useful, to refer only to the measured responses. In this paper we describe a technique suitable for identifying the modal model of a spatial frame in the frequency domain when the seismic input is unknown both in time contents and direction. In some previous theoretical works we established that this identification problem has a unique solution when at least three time‐history responses are known. Here numerical techniques are developed which allow the evaluation of the modal quantities in practice. Numerical applications are carried out on plane and spatial framed structures by using a modal model which may be complete, including all the structure's modes, or incomplete, including only the lowest modes. In most cases the obtained results are satisfactory. Copyright © 2005 John Wiley & Sons, Ltd.     

A bidirectional path-finding algorithm and data structure for maritime routing

Route planning is an important problem for many real-time applications in open and complex environments. The maritime domain is a relevant example of such environments where dynamic phenomena and navigation constraints generate difficult route finding problems. This paper develops a spatial data structure that supports the search for an optimal route between two locations while minimizing a cost function. Although various search algorithms have been proposed so far (e.g. breadth-first search, bidirectional breadth-first search, Dijkstra’s algorithm, A^*, etc.), this approach provides a bidirectional dynamic routing algorithm which is based on hexagonal meshes and an iterative deepening A^* (IDA^*) algorithm, and a front to front strategy using a dynamic graph that facilitates data accessibility. The whole approach is applied to the context of maritime navigation, taking into account navigation hazards and restricted areas. The algorithm developed searches for optimal routes while minimizing distance and computational time.     

用于气象站的地表气温观测仪器的设计与实验研究北大核心

大气科学研究对地表气温观测精度有高达0.1℃甚至0.05℃的需求。然而,现有的地表气温观测仪器受到太阳直接辐射、下垫面反射辐射、长波辐射和散射辐射等影响,辐射误差可达1℃。本文设计了一种基于导流装置的地表气温观测仪器。首先,利用计算流体动力学(Computational Fluid Dynamics,CFD)方法量化该仪器在各种环境条件下的辐射误差;然后,在此基础上,利用极限学习机(Extreme Learning Machine,ELM)方法拟合可针对多变量变化的辐射误差订正方程;最后,为验证该仪器的观测精度,进行了外场比对实验。在实验过程中,以076B型强制通风气温观测仪器的测量值作为温度基准。实验结果表明,该仪器的平均辐射误差和最大辐射误差分别为0.07℃和0.15℃。该仪器辐射误差的实验测量值与订正方程提供的辐射误差订正值之间的平均偏移量、均方根误差和相关系数分别为0.033℃、0.028℃和0.703。