Numerical simulation of wave-induced currents combined with parabolic mild-slope equation in curvilinear coordinates

Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

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.     

A Time-Domain Coupled Model for Nonlinear Wave Forces on A Fixed Box-Shaped Ship in A Harbor

A 2-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear wave forces on a fixed box-shaped ship in a harbor.The domain is divided into an inner domain and an outer domain.The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations.The other area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions.Along the interface boundaries between the inner domain and the outer domain,the volume flux is assumed to be continuous and the wave pressures are equal.Relevant physical experiment is conducted to validate the present model.It is shown that the numerical results agree with the experimental data.Compared with the coupled model with the flow in the inner domain governed by the Laplace equation,the present coupled model is more efficient and its solution procedure is more simple,which is particularly useful for the study on the effect of the nonlinear wave forces on a fixed box-shaped ship in a large harbor.     

��Ȩ������С���˵ĵ����ⷨ

?????????????С???????????????????????С???????????????С???????????????????????????????????????С?????????????????????????С??????????????????й?????????????????????????y?????????     

The tripole: A new coherent vortex structure of incompressible two-dimensional flows

Abstract

Using a contour dynamical algorithm, we have found rotating tripolar V-state solutions for the inviscid Euler equations in two-dimensions. We have studied their geometry as a function of their physical parameters. Their stability was investigated with the aid of contour surgery, and most of the states were found to be stable. Under finite-amplitude perturbations, tripoles are shown to either fission into two asymmetric dipoles or to evolve into a shielded axisymmetric vortex, demonstrating the existence of two new ‘‘reversible transitions'’ between topologically distinct coherent vortex structures. These dynamical results are confirmed by pseudo-spectral simulations, with which we also show how continuous tripolar long-lived coherent vortex structures can be generated in a variety of ways.     

广义相对论中静态荷电球体的内解

本文在假设荷电球体内部物质密度为ρ_m=μγ~α,电荷密度为ρ_e=ρ_0γ~(b_e-λ/2)的情况下,严格求解Einstein-Maxwell场方程,求得了静态荷电球体的一个较为普遍的内部解。这个解可以看作是Wyman得到的内部解推广到荷电情况下的结果,并且将Wilson,李鉴增,Santos,Pant和Sah等人给出的许多结果作为自己的特例包括在内。     

双时间步法在二维浅水方程求解中的应用

针对传统算法效率低的问题,将隐式双时间步法应用于求解二维浅水方程,建立了非结构网格下高效的有限体积模型。在应用双时间步法时,虚拟时间层中的定常问题采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法进行迭代求解。通过模拟计算4个典型算例以及与传统显式算法进行比较,对模型精度、效率及处理实际问题能力进行检验,分析了时间步长、内迭代次数对模型性能的影响。结果表明,双时间步法放宽了稳定性对时间步长的限制,时间步长可取到显式格式10倍以上,计算耗时减少了50%以上,模型具有良好精度与适应性,具有较好的推广应用价值。     

辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。