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.

Three-Dimensional Space Interpolation of Grey / Depth Image Sequence-A New Technique of Computer Graphics Synthesis

This paper advances a three-dimensional space interpolation method of grey / depth image sequence, which breaks free from the limit of original practical photographing route. Pictures can cruise at will in space. By using space sparse sampling, great memorial capacity can be saved and reproduced scenes can be controlled. To solve time consuming and complex computations in three-dimensional interpolation algorithm, we have studied a fast and practical algorithm of scattered space lattice and that of 'Warp' algorithm with proper depth. By several simple aspects of three dimensional space interpolation, we succeed in developing some simple and practical algorithms. Some results of simulated experiments with computers have shown that the new method is absolutely feasible.

Accuracy and Computational Cost of Interpolation Schemes While Performing <i>N</i>-Body Simulations

The continuous approximations play a vital role in N-body simulations. We constructed three different types, namely, one-step (cubic and quintic Hermite), two-step, and three-step Hermite interpolation schemes. The continuous approximations obtained by Hermite interpolation schemes and interpolants for ODEX2 and ERKN integrators are discussed in this paper. The primary focus of this paper is to measure the accuracy and computational cost of different types of interpolation schemes for a variety of gravitational problems. The gravitational problems consist of Kepler’s two-body problem and the more realistic problem involving the Sun and four gas-giants—Jupiter, Saturn, Uranus, and Neptune. The numerical experiments are performed for the different integrators together with one-step, two-step, and three-step Hermite interpolation schemes, as well as the interpolants.

Efficient Difference Schemes for the Caputo-Tempered Fractional Diffusion Equations Based on Polynomial Interpolation

The tempered fractional calculus has been successfully applied for depicting the time evolution of a system describing non-Markovian diffusion particles.The related governing equations are a series of partial differential equations with tempered fractional derivatives.Using the polynomial interpolation technique,in this paper,we present three efficient numerical formulas,namely the tempered L1 formula,the tempered L1-2 formula,and the tempered L2-1_(σ)formula,to approximate the Caputo-tempered fractional derivative of orderα∈(0,1).The truncation error of the tempered L1 formula is of order 2-α,and the tempered L1-2 formula and L2-1_(σ)formula are of order 3-α.As an application,we construct implicit schemes and implicit ADI schemes for one-dimensional and two-dimensional time-tempered fractional diffusion equations,respectively.Furthermore,the unconditional stability and convergence of two developed difference schemes with tempered L1 and L2-1_(σ)formulas are proved by the Fourier analysis method.Finally,we provide several numerical examples to demonstrate the correctness and effectiveness of the theoretical analysis.

Splitting Based Scheme for Three-dimensional MHD with Dual Time Stepping

A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.

Construction and Application of Subdivision Surface Scheme Using Lagrange Interpolation Polynomial

This paper offers a general formula for surface subdivision rules for quad meshes by using 2-D Lagrange interpolating polynomial [1]. We also see that the result obtained is equivalent to the tensor product of (2N + 4)-point n-ary interpolating curve scheme for N ≥ 0 and n ≥ 2. The simple interpolatory subdivision scheme for quadrilateral nets with arbitrary topology is presented by L. Kobbelt [2], which can be directly calculated from the proposed formula. Furthermore, some characteristics and applications of the proposed work are also discussed.

Cures for expansion shock and shock instability of Roe scheme based on momentum interpolation mechanism

The common defects of the Roe scheme are the non-physical expansion shock and shock instability. By removing the momentum interpolation mechanism(MIM), an improved method with several advantages has been presented to suppress the shock instability. However, it cannot prevent the expansion shock and is incompatible with the traditional curing method for expansion shock. To solve the problem, the traditional curing mechanism is analyzed. Effectiveness of the traditional curing method is discussed,and several defects are identified, one of which leads to incompatibility between curing shock instability and expansion shock. Consequently, an improved Roe scheme is proposed, which is with low computational costs, concise, easy to implement, and robust.More importantly, the proposed scheme can simultaneously solve the problem of shock instability and expansion shock without additional costs.

A Multi-Secret Sharing Scheme with Many Keys Based on Hermite Interpolation

A secret sharing scheme is one of cryptographies. A threshold scheme, which is introduced by Shamir in 1979, is very famous as a secret sharing scheme. We can consider that this scheme is based on Lagrange's interpolation formula. A secret sharing scheme has one key. On the other hand, a multi-secret sharing scheme has more than one keys;that is, a multi-secret sharing scheme has p (≥2) keys. Dealers distribute shares of keys among n participants. Gathering t (≤n) participants, keys can be reconstructed. In this paper, we give a scheme of a (t,n) multi-secret sharing based on Hermite interpolation, in the case of p≤t.

Using Three-Dimensional Lorenz Scatter Plots to Detect Patients with Atrioventricular Node Double Path Caused by Interpolated Ventricular Premature Systoles: A Case Study

A series of related electrophysiology phenomena can be caused by the occurrence of interpolated ventricular premature contraction.In our recent three-dimensional Lorenz R-R scatter plot research,we found that atrioventricular node double path caused by interpolated ventricular premature contraction imprints a specifi c pattern on three-dimensional Lorenz plots generated from 24-hour Holter recordings.We found two independent subclusters separated from the interpolated premature beat precluster,the interpolated premature beat cluster,and the interpolated premature beat postcluster,respectively.Combined with use of the trajectory tracking function and the leap phenomenon,our results reveal the presence of the atrioventricular node double conduction path.

Some Experiments with Multivariate Objective Analysis Scheme of Heights and Winds Using Optimum Interpolation

A two-dimensional, multitvariate objective analysis scheme for simultaneous analysis of geopotential height and wind fields has been developed over Indian and adjoining region for use in numerical weather prediction. The height-height correlations calculated using daily data of four July months (1976-1979), are used to derive the other autocorrelations and cross-correlations assuming geostropic relationship. A Gaussian function is used to model the autocorrelation function. Since the scheme is multivariate the regression coefficients (weights) are matrix.Near the equator, the geostrophic approximation relating mass and wind is decoupled in a way similar to Bergman (1979). The objective analyses were made over Indian and adjoining region for 850, 700, 500, 300 and 200 hPa levels for the period from 4 July to 8 July 1979, 12 GMT. The analyses obtained using multivariate optimum interpolation scheme depict the synoptic situations satisfactorily. The analyses were also compared with the FGGE analyses (from ECMWF) and also with the station observations by computing the root mean square (RMS) errors and the RMS errors are comparable with those obtained in other similar studies.

A Verifiable Multi-Secret Sharing Scheme Based on Hermite Interpolation

A threshold scheme, which is introduced by Shamir in 1979, is very famous as a secret sharing scheme. We can consider that this scheme is based on Lagrange＇s interpolation formula. A secret sharing scheme has one key. On the other hand, a multi-secret sharing scheme has more than one key, that is, a multi-secret sharing scheme has p （〉_ 2） keys. Dealer distribute shares of keys among n participants. Gathering t （〈 n） participants, keys can be reconstructed. Yang et al. （2004） gave a scheme of a （t, n） multi-secret sharing based on Lagrange＇s interpolation. Zhao et al. （2007） gave a scheme of a （t, n） verifiable multi-secret sharing based on Lagrange＇s interpolation. Recently, Adachi and Okazaki give a scheme of a （t, n） multi-secret sharing based on Hermite interpolation, in the case ofp 〈 t. In this paper, we give a scheme ofa （t, n） verifiable multi-secret sharing based on Hermite interpolation.

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).

Three-dimensional hydrodynamic model of Xiamen waters

A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.

A CLASS OF TWO-LEVEL EXPLICIT DIFFERENCE SCHEMES FOR SOLVING THREE DIMENSIONAL HEAT CONDUCTION EQUATION

A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.

A Unified Interpolating Subdivision Scheme for Curves/Surfaces by Using Newton Interpolating Polynomial

This paper presents a general formula for (2m + 2)-point n-ary interpolating subdivision scheme for curves for any?integer m ≥ 0 and n ≥ 2 by using Newton interpolating polynomial. As a consequence, the proposed work is extended for surface case, which is equivalent to the tensor product of above proposed curve case. These formulas merge several notorious curve/surface schemes. Furthermore, visual performance of the subdivision schemes is also presented.

Proficiency of Second Derivative Schemes for the Numerical Solution of Stiff Systems

This paper presents a study on the development and implementation of a second derivative method for the solution of stiff first order initial value problems of ordinary differential equations using method of interpolation and collocation of polynomial approximate solution. The results of this paper bring some useful information. The constructed methods are A-stable up to order 8. As it is shown in the numerical examples, the new methods are superior for stiff systems.

Studies on Non-interpolating Semi-Lagrangian Scheme and Numerical Solution to KdV Equation

A new non-interpolating semi-Lagrangian scheme has been proposed, which can eliminate any interpolation,and consequently numerical smoothing of forecast fields. Here the new scheme is applied to KdV equation and its performance is assessed by comparing the numerical results with those produced by Ritchie＇s scheme (1986).The comparison shows that the non-interpolating semi-Lagrangian scheme appears to have efficiency advantages.

Improvement of the Semi-Lagrangian Advection Scheme in the GRAPES Model:Theoretical Analysis and Idealized Tests

ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.

Optimised Spatial Sampling Scheme for Soil Electriclal Conductivity Based on Variance Quad-Tree (VQT) Method

The acquisition of precise soil data representative of the entire survey area, is a critical issue for many treatments such as irrigation or fertilization in precision agriculture. The aim of this study was to investigate the spatial variability of soil bulk electrical conductivity （ECb） in a coastal saline field and design an optimized spatial sampling scheme of ECb based on a sampling design algorithm, the variance quad-tree （VQT） method. Soil ECb data were collected from the field at 20 m interval in a regular grid scheme. The smooth contour map of the whole field was obtained by ordinary kriging interpolation, VQT algorithm was then used to split the smooth contour map into strata of different number desired, the sampling locations can be selected within each stratum in subsequent sampling. The result indicated that the probability of choosing representative sampling sites was increased significantly by using VQT method with the sampling number being greatly reduced compared to grid sampling design while retaining the same prediction accuracy. The advantage of the VQT method is that this scheme samples sparsely in fields where the spatial variability is relatively uniform and more intensive where the variability is large. Thus the sampling efficiency can be improved, hence facilitate an assessment methodology that can be applied in a rapid, practical and cost-effective manner.

Optimal Evacuation Scheme Based on Dam-Break Flood Numerical Simulation

The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area I is flooded at 2 500 s, and after 5 500 s, most of industrial area II is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area I and II. Based on the above analysis, the optimal evacuation scheme is determined.