In this work, we study the unsteady axisymmetric mixed convection boundary layer flow and heat transfer of non-Newtonian power law fluid over a cylinder. Different from most classical works, the temperature dependent variable fluid viscosity and thermal conductivity are taken into account in highly coupled velocity and temperature fields. The motion of the fluid can be modeled by a time-dependent nonlinear parabolic system in cylindrical coordinates, which is solved numerically by using Chebyshev spectral method along with the strong stability preserving (SSP) third order Runge-Kutta time discretization. We apply the numerical solver to problems with different power law indices, viscosity parameter, thermal conductivity parameter and Richardson numbers, and compute up to the steady state. The numerical solver is checked by testing the spectral convergence of the numerical approximation to a smooth exact solution of the PDEs with source terms. Moreover, the combined effects of pertinent physical parameters on the flow and heat transfer characteristics are analyzed in detail.

The numerical simulation of the band structure of three-dimensional dispersive metallic photonic crystals with face-centered cubic lattices leads to large-scale nonlinear eigenvalue problems, which are very challenging due to a high dimensional subspace associated with the eigenvalue zero and the fact that the desired eigenvalues (with smallest real part) cluster and close to the zero eigenvalues. For the solution of the nonlinear eigenvalue problem, a Newton-type iterative method is proposed and the nullspace-free method is applied to exclude the zero eigenvalues from the associated generalized eigenvalue problem. To find the successive eigenvalue/eigenvector pairs, we propose a new non-equivalence deflation method to transform converged eigenvalues to infinity, while all other eigenvalues remain unchanged. The deflated problem is then solved by the same Newton-type method, which is used as a hybrid method that combines with the Jacobi-Davidson and the nonlinear Arnoldi methods to compute the clustering eigenvalues. Numerical results illustrate that the proposed method is robust even for the case of computing many clustering eigenvalues in very large problems.

No abstract uploaded!

In this paper we consider a discontinuous Galerkin discretization of the ideal magnetohydrodynamics (MHD) equations on unstructured meshes, and the divergence free constraint ($\nabla \cdot \B = 0$) of its magnetic field $\B$. We first present two approaches for maintaining the divergence free constraint, namely the approach of a locally divergence free projection inspired by locally divergence free elements \cite{Li2005}, and another approach of the divergence cleaning technique given by Dedner et al. \cite{Dedner2002}. By combining these two approaches we obtain a efficient method at the almost same numerical cost. Finally, numerical experiments are performed to show the capacity and efficiency of the scheme.

This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal magnetohydrodynamics (MHD) on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of numerical MHD schemes with a Harten–Lax–van Leer (HLL) type flux on polytopal meshes in any space dimension. The main challenges overcome here include establishing certain relation between the PP property and a discrete divergence of magnetic field on general meshes, and estimating proper wave speeds in the HLL flux to ensure the PP property. In the 1D case, we prove that the standard DG and finite volume methods with the proposed HLL flux are PP, under a condition accessible by a PP limiter. For the multidimensional conservative MHD system, the standard DG methods with a PP limiter are not PP in general, due to the effect of unavoidable divergence error in the magnetic field. We construct provably PP high-order DG and finite volume schemes by proper discretization of the symmetrizable MHD system, with two divergence-controlling techniques: the locally divergence-free elements and suitably discretized Godunov–Powell source term. The former technique leads to zero divergence within each cell, while the latter controls the divergence error across cell interfaces. Our analysis reveals in theory that a coupling of these two techniques is very important for positivity preservation, as they exactly contribute the discrete divergence terms which are absent in standard multidimensional DG schemes but crucial for ensuring the PP property. Several numerical tests further confirm the PP property and the effectiveness of the proposed PP schemes. Unlike the conservative MHD system, the exact smooth solutions of the symmetrizable MHD system are proved to retain the positivity even if the divergence-free condition is not satisfied. Our analysis and findings further the understanding, at both discrete and continuous levels, of the relation between the PP property and the divergence-free constraint.