A primitive-variable Riemann method for solution of the shallow water equations with wetting and drying

A Riemann flux that uses primitive variables rather than conserved variables is developed for the shallow water equations with nonuniform bathymetry. This primitive-variable flux is both conservative and well behaved at zero depth. The unstructured finite-volume discretization used is suitable for highly nonuniform grids that provide resolution of complex geometries and localized flow structures. A source-term discretization is derived for nonuniform bottom that balances the discrete flux integral both for still water and in dry regions. This primitive-variable formulation is uniformly valid in wet and dry regions with embedded wetting and drying fronts. A fully nonlinear implicit scheme and both nonlinear and time-linearized explicit schemes are developed for the time integration. The implicit scheme is solved by a parallel Newton-iterative algorithm with numerically computed flux Jacobians. A concise treatment of characteristic-variable boundary conditions with source terms is also given. Computed results obtained for the one-dimensional dam break on wet and dry beds and for normal-mode oscillations in a circular parabolic basin are in very close agreement with the analytical solutions. Other results for a forced breaking wave with friction interacting with a sloped bottom demonstrate a complex wave motion with wetting, drying and multiple interacting wave fronts. Finally, a highly nonuniform, coastline-conforming unstructured grid is used to demonstrate an unsteady simulation that models an artificial coastal flooding due to a forced wave entering the Gulf of Mexico.     

基于小波变换和数学形态学的运动物体检测

提出了一种基于帧差法和小波变换相结合的运动目标检测方法,充分利用帧差法计算简单和小波变换的多尺度特性。实验表明,这种方法可以有效地从复杂自然场景的图像序列中检测出完整的运动目标。而且能够有效的抑制噪声。同时减少计算时间,满足检测的实时性要求。     

海态雷达的一种矢量海流合成方法及其计算机模拟

根据海洋表面流在相邻海域连续的特性和矢量合成原理,提出了一种生成海洋表面矢量流图的算法。该算法只在合成有效区域内才定义网格。这种算法减少了运算量,节省了存储空间,加快了运算速度。并通过计算机模拟显示,其结果能正确反映海流的实时状态。     

一个带不确定权的积分方程组解的对称性

结合积分形式移动平面法的思想,讨论Rn上积分方程组u(x)=∫Rn|x-y|α-na(y)v(y)qdy,v(x)=∫Rn|x-y|α-nb(y)u(y)pdy的正解关于某一点的对称性和单调性,其中0αn,p,q1,p+11+q+11=n n-α,a(x)和b(x)满足一些对称性、单调性.     

A new method of moving asymptotes for large-scale linearly equality-constrained minimization

A new method of moving asymptotes for large-scale minimization subject to linear equality constraints is discussed.In this method,linear equality constraints are deleted with null space technique and the descending direction is obtained by solving a convex separable subproblem of moving asymptotes in each iteration. New rules for controlling the asymptotes parameters are designed and the global convergence of the method under some reasonable conditions is established and proved.The numerical results show that the new method may be capable of processing some large scale problems.     

Unstructured Adaptive Grid Flow Simulations of Inert and Reactive Gas Mixtures

Unstructured adaptive grid flow simulation is applied to the calculation of high-speed compressible flows of inert and reactive gas mixtures. In the present case, the flowfield is simulated using the 2-D Euler equations, which are discretized in a cell-centered finite volume procedure on unstructured triangular meshes. Interface fluxes are calculated by a Liou flux vector splitting scheme which has been adapted to an unstructured grid context by the authors. Physicochemical properties are functions of the local mixture composition, temperature, and pressure, which are computed using the CHEMKIN-II subroutines. Computational results are presented for the case of premixed hydrogen–air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave fronts. A parametric analysis of the influence of the adaptation parameters on the computed solution is performed.     

Radial symmetry and regularity of solutions for poly-harmonic Dirichlet problems

Let B=B₁(0) be the unit ball in Rn and r=|x|. We study the poly-harmonic Dirichlet problem

     

基于ARIMA与神经网络集成的GDP时间序列预测研究

本文深入分析了单整自回归移动平均(ARIMA)模型与神经网络(NN)模型的预测特性和优劣,并在此基础上建立了由ARIMA模型和NN模型集成的GDP时间序列预测模型与算法。其基本思想是充分发挥两种模型在线性空间和非线性空间的预测优势,据此将GDP时间序列的数据结构分解为线性自相关主体和非线性残差两部分,首先用ARIMA模型预测序列的线性主体,然后用NN模型对其非线性残差进行估计,最终集成为整个序列的预测结果。仿真实验表明:集成模型的预测准确率显著高于单一模型的预测准确率,从而证实了集成模型用于GDP预测的有效性。     

带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

Local and parallel finite element algorithm for stationary incompressible magnetohydrodynamics

This article presents a local and parallel finite element method for the stationary incompressible magnetohydrodynamics problem. The key idea of this algorithm comes from the two‐grid discretization technique. Specifically, we solve the nonlinear system on a global coarse mesh, and then solve a series of linear problems on several subdomains in parallel. Furthermore, local a priori estimates are obtained on a general shape regular grid. The efficiency of the algorithm is also illustrated by some numerical experiments.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1513–1539, 2017