Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

Dynamic release process of pollutants during suspended sediment transport in aquatic system

Pollutants release is highly consistent with suspended sediment concentration （SSC） in water column, especially during re-suspension and transport events. The present research focuses on pollutant dynamic release from re-suspended sediment, especially the vertical distribution relationship between them. The sediment erosion experiments on a series of uniform flow are conducted in a circulate flume. Reactive tracer （phosphorus） is used as the contaminant in fine-grained sediments to identify the release characteristic length and time. Experimental results show that the flow condition near-bed depends on the sediment surface roughness. The region with high turbulent intensities corresponds to a high concentration sediment layer. In addition, the SSC decreases with the distance, water depth, and particle grain size. The sediment in a smaller grain size takes much more time to reach equilibrium concentration. Total phosphorus （TP） concentration changes along the water depth as SSC in the initial re-suspension stage, appearing in two obvious concentration regimes： the upper low-concentration layer and the high-concentration near-bottom layer. This layered phenomenon remains for about 3 hours until SSC distri- bution tends to be uniform. Longitudinal desorption plays an important role in long-way transport to reduce the amount of suspended sediment in water column.     

Likelihood Inference under Generalized Hybrid Censoring Scheme with Comp eting Risks

Statistical inference is developed for the analysis of generalized type-II hybrid censoring data under exponential competing risks model. In order to solve the problem that approximate methods make unsatisfactory performances in the case of small sample size, we establish the exact conditional distributions of estimators for parameters by conditional moment generating function(CMGF). Furthermore, confidence intervals(CIs) are constructed by exact distributions, approximate distributions as well as bootstrap method respectively, and their performances are evaluated by Monte Carlo simulations. And finally, a real data set is analyzed to illustrate all the methods developed here.     

A Finite Volume Implicit Euler Scheme for the Linearized Shallow Water Equations: Stability and Convergence

The linearized shallow water equations are discretized in space by a finite volume method and in time by an implicit Euler scheme. Stability and convergence of the scheme are proved.     

Numerical simulation and convergence analysis for a system of nonlinear singularly perturbed differential equations arising in population dynamics

In this article, we consider a system of nonlinear singularly perturbed differential equations with two different parameters. To solve this system, we develop a weighted monotone hybrid scheme on a nonuniform mesh. The proposed scheme is a combination of the midpoint scheme and the upwind scheme involving the weight parameters. The weight parameters enable the method to switch automatically from the midpoint scheme to the upwind scheme as the nodal points start moving from the inner region to the outer region. The nonuniform mesh in particular the adaptive grid is constructed using the idea of equidistributing a positive monitor function involving the solution gradient. The method is shown to be second order convergent with respect to the small parameters. Numerical experiments are presented to show the robustness of the proposed scheme and indicate that the estimate is optimal.     

A new coding scheme of QC-LDPC codes for optical transmission systems

Based on Galois Field (GF(q)) multiplicative group, a new coding scheme for Quasi-Cyclic Low-Density Parity-Check (QC-LDPC) codes is proposed, and the new coding scheme has some advantages such as the simpler construction, the easier implementation encoding, the lower complexity of the encoding and decoding, the more flexible adjustment of the code length as well as the code rate and so forth. Under the condition of considering the characteristics of optical transmission systems, an irregular QC-LDPC (3843,3603) code to be suitable for optical transmission systems is constructed by applying the proposed new coding scheme. The simulation result shows that the net coding gain (NCG) of the irregular QC-LDPC (3843,3603) code is respectively improved 2.14 dB, 1.19 dB, 0.24 dB and 0.14 dB more than those of the classic RS (255,239) code in ITU-T G.975, the LDPC (32640,30592) code in ITU-T G.975.1, the regular SCG-LDPC (3969,3720) code constructed by the Systematically Constructed Gallager (SCG) coding scheme and the regular QC-LDPC (4221,3956) code at the bit error rate (BER) of 10-8. Furthermore, all the five codes have the same code rate of 93.7%. Therefore, the irregular QC-LDPC (3843,3603) code constructed by the proposed new coding scheme has the more excellent error-correction performance and can be better suitable for optical transmission systems.     

A new level‐dependent coarse grid correction scheme for indefinite Helmholtz problems

In this paper, we construct and analyze a level‐dependent coarse grid correction scheme for indefinite Helmholtz problems. This adapted multigrid (MG) method is capable of solving the Helmholtz equation on the finest grid using a series of MG cycles with a grid‐dependent complex shift, leading to a stable correction scheme on all levels. It is rigorously shown that the adaptation of the complex shift throughout the MG cycle maintains the functionality of the two‐grid correction scheme, as no smooth modes are amplified in or added to the error. In addition, a sufficiently smoothing relaxation scheme should be applied to ensure damping of the oscillatory error components. Numerical experiments on various benchmark problems show the method to be competitive with or even outperform the current state‐of‐the‐art MG‐preconditioned Krylov methods, for example, complex shifted Laplacian preconditioned flexible GMRES. Copyright © 2013 John Wiley & Sons, Ltd.     

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

A class of finite difference schemes for interface problems with an HOC approach

In this paper, we propose a new methodology for numerically solving elliptic and parabolic equations with discontinuous coefficients and singular source terms. This new scheme is obtained by clubbing a recently developed higher‐order compact methodology with special interface treatment for the points just next to the points of discontinuity. The overall order of accuracy of the scheme is at least second. We first formulate the scheme for one‐dimensional (1D) problems, and then extend it directly to two‐dimensional (2D) problems in polar coordinates. In the process, we also perform convergence and related analysis for both the cases. Finally, we show a new direction of implementing the methodology to 2D problems in cartesian coordinates. We then conduct numerous numerical studies on a number of problems, both for 1D and 2D cases, including the flow past circular cylinder governed by the incompressible Navier–Stokes equations. We compare our results with existing numerical and experimental results. In all the cases, our formulation is found to produce better results on coarser grids. For the circular cylinder problem, the scheme used is seen to capture all the flow characteristics including the famous von Kármán vortex street. Copyright © 2016 John Wiley & Sons, Ltd.