山地边界层中场变量的动力学研究

采用与实测较接近的二次函数来表达Ekman层中的湍流粘性系数K,在圆形气压场条件下,求得了山地上空边界层中的风速,进而求得散度、涡度和垂直速度等场变量随高度的分布。并作图分析了这些场变量的一些动力学特征。改进了以往在求解析解时,略去运动方程中湍流粘性力项中的关于高度的一阶导数项,以及取山坡面上风速为零作下边界条件等欠合理欠精确的做法。所求得的风速、散度、涡度和垂直速度均用简单的初等函数表示出来,有助于边界层参数化和深化对边界层动力学的认识。     

不同岩性组合下的斜坡地震动力响应研究

本文采用大型有限差分软件FLAC3D, 建立上硬下软, 上软下硬, 含软弱基座和不含软弱基座4种斜坡的数值分析模型, 输入阪神地震波, 分析总结斜坡的动力响应规律。计算结果表明:(1)从加速度动力响应规律来看, 无论是水平向还是竖直向, 上软斜坡均比上硬斜坡大。并且, 两种斜坡的水平向放大系数均比竖直向放大系数大, 基本在1.52.0倍之间。(2)从位移规律上来看, 软岩相比于硬岩, 其对地震波更加敏感, 上硬下软岩性组合在强震动力作用下更易失稳。(3)从速度、位移的动力响应规律来看, 有软基斜坡地震动响应程度明显大于无软基斜坡。但是从边坡的应力变化来看, 软弱基座的存在, 一定程度上起到了吸收入射地震波能量的作用, 反而有利于边坡稳定性。     

Solving the linearized forward-speed radiation problem using a high-order finite difference method on overlapping grids

The linearized potential flow approximation for the forward speed radiation problem is solved in the time domain using a high-order finite difference method. The finite-difference discretization is developed on overlapping, curvilinear body-fitted grids. To ensure numerical stability, the convective derivatives in the free-surface boundary conditions are treated using an upwind-biased stencil. Instead of solving for the radiation impulse response functions, a pseudo-impulsive Gaussian type displacement is employed in order to tailor the frequency-content to the discrete spatial resolution. Frequency-domain results are then obtained from a Fourier transform of the force and motion signals. In order to make a robust Fourier transform, and capture the response around the critical frequency, the tail of the force signal is asymptotically extrapolated assuming a linear decay rate. Fourth-order convergence of the calculations on simple geometries is demonstrated, along with a nearly linear scaling of the solution effort with increasing grid resolution. The code is validated by comparison with analytical and semi-analytical solutions using submerged and floating closed-form geometries. Calculations are also made for a modern bulk carrier, and good agreement is found with experimental measurements.     

洞庭湖泥沙淤积数值模拟模式

在确保沙量守恒的非恒定流、非均匀沙的二维水沙数学模型的基础上,通过洞庭湖水沙输运和河床变形计算实践,比较了两种冲泻质与床沙质的转化模式,提出了包括冲泻质在内的泥沙淤积模式和水沙动边界计算模式.从洞庭湖泥沙输运和河床变形计算的结果看,洞庭湖泥沙淤积数值模式具有较好的稳定性和较高的模拟精度,说明所建湖泊泥沙模型的良好性能和具有推广运用的前景.     

An Effective Control System of Thrust and Moment Modulation for a Dynamically Positioned Vessel

Dynamic positioning (DP) is an operation method whereby the position of a surface vessel is maintained in close proximity to a required position in the horizontal plane through the controlled application of forces and moments generated by purposely installed thrusters. When estimating thrust, this kind of conventional control system often uses many acceleration sensors, velocity sensors, environment sensors, and filters. Usually, these sensors have measured electrical errors. To reduce the number of sensors used and to decrease the measurement errors, this article presents an effective control system for estimating thrust and moment commands, which is based on energy and impulsemomentum principles. Donha and Brinati's example is followed to verify the feasibility of the present control system, which performs semisubmersible platform positioning using an LQG controller, and the results are feasible and economical. A simulated coring vessel marine positioning in southern Taiwan is presented, which can estimate the counterthrust and moment commands, and the complex environmental forces and moments are described. The results can provide a valuable control system for dynamically positioned vessels.     

On dynamic coupling effects between a spar and its mooring lines

The responses of a spar constrained by slack mooring lines to steep ocean waves and tensions in the mooring lines are simulated using two different numerical schemes: a quasi-static approach (SMACOS) and a coupled dynamic approach (COUPLE). The two approaches are the same in computing wave loads on the structure. Their difference is in modeling dynamic forces of mooring lines; that is the dynamic forces are included in the computation of COUPLE but neglected in SMACOS. The numerical simulation is examined against the laboratory measurements of the JIP Spar in a water depth of 318 m. The dynamic coupling effects between the JIP Spar and its mooring lines in different water depths (318, 618 and 1018 m) are investigated by the comparison of numerical simulations obtained using the quasi-static and coupled dynamic approaches. It is found that the damping of mooring lines reduces the slow-drift surge and pitch of the Spar, especially in deep water. The reduction in the amplitude of slow-drift surge can reach about 10% in a water depth of 1018 m. The tension in mooring lines may greatly increase in the wave frequency range when dynamic forces in mooring lines are considered. The mooring-line tension in the wave frequency range predicted by the coupled dynamic approach can be eight times as great as the corresponding prediction by the quasi-static approach in a water depth of 1018 m. This finding may have important implications for the estimation of the fatigue strength and life span of the mooring lines deployed in deep water oceans.     

A parallel dynamic programming algorithm for multi-reservoir system optimization

This paper develops a parallel dynamic programming algorithm to optimize the joint operation of a multi-reservoir system. First, a multi-dimensional dynamic programming (DP) model is formulated for a multi-reservoir system. Second, the DP algorithm is parallelized using a peer-to-peer parallel paradigm. The parallelization is based on the distributed memory architecture and the message passing interface (MPI) protocol. We consider both the distributed computing and distributed computer memory in the parallelization. The parallel paradigm aims at reducing the computation time as well as alleviating the computer memory requirement associated with running a multi-dimensional DP model. Next, we test the parallel DP algorithm on the classic, benchmark four-reservoir problem on a high-performance computing (HPC) system with up to 350 cores. Results indicate that the parallel DP algorithm exhibits good performance in parallel efficiency; the parallel DP algorithm is scalable and will not be restricted by the number of cores. Finally, the parallel DP algorithm is applied to a real-world, five-reservoir system in China. The results demonstrate the parallel efficiency and practical utility of the proposed methodology.     

Comparison of numerical solution and semi-empirical formulas to predict the effects of important design parameters on porpoising region of a planing vessel

Stability of the marine vessels in different conditions is one of the most important problems in the design of a planing vessel. In this research, the effects of some important design parameters (mass, longitudinal center of mass, deadrise angle, and length) of DTMB 62 model 4667-1 planing hull on the drag and also on the longitudinal dynamic stability (porpoising) are investigated in the velocity range of 2.12–8.486 m/s in calm water. In this paper, both numerical simulation of Reynolds Average Naiver Stokes (RANS) equations and semi-empirical formulas of Savitsky are used to analyze the motion of a 4667 planing vessel in calm water with two degrees of freedom (2DOF). For this purpose a finite volume, ANSYS-FLUENT, code is used to solve the Navier-Stokes equations for the simulation of the flow field around the vessel. In addition, an explicit VOF scheme and SST- Kω model is used with dynamic mesh scheme to capture the interface of a two-phase flow and to model the turbulence respectively, in 2DOF model (heave and pitch). Also, the results of both methods are compared with each other. According to the present results, changing the aspect ratio of the vessel and also the longitudinal center of gravity have the most effect on the porpoising region.     

Fragility and risk assessment of freestanding building contents

Multistorey buildings often have a valuable inventory consisting of objects that their possible damage during an earthquake will cause unacceptable losses. The paper presents a novel, fully performance-based seismic reliability and risk assessment framework for freestanding structural components and contents that can be modelled as rocking rigid blocks. The seismic response of building contents depends on several parameters such as the geometry of the object, the dynamic characteristics of the building and the storey that the object is located. The demand at the storey level is first obtained, and then the response of the contents is calculated using the storey acceleration response history. The demand of the structure is obtained with the aid of a modified version of the Incremental Dynamic Analysis method and subsequently the fragility curves of the rocking building contents are derived for every storey of interest. Different options for fragility assessment are discussed, and the underlying details of the problem are investigated. A simplified approach, where the fragility of the freestanding components and the structure are derived separately, is also presented. The method combines existing fragility curves and thus is suitable for quickly assessing the reliability of a building's inventory, offering sufficient risk estimates.     

Stream-coordinate structure of oceanic jets based on merged altimeter data

The jet structure of the Southern Ocean front south of Australia is studied in stream-coordinate with a new altimeter product—Absolute Dynamic Topography (ADT) from AVISO. The accuracy of the ADT data is validated with the mooring data from a two-year subantarctic-front experiment. It is demonstrated that the ADT is consistent with in-situ measurements and captures the meso-scale activity of the Antarctic Circumpolar Current (ACC). Stream-coordinate analysis of ADT surface geostrophic flows finds that ACC jets exhibit large spatio-temporal variability and do not correspond to particular streamfunction values. In the circumpolar scope ACC jets display a transient fragmented pattern controlled by topographic features. The poleward shift of jet in streamfunction space, as revealed by a streamwise correlation method, indicates the presence of meridional fluxes of zonal momentum. Such cross-stream eddy fluxes concentrate the broad ACC baroclinic flow into narrow jets. Combined with a recent discovery of gravest empirical mode (GEM) in the thermohaline fields, the study clarifies the interrelationship among front, jet and streamfunction in the Southern Ocean.