Sound speed profiles in high spatiotemporal resolution using multigrid three-dimensional variational method:a coastal experiment off northern Shandong Peninsula

It is essential to ac quire sound speed profiles(SSPs)in high-precision spatiotemporal resolution for undersea acoustic activities.However,conventional observation methods cannot obtain high-resolution SSPs.Besides,S SPs are complex and changeable in time and space,especially in coastal areas.We proposed a new space-time multigrid three-dimensional variational method with weak constraint term(referred to as STC-MG3DVar)to construct high-precision spatiotemporal resolution SSPs in coastal areas,in which sound velocity is defined as the analytical variable,and the Chen-Millero sound velocity empirical formula is introduced as a weak constraint term into the cost function of the STC-MG3DVar.The spatiotemporal correlation of sound velocity observations is taken into account in the STC-MG3DVar method,and the multi-scale information of sound velocity observations from long waves to short waves can be successively extracted.The weak constraint term can optimize sound velocity by the physical relationship between sound velocity and temperature-salinity to obtain more reasonable and accurate SSPs.To verify the accuracy of the STC-MG3DVar,SSPs observations and CTD observations(temperature observations,salinity observations)are obtained from field experiments in the northern coastal area of the Shandong Peninsula.The average root mean square error(RMSE)of the STC-MG3DVar-constructed SSPs is 0.132 m/s,and the STC-MG3DVar method can improve the SSPs construction accuracy over the space-time multigrid 3DVar without weak constraint term(ST-MG3DVar)by 10.14%and over the spatial multigrid 3DVar with weak constraint term(SC-MG3DVar)by 44.19%.With the advantage of the constraint term and the spatiotemporal correlation information,the proposed STC-MG3DVar method works better than the ST-MG3DVar and the SCMG3DVar in constructing high-precision spatiotemporal re solution SSPs.

A New Extrapolation Economy Cascadic Multigrid Method for Image Restoration Problems

In this paper, a new extrapolation economy cascadic multigrid method is proposed to solve the image restoration model. The new method combines the new extrapolation formula and quadratic interpolation to design a nonlinear prolongation operator, which provides more accurate initial values for the fine grid level. An edge preserving denoising operator is constructed to remove noise and preserve image edges. The local smoothing operator reduces the influence of staircase effect. The experiment results show that the new method not only improves the computational efficiency but also ensures good recovery quality.

水下多格弹性角反射体声反射特性及强度增强方法

单格角反射体只能反射一定方向的入射波,为了反射不同方向入射声波,需要构造多格角反射体并研究其声反射特性。以不同类型八格金属板角反射体为对象,基于结构有限元耦合间接边界元法,对水下多格弹性角反射体的散射声场进行数值仿真,分析不同频率、材料、厚度、入射角度等因素影响下不同类型多格角反射体声反射性能,从便于工程应用角度提出增强其声反射能力的方法。通过典型情况下的水池实验,对仿真计算与实验数据的一致性进行验证。实验结果表明:水下多格弹性角反射体目标强度随入射频率增大而增大,并趋向稳定;目标强度随入射角度改变波动幅度较大,不同类型角反射体目标强度值差异较大;目标强度随平板厚度的增加而增大,但不显著;目标强度随着边长的增加而明显增大;软介质材料构成的多格角反射体,水下声反射性能明显优于金属角反射体。因此,为增强水下多格角反射体声反射能力,增大反射体边长及选用低阻抗材料是更有效的方法。

基于多重网格算法的大地电磁三维反演研究

为了推动大地电磁三维反演实用化的进程,本文融合了多重网格法(MG)与有限内存BFGS(L-BFGS)拟牛顿法,提出了一种高效的多尺度三维反演方法.该方法首先从细网格开始,逐步粗化生成最粗网格,然后从最粗网格开始反演,并逐级将粗网格反演结果为上一级细网格提供一个相对合理的初始模型,以加速细网格的反演迭代过程,并同时得到由粗到细的多尺度反演结果.本文首先测试了棋盘模型、起伏地层和侵入带模型,从理论上验证了该算法的高效性.结果表明,本文算法计算效率是传统反演的2~3倍.然后,利用该算法对美国Cascadia地区实测数据进行反演,进一步验证了该算法的实用性,结果显示该算法计算效率是传统反演的近2倍.期待本研究可为大地电磁三维反演问题提供一种高效、精确的解决方案.

Steklov-Lame特征值问题自适应多网格方法的后验误差估计

建立Steklov-Lame特征值问题的一种基于移位反迭代的有限元多网格离散方案,并研究该方案基于残差型的后验误差估计.首先给出近似特征函数在L^(2)(∂Ω)范数意义下的误差估计,其次给出多网格方案近似解的后验误差指示子,并证明后验误差指示子的可靠性和有效性.最后利用后验误差指示子设计自适应多网格算法并用于求解Steklov-Lame特征值问题.

基于弹流润滑的滚子轴承参数研究

基于弹流润滑理论研究滚子轴承结构参数对其润滑特性的影响。结合弹性接触变形方程以及流体动力学润滑方程,建立适合于滚子轴承的弹流润滑模型,研究在不同椭圆率、载荷、卷吸速度以及黏度等因素作用下滚子轴承的摩擦学性能变化规律。结果表明:椭圆率、载荷、卷吸速度以及黏度会不同程度地影响压力峰值及二次压力峰等参数;随着椭圆率的增大,油膜厚度以及压力显著增大;随着载荷的增大,总体压力、压力峰值、二次压力峰及其尖锐度明显增大,但最小油膜厚度略有下降;黏度与最小油膜厚度以及压力存在着明显正相关关系;卷吸速度与油膜厚度存在着微弱正相关关系,与油膜压力存在着微弱负相关关系。因此,一定程度上增大椭圆率并减小载荷,有利于提高润滑性能。

基于余维裂缝多孔介质渗流的差分多重网格法

研究了余维(又称为降维)裂缝多孔介质中不可压单相流体渗流问题的一种有效数值求解方法。结合多重网格法和块中心有限差分法的优点,给出了一种块中心有限差分多重网格法求解裂缝多孔介质渗流问题,通过该方法提高了计算效率。经过数值实验,验证了该方法的有效性。

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点。本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、Pre_V和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性。

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

To develop an efficient and robust aerodynamic analysis method for numerical optimization designs of wing and complex configuration, a combination of matrix preconditioning and multigrid method is presented and investigated. The time derivatives of three-dimensional Navier-Stokes equations are preconditioned by Choi-Merkle preconditioning matrix that is originally designed for two-dimensional low Mach number viscous flows. An extension to three-dimensional viscous flow is implemented, and a method improving the convergence for transonic flow is proposed. The space discretizaition is performed by employing a finite-volume cell-centered scheme and using a central difference. The time marching is based on an explicit Rtmge-Kutta scheme proposed by Jameson. An efficient FAS multigrid method is used to accelerate the convergence to steady-state solutions. Viscous flows over ONERA M6 wing and M100 wing are numerically simulated with Mach numbers ranging from 0.010 to 0.839. The inviscid flow over the DLR-F4 wing-body configuration is also calculated to preliminarily examine the performance of the presented method for complex configuration. The computed results are compared with the experimental data and good agreement is achieved. It is shown that the presented method is efficient and robust for both compressible and incompressible flows and is very attractive for aerodynamic optimization designs of wing and complex configuration.

A theoretical study of the multigrid three-dimensional variational data assimilation scheme using a simple bilinear interpolation algorithm

In order to solve the so-called ＂bull-eye＂ problem caused by using a simple bilinear interpolation as an observational mapping operator in the cost function in the multigrid three-dimensional variational （3DVAR） data assimilation scheme, a smoothing term, equivalent to a penalty term, is introduced into the cost function to serve as a means of troubleshooting. A theoretical analysis is first performed to figure out what on earth results in the issue of ＂bull-eye＂, and then the meaning of such smoothing term is elucidated and the uniqueness of solution of the multigrid 3DVAR with the smoothing term added is discussed through the theoretical deduction for one-dimensional （1D） case, and two idealized data assimilation experiments （one- and two-dimensional （2D） cases）. By exploring the relationship between the smoothing term and the recursive filter theoretically and practically, it is revealed why satisfied analysis results can be achieved by using such proposed solution for the issue of the multigrid 3DVAR.

MULTIGRID TECHNIQUE APPLIED TO LINEAR TIME-DEPENDENT HYPERBOLIC SYSTEMS

A system of linear time-dependent hyperbolic partial differential equations in the form of the time-domain Maxwell′s equations is numerically solved using ageometric multigrid method.The multilevel method is an adaptation of Ni′s cell-vertex based multigrid technique,originally proposed for accelerating steady state convergence of nonlinear time-dependent Euler equations of gas dynamics.We discuss issues pertaining to the application of the geometric multigrid method to a system of equations where the major issue is of accurately propagating linear waves over large distances leading to major constraints on the required grid resolution in terms of points-perwavelength.

L^2-ERROR OF EXTRAPOLATION CASCADIC MULTIGRID (EXCMG)

Based on an asymptotic expansion of finite element, an extrapolation cascadic multigrid method （EXCMG） is proposed, in which the new extrapolation and quadratic interpolation are used to provide a better initial value on refined grid. In the case of multiple grids, both superconvergence error in H^1-norm and the optimal error in l2-norm are analyzed. The numerical experiment shows the advantage of EXCMG in comparison with CMG.

Biorthogonal Wavelet Based Algebraic Multigrid Preconditioners for Large Sparse Linear Systems

In this article algebraic multigrid as preconditioners are designed, with biorthogonal wavelets, as intergrid operators for the Krylov subspace iterative methods. Construction of hierarchy of matrices in algebraic multigrid context is based on lowpass filter version of Wavelet Transform. The robustness and efficiency of this new approach is tested by applying it to large sparse, unsymmetric and ill-conditioned matrices from Tim Davis collection of sparse matrices. Proposed preconditioners have potential in reducing cputime, operator complexity and storage space of algebraic multigrid V-cycle and meet the desired accuracy of solution compared with that of orthogonal wavelets.

Standard and Economical Cascadic Multigrid Methods for the Mortar Finite Element Methods

In this paper,standard and economical cascadic multigrid methods are considered for solving the algebraic systems resulting from the mortar finite element methods. Both cascadic multigrid methods do not need full elliptic regularity,so they can be used to tackle more general elliptic problems.Numerical experiments are reported to support our theory.

3-D direct current resistivity forward modeling by adaptive multigrid finite element method

Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.

System of Multigrid Nonlinear Least-squares Four-dimensional Variational Data Assimilation for Numerical Weather Prediction(SNAP):System Formulation and Preliminary Evaluation

A new forecasting system-the System of Multigrid Nonlinear Least-squares Four-dimensional Variational(NLS-4DVar)Data Assimilation for Numerical Weather Prediction(SNAP)-was established by building upon the multigrid NLS-4DVar data assimilation scheme,the operational Gridpoint Statistical Interpolation(GSI)−based data-processing and observation operators,and the widely used Weather Research and Forecasting numerical model.Drawing upon lessons learned from the superiority of the operational GSI analysis system,for its various observation operators and the ability to assimilate multiple-source observations,SNAP adopts GSI-based data-processing and observation operator modules to compute the observation innovations.The multigrid NLS-4DVar assimilation framework is used for the analysis,which can adequately correct errors from large to small scales and accelerate iteration solutions.The analysis variables are model state variables,rather than the control variables adopted in the conventional 4DVar system.Currently,we have achieved the assimilation of conventional observations,and we will continue to improve the assimilation of radar and satellite observations in the future.SNAP was evaluated by case evaluation experiments and one-week cycling assimilation experiments.In the case evaluation experiments,two six-hour time windows were established for assimilation experiments and precipitation forecasts were verified against hourly precipitation observations from more than 2400 national observation sites.This showed that SNAP can absorb observations and improve the initial field,thereby improving the precipitation forecast.In the one-week cycling assimilation experiments,six-hourly assimilation cycles were run in one week.SNAP produced slightly lower forecast RMSEs than the GSI 4DEnVar(Four-dimensional Ensemble Variational)as a whole and the threat scores of precipitation forecasts initialized from the analysis of SNAP were higher than those obtained from the analysis of GSI 4DEnVar.

A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids

We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.

An efficient multigrid-DEIM semi-reduced-order model for simulation of single-phase compressible flow in porous media

In this paper,an efcient multigrid-DEIM semi-reduced-order model is developed to accelerate the simulation of unsteady single-phase compressible fow in porous media.The cornerstone of the proposed model is that the full approximate storage multigrid method is used to accelerate the solution of fow equation in original full-order space,and the discrete empirical interpolation method(DEIM)is applied to speed up the solution of Peng-Robinson equation of state in reduced-order subspace.The multigrid-DEIM semi-reduced-order model combines the computation both in full-order space and in reducedorder subspace,which not only preserves good prediction accuracy of full-order model,but also gains dramatic computational acceleration by multigrid and DEIM.Numerical performances including accuracy and acceleration of the proposed model are carefully evaluated by comparing with that of the standard semi-implicit method.In addition,the selection of interpolation points for constructing the low-dimensional subspace for solving the Peng-Robinson equation of state is demonstrated and carried out in detail.Comparison results indicate that the multigrid-DEIM semi-reduced-order model can speed up the simulation substantially at the same time preserve good computational accuracy with negligible errors.The general acceleration is up to 50-60 times faster than that of standard semi-implicit method in two-dimensional simulations,but the average relative errors of numerical results between these two methods only have the order of magnitude 10^(−4)-10^(−6)%.

MULTIGRID FOR THE MORTAR ELEMENT METHOD WITH LOCALLY P_1 NONCONFORMING ELEMENTS

In this paper we study the theoretical properties of multigrid algorithm for discretization of the Poisson equation in 2D using a mortar element method under the assumption that the triangulations on every subdomain are uniform.We prove the convergence of the W-cycle with a sufficiently large number of smoothing steps.The variable V-cycle multigrid preconditioner are also available.

A New Multigrid 3D-VAR Optimization Method for Bottom Friction Using HF Radar Current Observation

This paper proposes a new method for data assimilation of the surface radial current observed by High Frequency ground wave radar and optimization of the bottom friction coefficient.In this method,the shallow water wave equation is introduced into the cost function of the multigrid three-dimensional variation data assimilation method as the weak constraint term,the surface current and the bottom friction coefficient are defined as the analytical variables,and the high spatiotemporal resolution surface radial flow observed by the high-frequency ground wave radar is used to optimize the surface current and bottom friction coefficient.This method can effectively consider the spatiotemporal correlation of radar data and extract multiscale information from surface radial flow data from long waves to short waves.Introducing the shallow water wave equation into the cost function as a weak constraint condition can adjust both the momentum and mass fields simultaneously to obtain more reasonable analysis information.The optimized bottom friction coefficient is introduced into the regional ocean numerical model to carry out numerical experiments.The test results show that the bottom friction coefficient obtained by this method can effectively improve the accuracy of the numerical simulation of sea surface height in the offshore area and reduce the simulation error.