Improvement on open boundaries on a time dependent numerical model of wave propagation in the nearshore region

An improvement on the simulation of outgoing waves on a time dependent numerical model for water wave propagation in the nearshore region is presented. The governing equations consist of a system of first order partial differential equations (PDEs), the equation of continuity and the equation of motion. A comparative study of first order radiation boundary conditions (BCs) and first order radiation BCs combined with sponge layers is presented for cases where outgoing waves leave the numerical domain of calculation through the open boundary. A reduction of spurious reflections from the numerical open boundaries can be obtained with an irrelevant increase in terms of computational cost.     

Third order contribution to the wave-making resistance of a ship at finite depth of water

This paper presents a potential based boundary element method for solving a nonlinear free surface flow problem for a ship moving with a uniform speed in finite depth of water. The free surface boundary condition is linearized by the systematic method of perturbation in terms of a small parameter up to third order. The surfaces are discretized into flat quadrilateral elements and the influence coefficients are calculated by Morino's analytical formula. Dawson's upstream finite difference operator is used in order to satisfy the radiation condition. The second order solution gives better result than the first or third order solution. So the present method with the second order solution can be adopted as a powerful tool for the hydrodynamic analysis of the thin ship in finite depth of water.     

On Radiation Boundary Conditions and Wave Transformation Across the Surf Zone

The purpose of this paper is to extend the validity of Li's parabolic model (1994) by incorporating a combined energy factor in the mild-slope equation and by improving the traditional radiation boundary conditions. With wave breaking and energy dissipation expressed in a direct form in the equation, the proposed model could provide an efficient numerical scheme and accurate predictions of wave transformation across the surf zone. The radiation boundary conditions are iterated in the model without use of approximations. The numerical predictions for wave height distributions across the surf zone are compared with experimental data over typical beach profiles. In addition, tests of waves scattering around a circular pile show that the proposed model could also provide reasonable improvement on the radiation boundary conditions for large incident angles of waves.     

Second-order Taylor expansion boundary element method for the second-order wave radiation problem

A novel Boundary Element Method (BEM) named the second-order Taylor Expansion Boundary Element Method (the 2nd order TEBEM) is developed for the solution of the second-order wave radiation velocity potential and sum-frequency wave loads for floating bodies. The radiation condition is enforced by a hybrid method of the multi-transmitting formula and damping zone. For the interior domain problem of a cube and a sphere, numerical results demonstrate that the 2nd order TEBEM can accurately solve the first and second-order gradients of velocity potential on the no-smoothed and smoothed boundary compared to the low-order BEM. The double frequency forces acting on the truncated cylinder are calculated under finite water depth. The agreement between the 2nd order TEBEM and others' numerical results is good. Moreover, all of the singular integrals in the 2nd order TEBEM can be solved analytically, so its implementation is much easier compared to the high-order BEM.     

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.     

The prediction of the dynamic and structural motions of a floating-pier system in waves

In this study, a two-dimensional floating pier consists of single rectangular impermeable pontoon with side supporting pile-columns is studied. The purpose of this study is to present a theoretical solution for the linearized problem of incident waves exerting on a floating pier with pile-restrained. All boundary conditions are linearized in the problem, which is incorporated into a scattering problem and radiation problem with unit displacement. The method of separation of variables is used to solve for velocity potentials. For the radiation problem with unit heave and pitch amplitude, the boundary value problem with non-homogeneous boundary condition beneath the structure is solved by using a solution scheme. By calculating the wave force from velocity potential and solving the equation of motion of the floating structure simultaneously a close form theoretical solution for the problem is developed. The finite element method was also applied to calculate the dynamic responses on the supporting piles subjected to the pontoon motions and incident waves.     

Hydrodynamics of a submerged hydrofoil advancing in waves

The hydrodynamic problem of a hydrofoil travelling at constant speed in water waves has been investigated through velocity potential theory. The boundary conditions on the free surface have been linearized, and the effects are accounted for through the Green function. The overall problem is decomposed into the steady forward speed problem and periodic wave radiation and diffraction problems. Each of these problems is solved using the boundary integral equation over the hydrofoil surface together with a vortex sheet behind the trailing edge. The body surface boundary condition is imposed on its mean position. As a result the steady potential will contribute a well-known m_j term to the body surface boundary condition on the radiation problem. The numerical difficulty in dealing with this term is effectively resolved through a difference method. The effects of the thickness on the wave radiation and diffraction are investigated. The applicability of various reciprocity relationships in this problem is discussed.     

Modeling the water circulation in the North Atlantic in the scope of the CORE-II experiment

A numerical experiment on the reproduction of the variability in the state of North Atlantic water in 1948–2007 with a spatial resolution of 0.25° has been performed using the global ocean model developed at Institute of Numerical Mathematics, Russian Academy of Sciences (INM RAS), and the Shirshov Institute of Oceanology (IO RAS) (the INM–IO model). The data on the state of the atmosphere, radiation fluxes, and bulk formulas of the CORE-II protocol are used as boundary conditions. Five successive 60-year calculation cycles have been performed in order to obtain the quasi-equilibrium state of a model ocean. For the last 20 years, the main elements of large-scale ocean circulation have been analyzed and compared with the WOA09 atlas data and the results of other models.     

Open boundaries in short wave simulations — A new approach

A new way of implementing radiation boundary conditions in finite difference schemes is reported. Instead of prescribing the incident field at the model boundary, waves are generated inside the model boundary. All outgoing waves are absorbed at open boundaries using so-called “sponge” layers.     

Estimation of underwater scattering radiances up to the third order

Vertical distribution of underwater scattering radiances of the first, the second and the third order is estimated by the Jerlov's method. In the estimation the volume attenuation coefficient and the volume scattering function are assumed to be constant in near surface layer, and the sky radiation separated from solar radiation to be a parallel beam and to fall into the sea in the vertical direction only.In the zenith direction the second scattering radiance is more intense than the first scattering radiance in near surface layer. As for the total radiance which is a sum of radiances up to the third order, the zenith radiation is largest, the nadir radiation is moderate and the horizontal radiation is least.     

Vibration analysis of rectangular Mindlin plates on elastic foundations and vertically in contact with stationary fluid by the Ritz method

In this study Free vibration analysis of vertical rectangular Mindlin plates resting on Pasternak elastic foundation and fully or partially in contact with fluid on their one side is investigated for different combinations of boundary conditions. The plate is assumed to be one of vertical rectangular walls of a container in contact with fluid. In order to analyze the interaction of the Mindlin plate with the elastic foundation and fluid system, three displacement components of the plate are expressed in the Ritz method by adopting a set of static Timoshenko beam functions satisfying geometric boundary conditions in a Cartesian co-ordinate system. The method of separation of variables and the method of Fourier series expansion is used to model fluid and to obtain the exact expression of the motion of fluid in the form of integral equations. The fluid domain is finite in depth and width but infinite in the length direction. To demonstrate the accuracy of the present solution, convergence study is first carried out and then a few comparison studies are carried out with the available data in the literature. Finally, natural frequencies of rectangular plates are presented in tabular and graphical forms for different fluid levels, foundation parameters, aspect ratios, thickness to width ratios and boundary conditions.     

东山湾波浪对悬沙浓度场影响的数值模拟研究

研究了波致底切应力和二维辐射应力对悬沙的作用.首先对东山湾的水动力进行了数值模拟并与实测资料进行了对比检验,然后对东山湾悬浮泥沙进行了考虑波浪和不考虑波浪两种情况下的数值模拟并与实测资料进行了相关对比分析.在模式建立过程中,依据东山湾独特的窄口型半封闭河口海湾的特点,基于ECOMSED模式（2002）建立了东山湾三维水动力模型,并通过第三代海浪模式MASNUM加入了波浪对底切应力及辐射应力的影响,通过ECOMSED中的底边界层模型考虑了波浪增强底摩擦的作用,综合分析了东山湾的水动力及泥沙状况.结果表明在东山湾数值模拟中,该模式能较好地模拟这类海域的水动力及泥沙输运状况.在东山湾模拟计算中,潮流的作用强于波浪的作用,但考虑波浪因素后,泥沙模拟结果更好.在波浪的作用中,底切应力相比于二维辐射应力占有绝对的优势,两者相差2个量级以上,因此可以不考虑二维辐射应力的影响.     

The concentration and distribution of dimethyl sulfide in the marine atmospheric boundary layer near the equator

Ｔｈｅｃｏｎｃｅｎｔｒａｔｉｏｎａｎｄｄｉｓｔｒｉｂｕｔｉｏｎｏｆｄｉｍｅｔｈｙｌｓｕｌｆｉｄｅｉｎｔｈｅｍａｒｉｎｅａｔｍｏｓｐｈｅｒｉｃｂｏｕｎｄａｒｙｌａｙｅｒｎｅａｒｔｈｅｅｑｕａｔｏｒ￥ＬｉＸｉｎｇｓｈｅｎｇ；ＬｉＺｈｅ；Ｆ．Ｐａｒｕｎｇｏ...     

三维二阶水波绕射问题的有限元时域计算

本文用有限元法配合时步处理来求解三维非线性水波的绕射问题,自由表面条件和物面条件都满足到二阶,采用人工阻尼区来吸收反射波,流场内的速度势通过求解有限元方程得到。对垂直圆柱体的绕射问题进行了计算,得到了自由表面波高时间历程和圆柱所受到的波浪力,计算结果和有关文献的理论计算结果进行了比较。     

Radiative and Thermal Impacts of Smoke Aerosol Longwave Absorption during Fires in the Moscow Region in Summer 2010

The aerosol longwave radiative forcing of the atmosphere and heating rate of the near-surface aerosol layer are estimated for the extreme smoke conditions in the Moscow region in summer 2010. Thermal radiation fluxes in the atmosphere are determined using the integral transmission function and semiempirical aerosol model developed on the basis of standard aerosol models and measurements at the Zvenigorod Scientific Station, Obukhov Institute of Atmospheric Physics, Russian Academy of Sciences. The aerosol radiative forcing reached 33 W/m² at the lower atmospheric boundary and ranged between–1.0 and 1.0 W/m² at the upper atmospheric boundary. The heating rate of the 10-m atmospheric layer near surface was up to 0.2 K/h during the maximum smoke conditions on August 7–9. The sensitivity of the aerosol longwave radiative forcing to the changes in the aerosol absorption coefficient and aerosol optical thickness are estimated.     

Recent improvements to the physical model of the Bohai Sea,the Yellow Sea and the East China Sea Operational Oceano-graphy Forecasting System

In order to satisfy the increasing demand for the marine forecasting capacity, the Bohai Sea, the Yellow Sea and the East China Sea Operational Oceanography Forecasting System (BYEOFS) has been upgraded and improved to Version 2.0. Based on the Regional Ocean Modeling System (ROMS), a series of comparative experiments were conducted during the improvement process, including correcting topography, changing sea surface atmospheric forcing mode, adjusting open boundary conditions, and considering atmospheric pressure correction. (1) After the topography correction, the volume transport and meridional velocity maximum of Yellow Sea Warm Current increase obviously and the unreasonable bending of its axis around 36.1°N, 123.5°E disappears. (2) After the change of sea surface forcing mode, an effective negative feedback mechanism is formed between predicted sea surface temperature (SST) by the ocean model and sea surface radiation fluxes fields. The simulation errors of SST decreased significantly, and the annual average of root-mean-square error (RMSE) decreased by about 18％. (3) The change of the eastern lateral boundary condition of baroclinic velocity from mixed Radiation-Nudging to Clamped makes the unreasonable westward current in Tsushima Strait disappear. (4) The adding of mean sea level pressure correction option which forms the mean sea level gradient from the Bohai Sea and the Yellow Sea to the western Pacific in winter and autumn is helpful to increasing the fluctuation of SLA and outflow of the Yellow Sea when the cold high air pressure system controls the Yellow Sea area.     

On the calculation of two-dimensional added mass and damping coefficients by simple Green's function technique

This paper describes a numerical method for calculating the two-dimensional hydrodynamic coefficients of one or two infinitely long, arbitrary cylinders forced to oscillate in or below the free water surface. The oscillation modes, amplitudes and phases of the cylinders may be different from one another. Finite water depth and a quay can be taken into account. Special consideration has been given to the radiation boundary conditions.The computer program developed has been tested in various two-dimensional boundary situations; it has produced results in good agreement with results obtained by other methods.     

Three-Dimensional Boundary Element Method Applied to Nonlinear Wave Transformation

For higher accuracy in simulating the transformation of three dimensional waves,in consid-eration of the advantages of constant panels and linear elements,a combined boundary elements is appliedin this research.The method can be used to remove the transverse vibration due to the accumulation ofcomputational errors.A combined boundary condition of sponge layer and Sommerfeld radiation condi-tion is used to remove the reflected waves from the computing domain.By following the water particle onthe water surface,the third order Stokes wave transform is simulated by the numerical wave flume tech-nique.The computed results are in good agreement with theoretical ones.     

海洋表层温度对台风"蔷薇"路径和强度预测精度的影响

基于中尺度大气模式WRF(Weather Research and Forecasting Model),首先对2007年3次船舶辐射通量观测进行模拟,以检验WRF对长波和短波辐射通量的模拟能力,结果表明使用中国近海海洋环境数值预报系统环流模式POM(Princeton Ocean Model)模拟的高时空分辨率的海洋表层温度能够显著改进短波辐射通量的模拟,而对长波辐射通量模拟的改进不明显。然后,将业务化运行的中国近海海洋环境数值预报系统后报的逐时海洋表面温度(SST)作为WRF底边界条件,对2008年15号强台风"蔷薇"(Jangmi)过程进行了数值后报试验。结果表明,与使用NCEP/NCAR的SST试验后报的台风中心位置偏差相比,使用高时空分辨率的SST能够较为显著地改善"蔷薇"的路径模拟,台风中心位置模拟偏差减少11%,尤其在台风减弱阶段,台风中心位置模拟偏差减少37%。台风强度在台风发展的不同阶段对下垫面SST的变化敏感性不同。台风路径附近的海表面温度下降会导致海洋向大气输送的热量减少从而减弱台风强度。