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.

Revisting High-Resolution Schemes with van Albada Slope Limiter

Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisit the question of non-oscillatory behavior of high-resolution central schemes in terms of the slope limiter proposed by van Albada et al.(Astron Astrophys 108:76–84,1982).The van Albada(vA)limiter is smoother near extrema,and consequently,in many cases,it outperforms the results obtained using the standard minmod limiter.In particular,we prove that the vA limiter ensures the one-dimensional Total-Variation Diminishing(TVD)stability and demonstrate that it yields noticeable improvement in computation of one-and two-dimensional systems.

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.

A New Class of Simple,General and Efficient Finite Volume Schemes for Overdetermined Thermodynamically Compatible Hyperbolic Systems

In this paper,a new efficient,and at the same time,very simple and general class of thermodynamically compatiblefinite volume schemes is introduced for the discretization of nonlinear,overdetermined,and thermodynamically compatiblefirst-order hyperbolic systems.By construction,the proposed semi-discrete method satisfies an entropy inequality and is nonlinearly stable in the energy norm.A very peculiar feature of our approach is that entropy is discretized directly,while total energy conservation is achieved as a mere consequence of the thermodynamically compatible discretization.The new schemes can be applied to a very general class of nonlinear systems of hyperbolic PDEs,including both,conservative and non-conservative products,as well as potentially stiff algebraic relaxation source terms,provided that the underlying system is overdetermined and therefore satisfies an additional extra conservation law,such as the conservation of total energy density.The proposed family offinite volume schemes is based on the seminal work of Abgrall[1],where for thefirst time a completely general methodology for the design of thermodynamically compatible numerical methods for overdetermined hyperbolic PDE was presented.We apply our new approach to three particular thermodynamically compatible systems:the equations of ideal magnetohydrodynamics(MHD)with thermodynamically compatible generalized Lagrangian multiplier(GLM)divergence cleaning,the unifiedfirst-order hyperbolic model of continuum mechanics proposed by Godunov,Peshkov,and Romenski(GPR model)and thefirst-order hyperbolic model for turbulent shallow waterflows of Gavrilyuk et al.In addition to formal mathematical proofs of the properties of our newfinite volume schemes,we also present a large set of numerical results in order to show their potential,efficiency,and practical applicability.

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.

A CLASS OF COMPACT UPWIND TVD DIFFERENCE SCHEMES

A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD （CTVD） schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.

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.

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.

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.

Fifth-Order A-WENO Schemes Based on the Adaptive Diffusion Central-Upwind Rankine-Hugoniot Fluxes

We construct new fifth-order alternative WENO(A-WENO)schemes for the Euler equations of gas dynamics.The new scheme is based on a new adaptive diffusion centralupwind Rankine-Hugoniot(CURH)numerical flux.The CURH numerical fluxes have been recently proposed in[Garg et al.J Comput Phys 428,2021]in the context of secondorder semi-discrete finite-volume methods.The proposed adaptive diffusion CURH flux contains a smaller amount of numerical dissipation compared with the adaptive diffusion central numerical flux,which was also developed with the help of the discrete RankineHugoniot conditions and used in the fifth-order A-WENO scheme recently introduced in[Wang et al.SIAM J Sci Comput 42,2020].As in that work,we here use the fifth-order characteristic-wise WENO-Z interpolations to evaluate the fifth-order point values required by the numerical fluxes.The resulting one-and two-dimensional schemes are tested on a number of numerical examples,which clearly demonstrate that the new schemes outperform the existing fifth-order A-WENO schemes without compromising the robustness.

Numerical Simulation of Modified Kortweg-de Vries Equation by Linearized Implicit Schemes

In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.

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.

The Vertical Structures of Atmospheric Temperature Anomalies Associated with El Nio Simulated by the LASG/IAPAGCM: Sensitivity to Convection Schemes

The vertical structures of atmospheric temperature anomalies associated with El Nio are simulated with a spectrum atmospheric general circulation model developed by LASG/IAP （SAMIL）. Sensitivity of the model’s response to convection scheme is discussed. Two convection schemes, i.e., the revised Zhang and Macfarlane （RZM） and Tiedtke （TDK） convection schemes, are employed in two sets of AMIP-type （Atmospheric Model Intercomparison Project） SAMIL simulations, respectively. Despite some deficiencies in the upper troposphere, the canonical El Nio-related temperature anomalies characterized by a prevailing warming throughout the tropical troposphere are well reproduced in both simulations. The performance of the model in reproducing temperature anomalies in ＂atypical＂ El Nio events is sensitive to the convection scheme. When employing the RZM scheme, the warming center over the central-eastern tropical Pacific and the strong cooling in the western tropical Pacific at sea surface level are underestimated. The quadru-pole temperature anomalies in the middle and upper troposphere are also obscured. The result of employing the TDK scheme resembles the reanalysis and hence shows a better performance. The simulated largescale circulations associated with atypical El Nio events are also sensitive to the convection schemes. When employing the RZM scheme, SAMIL failed in capturing the classical Southern Oscillation pattern. In accordance with the unrealistic anomalous Walker circulation and the upper tropospheric zonal wind changes, the deficiencies of the precipitation simulation are also evident. These results demonstrate the importance of convection schemes in simulating the vertical structure of atmospheric temperature anomalies associated with El Nio and should serve as a useful reference for future improvement of SAMIL.

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

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.

NUMERICAL STUDIES OF PLANETARY BOUNDARY LAYER PARAMETERIZATION SCHEMES ON SUPER TYPHOON SANBA(2012)DURING ITS INITIAL STAGE

The formation and development of typhoons are closely related to the disturbed low vortexes at the planetary boundary layer(PBL). The effects of five PBL parameterization schemes(PBL schemes hereinafter) on the trajectory,intensity, and distribution of physical quantities are studied using the mesoscale WRF model on Super Typhoon Sanba(2012) during its initial stage. Results show that the five PBL schemes exhibit significant different effects on the simulated intensity and path. The results simulated by QNSE and ACM2 without the Bogus method are close to the best track data in the numerical experiments. When the Bogus method is adopted, the simulated trajectories improve significantly because the initial field is close to the true data. Among the five PBL schemes, QNSE and ACM2 with the Bogus method present improved simulated path and intensity compared with the three other schemes. This finding indicates that the two schemes deal with the initial PBL process satisfactorily, especially in the formation and development of disturbed low vortexes. The differences in the treatment methods of the five PBL schemes affect the surface layer physical quantities and the middle and upper atmospheres during the middle to late periods of the typhoon.Although QNSE and ACM2 present better simulation results than other schemes, they exhibit a few differences in the internal structure of the typhoon. The results simulated by MYJ are worse, and this method may be unsuitable for studying the formation and development of typhoons.

Optimization of dewatering schemes for a deep foundation pit near the Yangtze River, China

A deep foundation pit constructed for an underground transportation hub was excavated near the Yangtze River. Among the strata, there are two confined aquifers, between which lies an aquiclude that is partially missing. To guarantee the safety of pit excavation, the piezometric head of the upper confined aquifer, where the pit bottom is located, should be 1 m below the pit bottom, while that of the lower confined aquifer should be dewatered down to a safe water level to avoid uplift problem. The Yangtze River levee is notably close to the pit, and its deformation caused by dewatering should be controlled. A pumping test was performed to obtain the hydraulic conductivity of the upper confined aquifer. The average value of the hydraulic conductivity obtained from analytical calculation is 20.45 m/d, which is larger than the values from numerical simulation（horizontal hydraulic conductivity K_H = 16 m/d and vertical hydraulic conductivity K_V = S m/d）. The difference between K_H and K_V indicates the anisotropy of the aquifer. Two dewatering schemes were designed for the construction and simulated by the numerical models for comparison purposes. The results show that though the first scheme could meet the dewatering requirements, the largest accumulated settlement and differential settlement would be94.64 mm and 3.3‰, respectively, greatly exceeding the limited values. Meanwhile, the second scheme,in which the bottoms of the waterproof curtains in ramp B and the river side of ramp A are installed at a deeper elevation of-28 m above sea level, and 27 recharge wells are set along the levee, can control the deformation of the levee significantly.