On the Use of Monotonicity-Preserving Interpolatory Techniques in Multilevel Schemes for Balance Laws

Cost-effective multilevel techniques for homogeneous hyperbolic conservation laws are very successful in reducing the computational cost associated to high resolution shock capturing numerical schemes.Because they do not involve any special data structure,and do not induce savings in memory requirements,they are easily implemented on existing codes and are recommended for 1D and 2D simulations when intensive testing is required.The multilevel technique can also be applied to balance laws,but in this case,numerical errors may be induced by the technique.We present a series of numerical tests that point out that the use of monotonicity-preserving interpolatory techniques eliminates the numerical errors observed when using the usual 4-point centered Lagrange interpolation,and leads to a more robust multilevel code for balance laws,while maintaining the efficiency rates observed forhyperbolic conservation laws.

Deployment Dynamic Modeling and Driving Schemes for a Ring-Truss Deployable Antenna

Mesh reflector antennas are widely used in space tasks owing to their light weight,high surface accuracy,and large folding ratio.They are stowed during launch and then fully deployed in orbit to form a mesh reflector that transmits signals.Smooth deployment is essential for duty services;therefore,accurate and efficient dynamic modeling and analysis of the deployment process are essential.One major challenge is depicting time-varying resistance of the cable network and capturing the cable-truss coupling behavior during the deployment process.This paper proposes a general dynamic analysis methodology for cable-truss coupling.Considering the topological diversity and geometric nonlinearity,the cable network's equilibrium equation is derived,and an explicit expression of the time-varying tension of the boundary cables,which provides the main resistance in truss deployment,is obtained.The deployment dynamic model is established,which considers the coupling effect between the soft cables and deployable truss.The effects of the antenna's driving modes and parameters on the dynamic deployment performance were investigated.A scaled prototype was manufactured,and the deployment experiment was conducted to verify the accuracy of the proposed modeling method.The proposed methodology is suitable for general cable antennas with arbitrary topologies and parameters,providing theoretical guidance for the dynamic performance evaluation of antenna driving schemes.

Stability Analysis of Inverse Lax-Wendroff Procedure for a High order Compact Finite Difference Schemes

This paper considers the finite difference(FD)approximations of diffusion operators and the boundary treatments for different boundary conditions.The proposed schemes have the compact form and could achieve arbitrary even order of accuracy.The main idea is to make use of the lower order compact schemes recursively,so as to obtain the high order compact schemes formally.Moreover,the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform(FFT).With mathematical induction,the eigenvalues of the proposed differencing operators are shown to be bounded away from zero,which indicates the positive definiteness of the operators.To obtain numerical boundary conditions for the high order schemes,the simplified inverse Lax-Wendroff(SILW)procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method.Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.

High Order IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

Analysis and Comparison of Slope-cutting Widening Schemes in Highway Reconstruction and Expansion Project Based on MIDAS Software

In this paper,the geological condition of the right-side slope of the K114+694–K115+162 section of Yong-tai-wen Expressway is investigated and analyzed with the results showing that the strength of rock mass is the main contributor to the stability of the slope.Then,two widening schemes are proposed,which are the steep slope with strong support and the gentle slope with general support schemes.The static/slope module of MIDAS GTS finite element analysis software and the strength reduction method were used to compare the two schemes.The results show that the steep slope with a strong support scheme has obvious advantages in land requisition,environmental protection,and safety and is more suitable for reconstructing and expanding the highway slope.

基于半自治agent的profit-sharing增强学习方法研究

在基于半自治agent的系统中应用profit-sharing增强学习方法,并与基于动态规划的Q-learning增强学习方法进行比较,在不确定因素较多的动态环境中,当系统状态变化不是一个马尔科夫过程时profit-sharing方法具有很大优势。根据半自治agent中半自治的特性——受制性,提出了一种面向基于半自治agent的增强学习模型,以战场仿真中安全隐蔽的寻找模型为实例对基于半自治agent的profit-sharing增强学习模型进行了试验分析。

Stability of finite difference numerical simulations of acoustic logging-while-drilling with different perfectly matched layer schemes

In acoustic logging-while-drilling （ALWD） finite difference in time domain （FDTD） simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns （radius -27 mm）. Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole （the difference is 〉30 times）, the stability and efficiency of the perfectly matched layer （PML） scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain （FDTD） in the ALWD simulation, including field-splitting PML （SPML）, multiaxial PML（M- PML）, non-splitting PML （NPML）, and complex frequency-shifted PML （CFS-PML）. The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor （a0） and the scaling factor （β0） depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.

NUMERICAL SIMULATION OF SUPERSONIC AXISYMMETRIC FLOW OVER MISSILE AFTERBODY WITH JET EXHAUST USING POSITIVE SCHEMES

Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.

METHOD FOR CONSTRUCTING TAG-KEM SCHEMES WITH SHORT-MESSAGE PUBLIC-KEY ENCRYPTIONS

Tag key encapsulation mechanism （Tag-KEM）/data encapsulation mechanism （DEM） is a hybrid framework proposed in 2005. Tag-t（EM is one of its parts by using public-key encryption （PKE） technique to encapsulate a symmetric key. In hybrid encryptions, the long-raessage PKE is not desired due to its slow operation. A general method is presented for constructing Tag-KEM schemes with short-message PKEs. The chosen ciphertext security is proved in the random oracle model. In the method, the treatment of the tag part brings no additional ciphertext redundancy. Among all the methods for constructing Tag-KEM, the method is the first one without any validity checking on the tag part, thus showing that the Tag-KEM/DEM framework is superior to KEM＋DEM one.

A Family of High_order Accuracy Explicit Difference Schemes for Solving 2-D Parabolic Partial Differential Equation

A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).

Application of the two different salinity parameterization schemes in the sea ice model

In this study, we mainly introduce two salinity parameterization schemes used in Sea Ice Simulator （SIS）, that is, isosaline scheme and salinity profile scheme. Comparing the equation of isosaline scheme with that of salinity profile scheme, we found that there was one different term between the two schemes named the salinity different term. The thermodynamic effect of the salinity difference term on sea ice thickness and sea ice concentration showed that： in the freezing processes from November to next May, the sea ice temperature could rise on the influence of the salinity difference term and restrain sea ice freezing; at the first melting phase from June to August, the upper ice melting rate was faster than the lower ice melting rate. Then sea ice temperature could rise and accelerate the sea ice melting; at the second melting phase from September to October, the upper ice melting rate was slower than the lower ice melting rate, then sea ice temperature could decrease and restrain sea ice melting. However, the effect of the salinity difference term on the sea ice thickness and sea ice concentration was weak. To analyze the impacts of the salinity different term on Arctic sea ice thickness and sea ice concentration, we also designed several experiments by introducing the two salinity parameterizations to the ice-ocean coupled model, Modular Ocean Model （MOM4）, respectively. The simulated results confirmed the previous results of formula derivation.

HIGH RESOLUTION POSITIVITY-PRESERVING DIFFERENCE SCHEMES FOR TWO DIMENSIONAL EULER EQUATIONS

A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.

Parallel EOI algorithm with different insertion schemes

A parallel embedding overlapped iterative (EOI) algorithm about classicimplicit equations with asymmetric Saul'yev schemes (CIS-EOI) to solve one-dimensional diffusionequations is discussed to improve the properties of the segment classic implicit iterative (SCII)algorithm. The structure of CIS-EOI method is given and the stability of scheme and convergence ofiteration are proved by matrix method. The property of gradual-approach convergence is alsodiscussed. It has been shown that the convergent rate is faster and the property of gradual-approachconvergence also becomes better with the increasing of the net point in subsystems than with theSCII algorithm. The simulation examples show that the parallel iterative algorithm with a differentinsertion scheme CIS-EOI is more effective.

The Defining Schemes of Semantic Prosody

Sinclair is regarded as the father of semantic prosody study,however,the term "semantic prosody" was proposed by another famous linguist——Louw in 1993.From then on,several linguists proposed the definitions for this term and thus the features of semantic prosody become progressively clear.

ON DISPERSION-CONTROLLED PRINCIPLES FOR NON-OSCILLATORY SHOCK-CAPTURING SCHEMES

The role of dispersions in the numerical solutions of hydrodynamic equation systems has been realized for long time.It is only during the last two decades that extensive studies on the dispersion-controlled dissipative(DCD)schemes were reported.The studies have demonstrated that this kind of the schemes is distinct from conventional dissipation-based schemes in which the dispersion term of the modified equation is not considered in scheme construction to avoid nonphysical oscillation occurring in shock wave simulations.The principle of the dispersion controlled aims at removing nonphysical oscillations by making use of dispersion characteristics instead of adding artificial viscosity to dissipate the oscillation as the conventional schemes do.Research progresses on the dispersion- controlled principles are reviewed in this paper,including the exploration of the role of dispersions in numerical simulations,the development of the dispersion-controlled principles,efforts devoted to high-order dispersion-controlled dissipative schemes,the extension to both the finite volume and the finite element methods,scheme verification and solution validation,and comments on several aspects of the schemes from author's viewpoint.

SOME WEIGHT-TYPE HIGH-RESOLUTION DIFFERENCE SCHEMES AND THEIR APPLICATIONS

By the aid of an idea of the weighted ENO schemes, some weight-type high-resolution difference schemes with different orders of accuracy are presented in this paper by using suitable weights instead of the minmod functions appearing in various TVD schemes. Numerical comparisons between the weighted schemes and the non-weighted schemes have been done for scalar equation, one-dimensional Euler equations, two-dimensional Navier-Stokes equations and parabolized Navier-Stokes equations.

Water Vapor and Cloud Radiative Forcings over the Pacific Ocean Simulated by the LASG/IAP AGCM:Sensitivity to Convection Schemes

Characteristics of the total clear-sky greenhouse effect （GA） and cloud radiative forcings （CRFs）, along with the radiative-related water vapor and cloud properties simulated by the Spectral Atmospheric Model developed by LASGIAP （SAMIL） are evaluated. Impacts of the convection scheme on the simulation of CRFs are discussed by using two AMIP （Atmospheric Model Inter-comparison Project） type simulations employing different convection schemes： the new Zhang-McFarlane （NZH） and Tiedtke （TDK） convection schemes. It shows that both the climatological GA and its response to El Nio warming are simulated well, both in terms of spatial pattern and magnitude. The impact of the convection scheme on GA is not significant. The climatological longwave CRF （LWCRF） and its response to El Nio warming are simulated well, but with a prominently weaker magnitude. The simulation of the climatology （response） of LWCRF in the NZH （TDK） run is slightly more realistic than in the TDK （NZH） simulation, indicating significant impacts of the convection scheme. The shortwave CRF （SWCRF） shows large biases in both spatial pattern and magnitude, and the results from the TDK run are better than those from the NZH run. A spuriously excessive negative climatological SWCRF over the southeastern Pacific and an insufficient response of SWCRF to El Nio warming over the tropical Pacific are seen in the NZH run. These two biases are alleviated in the TDK run, since it produces vigorous convection, which is related to the low threshold for convection to take place. Also, impacts of the convection scheme on the cloud profile are discussed.

Stability of Difference Schemes with Intrinsic Parallelism for Quasilinear Parabolic Systems

In this paper, the first boundary problem of quasilinear parabolic system of second order is studied by the finite difference method with intrinsic parallelism. for the problem, the stability of the difference schemes with intrinsic parallelism are justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete data of the original problem, without assuming the existence of the smooth solutions for the origillal problem.

SCHEMES OF GAS PRODUCTION FROM NATURAL GAS HYDRATES

Natural gas hydrates are a kind of nonpolluting and high quality energy resources for future, the reserves of which are about twice of the carbon of the current fossil energy (petroleum, natural gas and coal) on the earth. And it will be the most important energy for the 21st century. The energy balance and numerical simulation are applied to study the schemes of the natural gas hydrates production in this paper,and it is considered that both depressurization and thermal stimulation are effective methods for exploiting natural gas hydrates, and that the gas production of the thermal stimulation is higher than that of the depressurization. But thermal stimulation is non-economic because it requires large amounts of energy.Therefore the combination of the two methods is a preferable method for the current development of the natural gas hydrates. The main factors which influence the production of natural gas hydrates are: the temperature of injected water, the injection rate, the initial saturation of the hydrates and the initial temperature of the reservoir which is the most important factor.