Insight into the accuracy of steady shock-capturing CFD methods is obtained through analysis of a simple problem involving steady transonic ow in a quasi-1D diverging duct. It is proved that the discrete solution error on either side of the shock is O(h n) where n is the order of accuracy of the conservative nite volume discretisation. Furthermore, it is shown that provided that n 2 then the er...

This paper focuses on three interconnected topics: (1) quantitative relationship between surface shortwave cloud radiative forcing, cloud fraction, and cloud albedo; (2) surface-based approach for measuring cloud albedo; (3) multiscale (diurnal, annual and inter-annual) variations and covariations of surface shortwave cloud radiative forcing, cloud fraction, and cloud albedo. An analytical expr...

We describe multi-material (more than two materials) interface reconstruction methods for 3D mesh of generalized polyhedrons. Basic information used in interface reconstruction is volume fractions of each material in mixed cell containing several materials. All methods subdivide mixed cell into set of pure non-overlapping sub-cells each containing just one material that have reference volume fr...

High-accuracy finite-difference schemes are used to solve the two-dimensional timedomain Maxwell equations for electromagnetic wave propagation and scattering. The high-accuracy schemes consist of a seven-point spatial operator coupled with a six-stage Runge–Kutta time-marching method. Two methods are studied, one of which produces the maximum order of accuracy and one of which is optimized for...

The hierarchical reconstruction (HR) [9, 20] is applied to a piecewise quadratic spectral volume method on two-dimensional unstructured grids as a limiting procedure to prevent spurious oscillations in numerical solutions. The key features of this HR are that the reconstruction on each control volume only uses adjacent control volumes, which forms a compact stencil set, and there is no truncati...

A postprocessing technique based on negative order norm estimates for the discontinuous Galerkin methods was previously introduced by Cockburn, Luskin, Shu, and Süli [Proceedings of the International Symposium on Discontinuous Galerkin Methods, Springer, New York, pp. 291– 300; Math. Comput., 72 (2003), pp. 577–606]. The postprocessor allows improvement in accuracy of the discontinuous Galerkin...

In many applications of ageand size-structured population models, there is an interest is obtaining good approximations of total population numbers rather than of their densities. Therefore, it is reasonable in such cases to solve numerically not the PDE model equations themselves, but rather their integral equivalents. For this purpose quadrature formulae are used in place of the integrals. Si...

In this paper, we further analyze, improve and develop the conservative fronttracking method developed in Mao(1995, SIAM. J. Numer. Anal., 32, 1677– 1703). The method in Mao(1995, SIAM. J. Numer. Anal., 32, 1677–1703) is only first-order accurate near the tracked discontinuities when applied to systems of conservation laws. This is because in the system case near tracked discontinuities there a...

We present a scheme for solving two-dimensional, nonlinear reaction-diiusion equations, s @p @t ? r (Krp) = f (p); using a mixed nite-element method. To linearize the mixed-method equations, we use a two grid scheme that relegates all of the Newton-like iterations to a grid 4H much coarser than the original one 4h, with no loss in order of accuracy so long as the mesh sizes obey H = O(p h). The...

In this paper, an efficient and unconditionally stable rectangular domain decomposition algorithm is proposed. The order of accuracy of the prediction scheme of the new algorithm is second. Numerical experiments support efficiency, accuracy, and unconditional stability of the method.