We construct a nonlinear monotone finite volume scheme for three-dimensional diffusion equation on tetrahedral meshes. Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme, we present a new efficient eliminating method. The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously. The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes, and also show that our scheme appears to be approximate second-order accuracy for solution.

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.

Currently large sky area spectral surveys like SDSS, 2dF, and LAMOST, using the new generation of telescopes and observatories, have provided massive spectral data sets for astronomical research. Most of the data can be automatically handled with pipelines, but visually inspection by human eyes is still necessary in several situations, like low SNR spectra, QSO recognition and peculiar spectra mining. Using ASERA, A Spectrum Eye Recognition Assistant, we can set up a team spectral inspection platform. On a preselected spectral data set, members of a team can individually view spectra one by one, find the best match template and estimate the redshift. Results from different members will be gathered and merged to raise the team work efficiency. ASERA mainly targets the spectra of SDSS and LAMOST fits data formats. Other formats can be supported with some conversion. Spectral templates from SDSS and LAMOST pipelines are embedded and users can easily add their own templates. Convenient cross identification interfaces with SDSS, SIMBAD, VIZIER, NED and DSS are also provided. An application example targeting finding strong emission line spectra from LAMOST DR2 is presented.

The extension of diamond scheme for diffusion equation to three dimensions is presented. The discrete normal flux is constructed by a linear combination of the directional flux along the line connecting cell-centers and the tangent flux along the cell-faces. In addition, it treats material discontinuities by a new iterative method. The stability and first-order convergence of the method is proved on distorted meshes. The numerical results illustrate that the method appears to be approximate second-order accuracy for solution.

In this paper, we are concerned with the constrained finite element method based on domain decomposition satisfying the discrete maximum principle for diffusion problems with discontinuous coefficients on distorted meshes. The basic idea of domain decomposition methods is used to deal with the discontinuous coefficients. To get the information on the interface, we generalize the traditional Neumann-Neumann method to the discontinuous diffusion tensors case. Then, the constrained finite element method is used in each subdomain. Comparing with the method of using the constrained finite element method on the global domain, the numerical experiments show that not only the convergence order is improved, but also the nonlinear iteration time is reduced remarkably in our method.

In this paper, we construct a global repair technique for the finite element scheme of anisotropic diffusion equations to enforce the repaired solutions satisfying the discrete maximum principle. It is an extension of the existing local repair technique. Both of the repair techniques preserve the total energy and are easy to be implemented. The numerical experiments show that these repair techniques do not destroy the accuracy of the finite element scheme, and the computational cost of the global repair technique is cheaper than the local repair technique when the diffusion tensors are highly anisotropic.

The carbuncle phenomenon has been regarded as a spurious solution produced by most of contact-preserving methods. The hybrid method of combining high resolution flux with more dissipative solver is an attractive attempt to cure this kind of non-physical phenomenon. In this paper, a matrix-based stability analysis for 2-D Euler equations is performed to explore the cause of instability of numerical schemes. By combining the Roe with HLL flux in different directions and different flux components, we give an interesting explanation to the linear numerical instability. Based on such analysis, some hybrid schemes are compared to illustrate different mechanisms in controlling shock instability. Numerical experiments are presented to verify our analysis results. The conclusion is that the scheme of restricting directly instability source is more stable than other hybrid schemes.

For a new nonlinear iterative method named as Picard-Newton (P-N) iterative method for the solution of the time-dependent reaction-diffusion systems, which arise in non-equilibrium radiation diffusion applications, two time step control methods are investigated and a study of temporal accuracy of a first order time integration is presented. The non-equilibrium radiation diffusion problems with flux limiter are considered, which appends pesky complexity and nonlinearity to the diffusion coefficient. Numerical results are presented to demonstrate that compared with Picard method, for a desired accuracy, significant increase in solution efficiency can be obtained by Picard-Newton method with the suitable time step size selection.

In this paper, the existence and uniqueness of the global smooth solution are proved for an evolutionary Ginzburg–Landau model for superconductivity under the Coulomb and Lorentz gauge.

In this paper the initial value problem for a class of Zakharov equations arising from ion-acoustic modes is discussed. Without assuming the Cauchy data are small, we prove the existence and uniqueness of the global smooth solution for the problem via the so-called continuous method and delicate a priori estimates.