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 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.     

无结构网格上二维浅水流动的数值模拟

提出在无结构网格上建立有限体积高性能格式族的统一框架,通过引入跨单元界面法向数值通量的逆风分解,将一维Osher、TVD两种通量分裂格式自然地推广至二维浅水方程组。给出了各种情况下浅水方程组有限体积法边界处理的计算公式。最后,利用该格式对穹包溢流、陆地动边界和河口潮流三个问题进行计算。结果表明,这类守恒格式具有高精度、无振荡性以及处理复杂流态过渡、自动捕俘间断和模拟陆地动边界的功能。     

A mixed finite element method for the solution of the magnetostatic problem with highly discontinuous coefficients in 3D

This paper present a new a mixed finite element method for the simulation of magnetostatic problems with highly discontinuous permeability. The method is derived from the well studied mixed formulation for the div-grad system that is known to be accurate for very large discontinuities. The method robustness is demonstrated on a test model problem.     

Bragg resonant reflection of carrier waves composing wave groups

Numerical analyses for the Bragg resonant reflection of carrier waves associated long waves due to sinusoidally varying seabeds are performed by using a set of coupled ordinary differential equations derived from the Boussinesq equations. The Boussinesq equations are firstly approximated with the Fourier decomposition. The coupled governing equations are then derived and used to simulate evolution of both short and long wave components. It is also found that wave groups are generated by two carrier waves with slightly different frequencies. The wave energy of the initial wave components is transferred to other harmonic components during propagation over a long distance. Evolution and reflection of both short and long waves were largely affected by nonlinearity.     

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

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

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.     

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?????????     

关于二阶伴随模型的理论研究

Hesse矩阵-目标函数关于控制变量的二阶偏导数形成的矩阵,在变分数据同化过程中以及敏感性分析等方面起着重要的作用;它可以通过建立数学模型的一阶和二阶伴随模型求得.以浅水方程模型为例,利用泛函的Gâteaux微分和Hilbert空间上伴随算子的概念,导出了一阶和二阶伴随模型并由此得到Hesse矩阵.改进了Zhi Wang等(1992)建立的二阶伴随模型理论.