In the paper, a new time second-order mass-conserved implicit/explicit domain decomposition method (DDM) for the diffusion equations is proposed. In the scheme, firstly, we calculate the interface fluxes of sub-domains from the obtained solutions and fluxes at the previous time level, for which we apply high-order Taylor’s expansion and transfer the time derivatives to spatial derivatives to improve the accuracy. Secondly, the interior solutions and fluxes in sub-domains are computed by the implicit scheme and using the relations between solutions and fluxes, without any correction step. The mass conservation is proved and the convergence order of the numerical solutions is proved to be second-order in both time and space steps. The super-convergence of numerical fluxes is also proved to be second-order in both time and space steps. The scheme is stable under the stable condition r≤3/5. The important feature is that the proposed domain decomposition scheme is mass-conserved and is of second order convergence in time. Numerical experiments confirm the theoretical results.

In this paper, a conservative parallel iteration scheme is constructed to solve nonlinear diffusion equations on unstructured polygonal meshes. The design is based on two main ingredients: the first is that the parallelized domain decomposition is embedded into the nonlinear iteration; the second is that prediction and correction steps are applied at subdomain interfaces in the parallelized domain decomposition method. A new prediction approach is proposed to obtain an efficient conservative parallel finite volume scheme. The numerical experiments show that our parallel scheme is second-order accurate, unconditionally stable, conservative and has linear parallel speed-up.

In this paper, the high-order space-time discontinuous Galerkin cell vertex scheme (DG-CVS) developed by the authors for hyperbolic conservation laws is extended for time dependent diffusion equations. In the extension, the treatment of the diffusive flux is exactly the same as that for the advective flux. Thanks to the Riemann-solver-free and reconstruction-free features of DG-CVS, both the advective flux and the diffusive flux are evaluated using continuous information across the cell interface. As a result, the resulting formulation with diffusive fluxes present is still consistent and does not need any extra ad hoc techniques to cure the common “variational crime” problem when traditional DG methods are applied to diffusion problems. For this reason, DG-CVS is conceptually simpler than other existing DG-typed methods.

This note points out an error made in an earlier paper treating a diffusion equation for the probability distribution of gene frequencies under selection. The main interest lies in the orthogonality properties of certain eigenfunctions and the determination of the corresponding Fourier coefficients. Two particular cases of selection parameter values are considered. Approximate formulas are given for the two smallest eigenvalues for arbitrary values of the selection parameters.

A technique for solving a partial differential equation is presented. The technique is based upon the known solution of a similar equation. The method is used to attempt to solve the equation describing the change in the frequency of an allele in the presence of selection and random drift in a finite population. Two cases can be solved within a reasonable degree of approximation: (a) the viabilities are additive and (b) heterozygotes are symmetrically overdominant to the homozygotes. The solutions in both cases are compared with the exact discrete solutions found by powering the re evant transition matrix.